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The  magnetohydrodynamic  (MHD)  and  thermoacoustic  sound  generation  mechanisms 
were  investigated  analytically  and  experimentally  for  a  source  volume  region 
generated  by  a  time  harmonic  electric  current  distribution  interacting  in  a 
static  magnetic  field  using  sodium  chloride  and  water  as  an  electrolyte. 

The  MHD  scund  source  is  the  result  of  the  Lorentz  force  on  a  conducting  medium 
generated  by  the  interaction  of  electrical  current  and  a  magnetic  field.  The 
thermoacoustic  source  is  the  result  of  fluid  density  changes  produced  by  Joule 
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heating  when  electrical  power  is  dissipated  in  a  resistive  medium.  ,  MHD  and 
thermoacoustic  source  mechanisms  were  found  to  produce  two  distinct  acoustic 
radiation  moments:  a  dipole  field  and  monopole  field,  respectively. 

The  source  volume  region  is  contained  in  a  rigid  walled  waveguide  where  sound 
from  the  source  volume  region  radiates  into  the  free  medium  through  apertures 
at  either  end  of  the  waveguide.  Transmitting  sensitivities  for  both  mechanisms 
were  derived  analytically  for  the  waveguide  geometry  and  compared  to  experimental 
data  where  good  agreement  was  found  in  a  frequency  band  of  500  to  4000  Hz.  The 
source  directivity  for  the  dipole  moment  radiation  was  measured  and  found  to  be 
in  good  agreement  with  the  classic  cosine  radiation  pattern  at  low  frequencies. 

Expressions  fur  the  transduction  process  efficiency  and  input  electrical  imped¬ 
ance  were  also  derived  and  compared  to  measurements. 

r\ 

/, 


for 

'^TIS  CRA&i 
1  CllC  TAB 
U’lj'irour.f'pfj 

By  . . 

B'lit-ibjtio"/  1 

'dilib-itly  -Aides  "  1 

An,  !,l  .  j  I  Qr 
'■t  Olvll 


i 


UNCLASSIFIED 


security  classification  of  this  RACE(IWi#n  Dmit  Enttfd) 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES  v 

LIST  OF  IMPORTANT  SYMBOLS  ix 

Chapter  1  INTRODUCTION  1 

A.  Past  Research  in  the  Area  of  MHD  Transdi'ction  3 

B.  Present  Work  4 

Chapter  2  THEORY;  THE  MAGNETOHYDRODYNAMIC- 

THERMOACOUSTIC  WAVE  EQUATION 

A.  Ohm's  Law  for  an  Electrolyte  7 

B  Derivation  of  the  Continuity  Equation  9 

C.  Derivation  of  the  Momentum  Equation  f  5 

D.  Small  Signal  Acoustic  Approximatioi  1 7 

E.  One-Dimensional  Wave  Equation  20 

F.  Dispersion  of  the  Acoustic  Field  due  to  the  Induced  Current 

Density  in  a  High  Conductivity  Fluid  20 

Chapter  3  THEORY:  SOLUTION  OF  THE  WAVE  EQUATION  FOR 

THE  ACOUSTIC  FIELD  WITHIN  THE  WAVEGUIDE  23 

A.  Boundary  Conditions  of  the  Transducer  Waveguide  23 

B.  Justification  of  the  Boundary  Conditions  29 

C.  Acoustic  Wavefield  due  to  the  Magnetohydrodynamic  Source  35 

D.  Acoustic  Wavefield  due  to  the  Thermoacoustic  Source  47 

Chapter  4  THEORY:  ELECTROACOUSTIC  RELATIONS  53 

A.  Polarization  Impedance  of  the  MHD  Waveguide  Electrodes  53 

B.  Electrical  Impedance  of  the  MHD  Transduction  Process  59 

C.  MHD  Transduction  Process  Power  Efficiency  64 


Page 

D.  MHD  Transmitting  Current  Sensitivity  66 

E.  Electrical  Impedance  of  the  Thermoacoustic  Transduction 

Process  73 

F.  Thermoacoustic  Transduction  Process  Power  Efficiency  78 

G.  Transmitting  Sensitivity  of  the  Thermoacoustic  Process  83 

Chapter  5  TRANSMITTING  EXPERIMENT  87 

A.  Optimization  of  the  MHD  Transmitting  Current  Sensitivity  87 

B.  Acoustic  Speed  Measurement  of  the  Waveguide  94 

C.  Measurement  of  the  Electrical  Impedance  of  the  Transducer  97 

D.  Transmitting  Circuit  and  Acoustic  Signal  Measurement 

System  101 

E.  Measurement  of  the  MHD  Transmitting  Sensitivity  and 

Directivity  1 07 

F.  Measurement  of  the  Thermoacoustic  Transmitting  Sensitivity  1 1 4 

Chapter  6  CONCLUSIONS  119 

Appendix  A.  1  DERIVATION  OF  THE  MAGNETIC  DISPERSION 

RELATION  123 

Appendix  A.2  THERMOVISCOUS  PROCESSES  WITHIN  THE 

WAVEGUIDE  126 

Appendix  -B  PRESSURE  FIELD  SPHERICAL  HARMONIC 

FUNCTION  EXPANSION,  RADIATION  IMPEDANCES, 

AND  FARFIELD  PRESSURE  1 30 

Appendix  C  REFLECTION  COEFFICIENT  EQUATIONS  1 34 

Appendix  D  COMPUTER  PROGRAMS  136 

REFERENCES  153 


iv 


LIST  OF  FIGURES 


Fi.gui:g  Pagg 

1.1  Dipole  Radiation  Pattern  MHD  Transducer  2 

2.1  Integral  Continuity  Equation  Control  Volume  1 1 

2.2  Integral  Momentum  Equation  Control  Volume  16 

2.3  MHD  Force  and  Thermoacoustic  Volume  Expansion 

Acting  in  a  Spherical  Volume  1 9 

3.1  Fluid  "Piston"  Acoustic  Aperture  Model  25 

3.2  MHD  Source  Geometry  and  Coordinate  Systems  27 

3.3  Self-Impedance  of  the  MHD  Waveguide  Apertures  30 

3.4  Mutual  Impedance  of  the  MHD  Waveguide  Apertures  31 

3.5  Ratio  of  Mutual  to  Self-Impedance  Magnitudes  33 

3.6  Ratio  of  Real  Components  of  the  Self-  and  Mutual 

Impedance  34 

3.7  Delta  Function  Source  Representation  for  the  Plane  Wave 

Green's  Function  Solution  37 

3.8  Magnitude  of  the  Reflection  Coefficient  42 


3.9  Magnitude  of  the  Modal  Particle  Velocity  and  Modal 
Pressure,  6%  NaCI-HgO  Solution 

ilgr:  0.10  m,  il  -  0.10  m,  Cq=  1536  m/c,  R  =  0.09  m  44 

3.10  Magnitude  of  the  Modal  Particle  Velocity  and  Modal 
Pressure,  6%  NaCI-HgO  Solution 
fi.g  =  0.10  m,  Jl  =  0.10  m,  Cq=  650  m/s,  R  =  0.09  m 


V 


45 


Borne  Ease 

3.1 1  Magnitude  of  Modal  Particle  Velocity  and  Modal 
Pressure,  6%  NaCI-H20  Solution 

ilj,  =  0.05  m,  il  =  0.10  m.  c^=  650  m/s,  R  =  0.09  m  46 

4.1  Equivalent  Circuit  of  the  MHD  Input  Impedance  60 

4.2  MHD  Mechanism  Impedance  for  6%  NaCI-H20  Solution  63 

4.3  MHD  Transduction  Process  Power  Efficiency  65 

4.4  Transmitting  Current  Sensitivity  Level  as  a  Function  of 

Sound  Speed  and  Frequency  69 

4.5  Transmitting  Current  Sensitivity  Level  as  a  Function  of 

Baffle  Radius  and  Frequency  70 

4.6  Transmitting  Current  Sensitivity  Level  as  a  Function  of 

Source  Length  and  Frequency  71 

4.7  Quadratic  Current  Transfer  Function, 

Single  Frequency  Input  for  6%  NaCI  Solution  77 

4.8  Quadratic  Current  Transfer  Function, 

Single  Frequency  Input  vor  Mercury  79 

4.9  Equivalent  Circuit  of  the  Combined  MHD  and 

Thermoacoustic  Input  Impedance  80 

4.10  Thermoacoustic  Transduction  Process  Efficiency  82 

4.1 1  Quadratic  Voltage  Transmitting  Sensitivity, 

Single  Frequency  Input  (QVTS)  85 

5.1  Magnetohydrodynamic  Transducer,  Aperture  View  88 


VI 


Eigurg 


Eagg. 


5.2 

5.3 

5.4 

5.5 

5.6 

5.7 

5.8 

5.9 

5.10 

5.11 

5.12 


Transmitting  Current  Sensitivity  as  a  Function  of 
Aperture  Radius,  for  =  1550  m/s,  R  =  0.05  m,  B  =  1  T, 
f  =  1000  Hz  90 

Transmitting  Current  Sensitivity  Level  (TSCL)  as  a  Function 
of  Aperture  Radius  and  Frequency,  for  c^  =  1 550  m/s, 

R  =  0.05m,  B  =  1T  91 

Magnetic  Induction  as  a  Function  of  the  Transverse 

Waveguide  Dimension,  (Field  Gap)  93 

Magnetohydrodynamic  Transducer  Waveguide  Construction 

Details  95 

Measured  Magnetic  Induction  Levels  (Gauss)  in  the  MHD 

Transducer  Waveguide  96 

Phase  Speed  Measurement  System  98 

Block  Diagram  of  Transducer  Electrical  Impedance 

Measurement  System  99 

Measured  Input  Electrical  Impedance  of  the  MHD 

Transducer  Containing  a  6%  NaCI-HgO  Solution  100 

Theoretical  and  Measured  Electrical  Input  Impedance 

of  the  MHD  Transducer  1 02 

Block  Diagram  of  MHD  Acoustic  Signal  Measurement 

System  1 04 

MHD  Transducer  Mounted  on  the  Rotator  Column  with 

the  H56  Hydrophone  106 


VII 


Figure  Bags 

5.13  Constructive  Interference  of  the  MHD  Pressure  Signal 

on  the  Waveguide  Axis,  \j/  =  0°,  at  6900  Hz  1 08 

5.14  Magnetohydrodynamic  Transmitting  Sensitivity 
Experiment  and  Theory,  6%  NaCI-HgO  Solution, 

Measurements  Made  on  the  Waveguide  x  Axis,  \j/  =  0°  109 

5.15  Magnetohydrodynamic  Source  Directivity, 

Farfield  Theory  and  Experiment  1 1 1 

5. 1 6  Effective  Baffle  Radii  of  the  MHD  T ransducer  1 1 3 

5.17  Second  Harmonic  Thermoacoustic  Pressure  Signal  115 

5. 1 8  Thermoacoustic  Quadratic  Voltage  T ransmitting 

Sensitivity  (QVTS),  Theoretical  and  Experimental  Data 
Comparison  116 

A. 2.1  Peak  Modal  Pressure  in  the  MHD  Waveguide  as  a  Function 

of  Frequency,  Q  Determination  Plot,  6%  NaCI-HgO  Solution  1 29 

B. 1  Aperture  Geometry  131 


LIST  OF  IMPORTANT  SYMBOLS 


Note  that  the  bold  symbols  imply  vector  fields. 


a,  b 

X  location  of  the  left  and  right 

MHD  source  volume  boundaries 

A 

transducer  electrical  terminal 
admittance 

B 

magnetic  induction  vector  field 

Cfm 

adiabatic  sound  speed  in  the  free 
medium 

Co 

adiabatic  sound  speed  in  the 
waveguide 

Cd 

diffusional  capacitance 

Co 

electric  double  layer  capacitance 

specific  heat  at  constant 
pressure 

D 

diffusion  constant,  1 .0  X  10'^  m^/s 

e.Qo 

relative  permittivity  and 
permittivity  of  free  space, 

8.85  XI 0-12  F/m 

■'I 

E 

electric  intensity  vector  field 

f  =  JXB 

magnetohydrodyamic  fluid 
volume  force 

F 

Faraday  constant, 

96524  coulomb/mole 

G(x,Xo) 

plane  wave  Green's  function  for 
the  waveguide  for  a  Dirac 
source  at  x^ 

GJx.a),  GR(x,b) 

plane  wave  Green's  function  valid 
in  the  waveguide  to  the  left  and 
right  of  the  MHD  source  volume, 
respectively 

GxL(x,a), 

derivatives  of  the  plane  wave 
Green's  function  with  respect  to 
the  X  coordinate 

h,Ji) 

Hankel  function  of  the  first  kind 
of  order  n 

H 

Heavyside's  unit  step  function 

I 

terminal  electric  current 
magnitude 

j 

J|  Jq 

current  density  vector  field  and 
the  step  function  amplitude 

k(m  --- 

wave  number  of  the  free  medium 

o 

O 

3 

II 

wave  number  of  the  medium  in 
the  waveguide 

-.d,  iloriil 

notation  for  indicating  the  x 
location  of  the  left  and  right 
aperture,  respectively 

ily,  ilz 

the  transverse  y  direction 
dimension  and  z  direction 
dimension  of  the  rectangular 
waveguide 

I 

path  integration  vector 

n 

number  of  moles  of  electrons 
required  per  mole  of  reaction 
product 

N 

volumetric  concentration, 
moles/m^ 

P 

waveguide  acoustic  pressure 
field 

P±lfm 

free  space  pressure  field  due  to 
the  aperture  at  -il  or  il 

P 

total  pressure  P  =  Pq  +  p 

X 


^n’  ^n-1’  ^n+1 

ordinary  Legendre  polynon^ial  of 
order  n,  n-1,  n+1 

Po 

ambient  pressure 

q 

unit  electronic 
charge,1.602X  10-'-C 

QCTF 

quadratic  current  transfer 
function  of  the  thcrmoacoustic 
source 

QVTS 

quadratic  voltage  transmitting 
sensitivity  of  the 
thermoacoustic  source 

r 

free  medium  radial  coordinate 

R 

ideal  gas  constant 

8.31  joules/(Kelvin  mole) 

Rg  =  V 

effective  aperture  radius  of  the 
rectangular  waveguide 

Rd 

real  component  of  the  diffusion 
impedance 

R±i 

aperture  reflection  coefficient 
evaluated  at  aperture  -il  or  il 

Rs 

spherical  baffle  radius 

s 

unit  area  normal  electrode 
surface  vector 

t 

time 

T 

absolute  temperature,  Kelvin 
scale 

TCS 

MHD  mechanism  transmitting 
current  sensitivity 

u 

waveguide  small  signal  particle 
velocity 

V 

volume 

terminal  potential 

radiated  acoustic  power 

input  electrical  po’’  .-rto 
transducer 

waveguide  longitudinal 
coordinate 

waveguide  transverse 
coordinates 

MHD  transducer  electrolyte 
impedance 

electrical  diffusion  impedance 

aperture  mutual  radiation 
impedance 

polarization  impedance  of  the 
electrodes 

aperture  self-radiation 
impedance 

transducer  electrical  terminal 
impedance 

coefficient  of  thermal  e).pansion, 
1/T 

Dirac  delta  function 

small  X  coordinate  increment 

transducer  power  efficiency 

aperture  half  angle 

phase  of  the  aperture  particle 
velocities 

number  of  moles  of  ions  required 
per  mole  of  reaction  product 


concentration  spatial  coordinate, 
direction  normal  to  electrode 
surface 

total  fluid  mass  density  p  =  po+6p 

acoustic  fluid  mass  density 
perturbation 

waveguide  fluid  mass  density 

free  medium  fluid  mass  density 

ionic  fluid  conductivity 

waveguide  acoustic  potential 
wavefield 

free  medium  angular  coordinate 
angular  frequency 
angular  frequency,  2co 


C!  lAPTER  1 


INTRODUCTION 

The  subject  matter  of  this  thesis  is  an  investigation  of  the 
underwater  acoustic  transmitting  characteristics  of  a  magnetohydrodynamic 
(MHD)  transducer.  The  transduction  mechanism  of  the  MHD  transducer  is 
the  Lorentz  force  generated  on  a  conducting  fluid  conveying  electrical 
charge  in  the  presence  of  a  magnetic  field.  The  transducer  consists  of  a 
sodium  chloride-v.'ater  solution  (electrolyte)  contained  in  a  rigid  wall 
rectangular  waveguide,  10  cm  long  and  3.8  cm  x  3.2  cm  cross  section,  with 
acoustic  apertures  at  either  end.  The  waveguide  is  located  between  the 
poles  of  two  large  permanent  magnets  (see  Fig  1.1). 

A  time  harmonic  electric  current  density,  Je‘i“*,  is  applied  to  the 
electrolyte  through  electrodes  on  two  facing  walls  which  are  parallel  to  the 
magnetic  field,  B.  The  current  density,  which  is  orthogonal  to  the  magnetic 
field  and  the  waveguide  axis,  produces  a  force  on  the  fluid,  f=JXBe‘i“',  in  the 
longitudinal  (axis)  direction  of  the  waveguide.  Assuming  the  J  and  B  fields  to 
be  uniform,  the  time  harmomic  force  generates  a  plane  wave  acoustic 
pressure  field  within  the  waveguide  which  radiates  into  the  free  medium  via 
the  apertures. 

The  electromechanical  analogy  of  this  transducer  is  the  moving 
coil  transducer  where  the  fluid  of  the  MHD  transducer  i.j  likened  to  the  coil  of 
the  armature.  But,  unlike  the  moving  coil  device  which  radiates  sound  by  the 
motion  of  a  surface,  the  MHD  transducer  has  no  radiating  surface.  The 
acoustic  wave  is  a  result  of  the  direct  application  of  the  electromagnetic 
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FIGURE  1.1 

DIPOLE  RADIATION  PATTERN  MHD  TRANSDUCER 


energy  to  the  medium  (more  appropriately,  the  radiating  region  of  the 
transducer  is  termed  the  MHD  source  volume). 

A.  PAST  RESEARCH  IN  THE  AREA  OF  MHD  TRANSDUCTION 

The  first  attempts  to  use  the  MHD  mechanism  in  fluids  for 

transmitting  or  receiving  acoustic  signals  are  not  readily  identifiable.  Various 
steady  fluid  flow  rate  measuring  devices  and  electromagnetic  pumps  have 
been  in  existence  since  at  least  1945."'  Acoustic  transmitting  investigations, 
in  the  context  that  they  will  be  discussed  in  this  thesis,  were  first  carried  out 
by  Campanella^  in  1955,  and  later  by  Ajisaka  and  Hixson^’"*  in  1975. 

The  investigations  of  Campanella  and  Ajisaka  and  Hixson 
primarily  addressed  both  theoretically  and  experimentally  the  generation  of 
plane  wave  signals  by  the  MHD  source  radiating  in  a  rigid  walled 
waveguide  of  semi-infinite  length  with  rigid  and  pressure  release  boundaries 
at  the  origin.  The  MHD  source  region  extented  from  the  boundary  to  a  finite 
distance,  ,-klso  the  infinite  length  waveguide  was  investigated  by  both 
investigators. 

Ajisaka  and  Hixson^  also  considered  the  finite  length  waveguide 
radiating  from  both  apertures  into  unbounded  half  spaces.  They  treated  the 
apertures  of  the  waveguide  as  pistons  on  an  infinite  baffle  and  derived  a 
pressure  transmitting  sensitivity  expression  which  is  flat  over  a  very  broad 
frequency  band.  However  the  sensitivity  expression  which  was  derived  was 
for  the  apenure  pressure  and  not  the  pressure  at  the  standard  1  m  distance. 
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B.  PRESENT  WORK 


This  thesis  is  a  continuation  of  the  work  begun  by  Ajisaka  and 
Hixson,  by  suggestion  of  Dr.  Elmer  Hixson,  with  the  intent  to  investigate  in  a 
more  precise  manner  the  radiating  characteristics  of  the  MHD  source  in  an 
infinite  free  space.  Since  Ajisaka  and  Hixson  did  not  experimentally 
investigate  the  free  medium  transmitting  characteristics  of  the  MHD  source, 
this  was  to  be  a  primary  task  of  this  investigation. 

Chapter  2  of  this  study  is  the  rigorous  derivation  of  the 
inhomogeneous  pressure  wave  equation  governing  the  application  of  the 
MHD  force  to  a  fluid  medium.  It  was  realized  from  considerations  of  the  first 
law  of  thermodynamics  that  when  electrical  current  is  passed  through  a 
medium  of  finite  conductivity  heat  is  generated.  The  generation  of  heat  in  the 
medium  results  in  the  thermal  expansion  of  the  fluid,  which  creates  an 
additional  source  of  sound.  This  thermal  source  mechanism  is  accounted  for 
in  the  derivation  of  the  inhomogeneous  wave  equation.  The  two  source 
terms  of  the  wave  equation  are  recognized  mathematically  as  being 
monopole  and  dipole  source  distributions,  corresponding  respectively  to  the 
thermal  and  MHD  source  mechanisms.  The  thermal  source  is  shown  to  be  a 
quadratically  nonlinear  source  in  the  frequency  domain;  the  spectral  nature 
is  discussed  in  Chapters  2  and  4. 

The  solution  of  the  wave  equation  for  the  acoustic  wavefield  due 
to  the  MHD  and  thermal  sources  in  the  waveguide  cannot  be  found  by  the 
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usual  convolution  of  the  free  space  Green's  function  with  the  source 
distribution  due  to  the  presence  of  the  waveguide  boundaries  and  magnets. 
Thus  finding  the  appropriate  Green’s  function  for  this  complicated  geometry 
was  not  attempted. 

In  Chapter  3  solutions  for  the  acoustic  fields  of  both  source 
mechanisms  are  found  by  separating  the  acoustic  fields  into  two  regions: 
inside  the  waveguide  and  the  exterior  free  medium  fields.  A  plane  wave 
solution  for  the  acoustic  field  is  proposed  for  the  waveguide  region  assuming 
impedance  boundary  conditions  at  the  apertures.  To  find  the  exterior 
acoustic  field  the  waveguide  and  magnet  assembly  were  assumed  to 
comprise  a  rigid  spherical  surface  with  two  radiating  apertures.  The 
spherical  surface  allows  for  the  exterior  pressure  field  to  be  solved 
analytically  in  terms  of  a  Hankel  function  expansion  of  a  spherical  wavefield 
given  the  velocity  distribution  of  the  apertures.  Using  the  Hankel  function 
expansion  of  the  exterior  wavefield,  self-  and  mutual  impedance  relations 
were  found  for  the  apertures,  thus  allowing  the  solution  of  the  interior  and 
therefore  the  exterior  acoustic  fields.  Geometrical  and  wave  number 
limitations  are  placed  on  the  domain  of  the  solutions  due  to  the  matching  of  a 
plane  wavefield  within  the  waveguide  to  the  spherical  field  on  the  exterior. 

In  Chapter  4  derivations  of  the  transmitting  electrical  impedance, 
efficiency,  and  sensitivity  expressions  are  provided.  The  electrochemistry 
and  the  electrical  impedance  nature  of  an  electrolyte  in  the  vicinity  of  an 
electrode  surface  are  discussed,  and  an  equivalent  circuit  representation  of 
the  these  processes  is  presented.  Included  in  an  appendix  is  a  computer 
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program  which  computes  the  wavefields  and  transducer  characteristics 
discussed  above. 

Chapter  5  discusses  the  design  of  the  transducer  used  in  the 
experiment  and  the  measurement  procedures,  and  compares  the  theoretical 
predictions  to  experimental  data.  As  will  be  shown,  the  analysis  is  in  good 
agreement  with  experimental  measurements  up  to  frequencies  at  which 
differences  in  assumed  and  actual  transducer  geometries  become  important. 

Chapter  6  presents  a  summary  of  the  research  and  conclusions 
which  may  be  drawn. 


6 


CHAPTER  2 


THEORY:  THE  MAGNETOHYDRODYNAMIC-THERMOACOUSTIC 

WAVE  EQUATION 

In  this  chapter  the  small  signal,  inhomogeneous  pressure  wave 
equation  governing  the  acoustic  wavefield  within  a  magnetohydrodynamic 
source  volume  is  derived.  The  MHD  acoustic  field  is  found  to  have  two 
source  mechanisms:  (1)  the  primary  machanism,  the  magnetohydrodynamic 
force  and  (2)  a  thermal  expansion  mechanism  due  to  joule  heating  of  the 
medium  by  the  electric  current.  The  second  source  is  known  as  a 
thermoacoustic  source. 

A.  OHM'S  LAW  FOR  AN  ELECTROLYTE 

In  the  derivation  of  the  continuity,  momentum,  and  wave 
equations  to  be  presented  in  this  chapter  it  is  assumed  that  Ohm's  law  is 
valid  for  the  conduction  of  electrical  charge  in  the  electrolyte  solution  of  the 
MHD  transducer.  The  generalized  Ohm's  law  is  stated  below  under  the 
additional  assumption  that  the  ohmic  conductivity  a  is  a  constant  (not 
dependent  on  frequency  or  space), 

J  =  a(E  +  uXB)  +  eepE,  (2.1) 

J,  E,  B,  and  u  in  Eq.  (2.1)  are  defined  as  the  current  density  vector  field 
(A/m^),  the  electric  vector  field  (V/m),  magnetic  induction  vector  field,  and  the 
electrolyte  fluid  particle  velocity  field,  e  and  e^  are,  respectivily,  the  relative 
and  free  space  permitivities. 

In  this  thesis  it  is  justifiable  to  neglect  the  effects  of  the  induced 
current  (au  x  B)  term  and  displacement  current  (ee^E,)  term  in  Eq.  (2.1).  The 
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grounds  on  which  these  approximations  are  valid  are  as  follows. 

(1)  The  displacement  current  term  can  be  rewritten  in  the  form 
-jcoeegE  for  a  time  harmonic  signal  e‘K  The  maximum  frequency  of  interest 
in  (his  investigation  is  approximately  15000  Hz,  and  the  relative  dielectric 
permitivity  of  water  is  88.  Therefore  the  coefficient  of  the  displacement 
current,  coeeg,  could  be  as  great  as  7.34  X  10'^  1/ohm  m.  The  measured 
value  of  the  ohmic  conductivity  of  a  6%  NaCI  solution  is  5.28  1/ohm  m. 
Therefore,  on  the  basis  of  ihe  dielectic  conductivity  being  very  small  relative 
to  the  ohmic  conductivity,  the  displacement  term  is  justifiably  neglected. 

(2)  The  small  signal  acoustic  approximation  will  be  nr.ade  in  the 
next  section  in  the  derivation  of  the  acoustic  wave  equation.  This  implies  that 
the  particle  velocity  u  must  be  much  less  in  magnitude  than  the  adiabatic 
sound  speed  c^.  For  the  6%  sodium  chloride  solution  im'estigated  in  the 
experiment  the  value  of  c^  is  1550  m/s.  The  amplitude  of  the  induction  field  B 
considered  is  of  the  order  of  1  tesla.  If  the  maximum  particle  velocity  is  150 
m/s,  then  the  induced  current  density,  auXB,  will  be  792  A/m^.  The  typical 
current  density  used  in  the  experiment  is  15000  A/m-,  so  the  projected 
maximum  induced  current  density  is  only  5%  of  the  typical  applied  current 
density.  Therefore  induced  current  is  justifiably  neglected,  but  as  a  practical 
matter  the  particle  velocities  generated  in  the  experiment  were  of  the  order 
lO’'*  m/s  so  the  approximation  is  very  good. 

With  the  above  approximations  applied  to  Eq.  (2.1),  Ohm's  law  is 

restated  as 

J  =  oE  .  (2.2) 


8 


More  will  be  said  concerning  the  electrolyte  conductivity  in  Chapter  4  in 
connection  with  the  input  impedance  of  the  transducer. 

B.  DERIVATION  OF  THE  CONTINUITY  EOUATION 

Temperature  rises  associated  with  a  time  varying  electrical 
current  flowing  in  a  resistive  fluid  produce  a  monopole  acoustic  source 
distribution  due  to  thermal  volume  expansions  of  the  fluid.  The 
hydrodynamic  equations  for  what  is  termed  in  this  thesis  an  ohmic 
thermoacoustic  source  have  not  been  derived  explicitly  in  any  of  the  popular 
literature  on  acoustics.  Didenkulov^  treats  the  problem  of  current  flowing  in  a 
cylindrical  volume  in  a  fluid  but  does  not  state  the  inhomogenous  wave 
equation  and  essentially  treats  the  acoustic  source  volume  as  a  radiating 
cylindrical  surface.  Morse  and  Ingard®  treat  the  subject  of  a  generic  heat 
source  but  do  not  address  the  case  of  ohmic  heating.  Westervelt  and 
Larson^  discuss  optoacoustic  laser  generated  sound  in  a  fluid,  which  is  an 
electromagnetic  heating  process  at  optical  frequencies.  The  derivation  to 
follow  accounts  for  the  ohmic  thermoacoustic  heat  source  through  the 
continuity  equation  as  an  added  mass  term. 

Assuming  a  constant  pressure  process,  a  change  in  density 
results  from  added  energy  when  a  fluid  is  heated.  This  change  in  density, 

Ap,  can  be  related  to  a  small  temperature  rise  AT  of  the  medium  by  the 

coefficient  of  thermal  expansion  p, 

Ap  =  p  p  AT.  (2.3) 

Here  p  is  the  total  fluid  mass  density. 
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The  derivation  of  the  continuity  equation  follows  from  the 
following  statement. 


Time  rate  of 
change  of  mass 
in  the  control 
volume 


Rate  of  mass  Rate  of  mass  Time  rate  of  change 

flow  into  the  -  flow  out  of  +  of  mass  production 
control  volume  the  control  within  the  control  volume 


The  incremental  mass  flow  rate  into  a  control  volume  V  across  its 
surface  S  is  -pu  •  dS,  the  minus  sign  resulting  fr.jm  the  unit  normal  vector  to 
the  surface  S  pointing  away  from  the  control  volume  (see  Fig.  2.1).  Net  mass 
flow  into  and  out  of  the  control  volume  is  expressed  as  the  integral  of  the 
incremental  mass  flow  rate  over  the  entire  closed  surface  of  the  control 
volurr'}. 


Net  mass  flow  rate  into  control  volume  =  - 


pu  •  dS. 


j 

S 


The  time  rate  of  change  of  mass  in  the  control  volume  due  to  heating  is 
defined  as  (see  Fig.  2.1) 


8  r  8 

Ap  dv  =  p  — 
8t 


8t 


J 

V 


p  AT  dv . 


(2.4) 


j 

V 


Thus  the  continuity  equation  takes  the  following  form: 


8 

^  J 


p  dv  = 


pu-  dS  +  p  — 
8t 


p  AT  dv . 


(2.5) 


J 

V 


Equation  (2.5)  is  valid  so  long  as  p  is  independent  of  time  and  temperature. 
These  assumptions  imply  small  temperature  changes.  ( In  general  p  is  a 


thermodynamic  variable  which  is  dependent  on  temperature  and  pressure.) 

The  temperature  chanop  can  now  be  related  to  the  Ohmic 
heating  mechanism  by  the  first  law  of  thermodynamics.  If  the  assumption  is 
made  that  the  heat  energy  added  to  the  fluid  is  an  isentiopic  process  and 
occurs  at  constant  pressure  then  the  temperature  rise  in  the  fluid  can  be 
related  to  the  electromagnetic  field  energy  input  by  the  coefficient  of  specific 
heat  Cp  (joules/(kg  K)).  On  a  per  unit  volume  basis  the  relation  is 


Energy  inputg^,  = 


P  Cp(T)  dT . 


j 


(2.6) 


For  small  changes  in  temperature,  AT./Tq  «  1,  Cp  is  approximately  constant 
so  Eq.  (2.6)  can  be  written  as 


Energy  Inputemf  =  P  Cp(T^)  (T,  -  T J  =  p  Cp(TJ  AT  .  (2.7) 


From  Jackson®  the  rate  of  conversion  of  electromagnetic  energy 
(power)  per  unit  volume  into  thermal  power  is  given  as 

Power  Inputg^,  =  J-E.  (2.8) 

Substituting  for  E  using  Ohm's  law,  Eq.  (2.2),  in  Eq.  (2.8),  then  taking  the 
time  derivative  of  Eq.  (2.7)  and  equating  the  result  with  Eq.  (2.8)  the 
relationship  between  the  time  rate  of  chanye  of  density  and  temperature  and 
current  density  is 

JJ 

(pAT),  = -  ■  (2.9) 
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Substituting  Eq.  (2.9)  into  Eq.  (2.5)  results  in  the  integral  form  of 
the  o’  Tiic  heating  continuity  equation, 


a 

at 


pdv  =  - 


pu-  dS 


V 


A  ^ 

aC  J 


JJ 


dv 


(2.10) 


The  surface  integral  of  Eq.  (2.10)  can  be  converted  to  a  volume 
integral  by  use  of  Gauss's  divergence  theorem.  This  results  in  the  differential 
form  of  the  thermoacoustic  continuity  equation, 


p,  +  V-  (p  u) 


(2.11) 


Note  that  J-J  in  all  the  preceding  equations  is  a  notational  form  for  the  real 
component  of  the  electrical  power  delivered  to  the  medium  since  it  is  the  real 
power  which  heats  the  medium  and  produces  thermoacoustic  radiation.  J*J 
is  meant  to  imply  the  following  operation  on  J,  assuming  a  time  harmonic 
dependence  of  e'i“‘. 

J(x,t)  =  Re[  J(x)  e-i"*]  =  (1/2)  [J(x)  e-j“‘  +  J(x)*  ei®']  (2.12) 

Then, 

J(x,t)-J(x,t)  =  (1/4)[J(x)  e-i®‘  +  J(x)  *  ei®*]-[  J(x)  e'i®'  +  J(x)  *  eh  (2.13) 
or, 

J(x,t)-J(x,t)  =  (1/2)  Re[J(x)*-J(x)  +  J(x)-J(x)  e-i^®'].  '2.14) 


Equation  (2.14)  implies  that  the  power  dissipated  in  the  medium 
has  a  tir  e  independent  component  plus  an  oscillating  component.  This 
ccnesponds  to  the  average  power  plus  the  time  harmonic  fluctuation.  For  a 
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;eal  valued  source  volume  function  (J{x)  real)  Eq.  (2.14)  reduces  to 


J(x,t)-J(x,t)  =  (1/2)  J2(x)  Re[  1  +  e-j2“‘]  .  (2.15) 

For  the  Cartesian  system, 

J2(X)  Hj^2(x)  +  Jy2(x)+J^2(x).  (2.16) 


The  physical  implication  of  the  constant  term  in  Eq.  (2.15)  is  that 
the  power  flow  into  the  medium  is  always  positive.  Since  the  power  input  is 
related  to  p,  by  the  continuity  equation,  Eq.  (2.11),  the  waves  produced  by 
ohmic  heating  will  be  compressional.  Intuitively  this  is  a  reasonable  result 
since  the  derivation  assumed  an  isentropic  process,  which  implies  that  the 
heat  added  to  the  medium  is  not  dissipated  but  raises  the  temperature  of  the 
medium,  and  therefore  causes  a  continuous  thermal  expansion  of  the 
medium. 

Integrating  Eq.  (2.9)  over  a  pulse  time  period,  x,  the  constraint  on 
energy  transport  into  the  fluid  such  that  the  small  temperature  rise 
assumption  is  valid  is 


1 


X 

JJ  dt 

0 


«  1 


(2.17) 


The  temperature  change  constraitit  of  Eq.  (2.17)  can  be  evaluated  using 
Eq.  (2.15)  if  the  change  in  density  is  assumed  small.  For  an  integer  number 
of  current  cycles  the  constraint  relation  for  thv^  current  density  pulse  length 
time,  X,  becomes 
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(2.18) 


V  '>:  «  2  oCp  \  . 


The  units  of  are  joules/metei^  whicn  is  the  volumetric  internal  energy 
density  change  of  the  electrolyte  by  the  applied  current  density  applied  by 
a  sinusoidal  signal  for  a  period  x.  The  right  side  of  Eq.  (2.18)  is  computed  to 
be,  for  a  S%  NaCI-water  solution  at  20°C,  2.5  X  10®  joules/meter®  or  2.5 
joules/cm^.  During  experimental  investigation,  described  in  Chapter  5,  the 
pulse  time  was  5  ms;  therefore  the  constraint  on  the  current  density  was 
«  1.62  X  10®  A/m2;  the  experimental  value  was  15000  A/m^. 

C.  DERIVATION  OF  THE  MOMENTUM  EQUATION 

The  momentum  conservation  la'*,  is  basically  a  statement  of 

Newton’s  second  law: 


time  rate  of  change  of  net  momentum  influx  sum  of  the  forces 

momentum  inside  the  =-  across  the  control  +  acting  on  the 
xntrol  volume  volume  surfaces  control  volume 


The  integral  form  of  the  three-dimensional  magnetohydrodynamic 
momentum  equation  follows  directly  from  the  verbal  statement  (see  Fig.  2.2). 

a  "  ‘  ‘  ' 


at 


p  u  dv  =  - 


j 

V 


(pu)  u-  dS 


P  ds  + 


JXB  dv  . 


(2.19) 


j 

V 


Gauss's  divergence  theorem  is  once  again  applied  to  the  surface  integral 
and  the  differential  form  of  Eq.  (2.19)  follows,^® 

p,u  +  pu,  +  V-  (pu)u  +  (pu-V)u  +  VP  =  JXB  .  (2.20) 

The  underlined  terms  of  Eq.  (2.20)  are  replaced  using  the  continuity 
equation,  Eq.  (2.1 1).  The  result  is  the  magnetohydrodynamic-thermoacoustic 
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FIGURE  2.2 

INTEGRAL  MOMENTUM  EQUATION  CONTROL  VOLUME 
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momentum  equation, 


pu,  +  (pu-V)u  +  VP  - (J-J)u  +  JXB.  (2.21) 

aCp 

Morse  and  Ingard^^  discuss  the  existence  of  electromagnetic 
forces  on  a  fluid  but  do  not  derive  an  explicit  expression  for  the  Lorentz  force 
in  connection  with  the  momentum  equation. 

D.  SMALL  SIGNAL  ACOUSTIC  APPROXIMATION 


The  small  signal  acoustic  approximation  (Mach  number  «  1)  is 
made  in  the  wave  equation  and  boundary  condition  derivations.  The  small 
signal  acoustic  perturbation  variables  of  first  order  magnitude  are  the 
perturbed  pressure  and  density, 

p  =  P  -  Po  .  (2.22a) 

and 

5p  =  P  *  Po  .  (2-22b) 


where  P^  and  p^  are  the  quiescent  pressure  and  density  respectively.  The 
small  signal  approximation  constrains  the  perturbation  variables  to  the 
following  limits,^ 


and, 


iPl  « 


|5p| «  Po 

The  equation  of  state  is  given  as 


e  =  constant,  P  =  P^,  p  =  p^ 


(2.23a) 

(2.23b) 


(2.23c) 


where  c^  is  the  isentropic  sound  speed,  defined  as  a  constant  for  the  fluid  at 
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constant  entropy,  e,  and  for  small  perturbations  in  total  pressure  P  and  total 
fluid  density  p  about  the  ambient  values. 

The  linearized  small  signal  equation  of  state  can  be  stated  in 
terms  of  the  perturbation  variables  as 

p  =  Co2  5p.  (2.24) 

Linearized  small  signal  versions  of  Eqs.  (2.12)  and  (2.20)  are,  respectively, 

P 

5p,  +  p^V-  u  =  - (J-J)  (2.25) 

and 

PqU,  +  Vp  =  JXB  .  (2.26) 

Performing  the  asual  mathematical  operations  on  Eqs.  (2.24),  (2.25),  and 
(2.26)  the  MHD-thermoacoustic  wave  equation  results, 

1  -P 

V2p  -  p„  =  -  (J-J),  -I-  V-  (JXB)  ,  (2.27) 

where  the  inhomogeneous  source  terms  are  the  thermoacoustic  source  and 
the  MHD  source,  respectively.  Pierce  identifies  the  first  source  term  as  a 
monopole  source  distribution  function  and  the  second  source  term  as  a 
dipole  source  distribution  function  where  JXB  is  the  dipole  moment 
vector.  ■'2 

The  simplest  control  volume  to  view  the  qualitative  nature  of  the 
two  sources  is  that  of  a  sphere,  such  as  shown  in  Fig.  2.3.  Thermal 
expansion  of  the  control  volume  due  to  a  uniform  current  density  flowing 
across  the  volume  will  generate  a  uniform  spherical  wave.  The  MHD  force 
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EXPANDING  AND  TRANSLATING 
CONTROL  VOLUMN 


(■PHASE) 


WAVEFRONT  PRODUCED  DUE 
TO  THE  MHD  FORCE  CAUSING 
A  TRANVERSE  MOTION  OF  THE 
CONTROL  VOLUME 
(+  PHASE) 


WAVEFRONT  OF  CONSTANT 
PHASE  FROM  THERMAL 
EXPANSION  OF  CONTROL  VOLUME 


FIGURE  2.3 

MHD  FORCE  AND  THERMOACOUSTIC  VOLUME  EXPANSION 
ACTING  IN  A  SPHERICAL  VOLUME 


ARL:UT 
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causes  a  translational  oscillation  of  the  control  volume  in  the  direction  of  the 
dipole  moment,  creating  a  dipole  wavefield. 

E.  ONE-DIMENSIONAL  WAVE  EQUATION 

In  the  one-dimer'sional  development  for  the  Cartesian  system  the 
vector  field  J  x  B  is  in  the  x  axis  direction.  The  bold  face  notation  is  thus 
suppressed  below  and  the  cross-product  implies  J  B  sin0,  where  0  is  the 
angle  between  the  vector  fields  in  the  y-z  plane. 

1  -p 

P«x--r  P„  =— (J-J),  +  (2.28) 

The  one-dimensional  wave  equation  is  the  basis  for  the  plane 
wave  field  assumed  in  the  MHD  transducer  waveguide.  Chapter  3  is 
devoted  to  the  derivation  of  the  boundary  conditions  and  solution  of  the 
waveguide  field. 

F.  DISPERSION  OF  THE  ACOUSTIC  FIELD  DUE  TO  THE  INDUCED 
CURRENT  DENSITY  IN  A  HIGH  CONDUCTIVITY  FLUID 

In  the  above  derivation  of  the  wave  equation,  induced  currents 
were  ignored.  When  one  considers  a  high  conductivity  fluid,  such  as 
mercury,  which  has  a  conductivity  on  the  order  of  10®  t/ohm  m  such  an 
assumption  must  be  carefully  examined.  As  in  the  NaCI-water  case,  the 
displacement  term  can  be  neglected,  and  Ohm’s  law  can  be  written  as 

J  =  a{E  +  uXB)  .  (2.29) 

However,  the  induced  current  density  term  cannot  be  eliminated.  If  again  we 
assume  that  the  current  density,  J,  is  15000  A/m^  and  the  particle  velocity  is 
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150  m/s,  E  is  then  found  to  be  150.5  V/m  for  the  me.^ury  medium.  The 
induced  electric  field,  uXB,  is  150  V/m  for  a  magnetic  field  of  1  T.  Therefore 
the  induced  term  should  not  in  general  be  neglected  in  the  derivation  of  the 
acoustic  wave  equation. 

When  Eq.  (2.29)  is  substituted  into  the  inhomogeneous  terms  of 
the  wave  equation,  Eq.  (2.27),  a  frequency  dependent  dispersion 
relationship  arises  vvhich  depends  on  the  magnetic  field  and  fluid 
conductivity.  The  dispersion  relation  is  derived  in  Appendix  A.1  for  plane 
waves  propagating  in  a  conducting  medium  in  a  direction  normal  to  an 
applied  static  magnetic  field.  The  dispersion  relation  is  stated  below. 


k  = 


jcoo 


\\%  J 


Bo^'l 


Po‘=o 


1/2 


(2.30) 


Physically  the  effect  of  the  magnetic  field  is  to  dampen  or  attenuate  the  wave 
by  converting  the  kinetic  energy  of  the  wave  to  electrical  energy  which  is 
then  dissipated  by  the  resistivity  of  the  medium.  From  Eq.  (2.30)  it  can  be 
seen  that  the  dispersion  effects  of  the  magnetic  field  can  be  neglected  when 

aB2 

—  «  1  (2.31) 

“Po 


Consider  the  6%  NaCI-water  solution  used  in  the  experiment, 
assuming  a  conductivity  of  5.28  mho/m  and  a  magnetic  field  of  1  T.  The 
dispersion  is  negligible  for  frequencies  much  greater  than  0.0007  Hz. 
Neglecting  dispersion  in  NaCI-water  solutions  is  thus  justifiable.  However, 
for  mercury  which  has  a  conductivity  of  ~10®  mho/m,  the  frequency  must  be 
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much  greater  than  74  Hz  before  dispersion  may  be  neglected. 

The  significance  of  the  dispersion  relation  arises  in  the  derivation 
of  the  transducer  impedance  (Chapter  4),  in  that  it  indicates  whether  the 
induced  electric  field,  and  thus  the  acoustic  coupling  contribution  to  the 
impedance,  will  be  important. 
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CHAPTER  3 


THEORY:  SOLUTION  OF  THE  WAVE  EQUATION  FOR  THE  ACOUSTIC 
FIELD  WITHIN  THE  WAVEGUIDE 


The  solution  of  Eq.  (2.28)  for  the  pressure  wavefield  excited  in  the 
waveguide  by  the  magnetohydrodynamic  and  thermoacoustic  sources  is 
found  by  solving  the  two  point  boundary  value  problem  for  this  equation  and 
the  related  boundary  conditions  at  the  apertures. 

The  inhomogeneous  terms  are  assumed  to  have  a  time 
harmonic  dependence  of  the  form 


(JXB),  =  (J(x)  X  B)^  e-i“'  (3.1) 

and 


-p(J-J)t  jcopJ'M  , 

-  =  -  e'i2“‘  .  (3.2) 

oCp  aCp 

Notice  that  the  thermoacoustic  source  responds  at  twice  the  dri\ 
frequency,  co,  and  is  therefore  twice  the  MHD  source  frequency. 

Due  to  the  linearity  of  the  wave  equation,  Eq.  (2.33),  the  solution 
of  the  wavefield  can  be  found  by  considering  the  two  source  mechanisms 
independeruly.  The  acoustic  fields  due  to  the  sources  can  then  be  added  to 
find  the  total  field. 


A.  BOUNDARY  CONDITIONS  OF  THE  TRANSDUCER  WAVEGUIDE 

The  boundary  conditions  of  the  waveguide  are  a  core  problem  of 

this  thesis.  By  specifying  the  boundary  conditions,  the  form  of  the  field 
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solution  is  specified  and  thus  the  transduction  nature  of  the  waveguide 
system  is  determined.  Hence  it  is  essential  that  the  boundary  conditions  of 
the  waveguide  be  accurately  stated. 

The  apertures  of  the  waveguide  are  assumed  to  be  radiating  into 
an  acoustic  free  medium.  Since  the  power  transmitted  by  the  apertures 
radiates  into  the  entire  medium  it  is  evident  that  some  of  the  power  radiated 
by  one  aperture  will  illuminate  the  other.  Thus  there  is  coupling  of  acoustic 
power  between  the  waveguide  apertures.  The  specification  of  the  boundary 
conditions  must  therefore  take  the  radiation  coupling  into  account.  The 
aperture  acoustic  radiation  coupling  is  described  by  the  mutual  acoustic 
radiation  impedance. 

The  acoustic  pressure  in  the  waveguide  at  an  aperture  is 
assumed  to  be  related  to  the  aperture  particle  velocities  by  attributing  a  self- 
and  mutual  radiation  impedance  to  the  aperture.  In  order  to  describe  the 
aperture  as  an  impedance  boundary  the  aperture  is  modeled  acoustically 
as  a  fluid  "piston  layer",  to  which  an  acoustic  radiation  impedance  can  be 
attributed.  The  piston  layer  is  assumed  to  transmit  only  the  acoustic  particle 
velocity. 

Two  pressures  are  present  on  the  free  medium  radiating  surface 
of  the  piston,  one  pressure  is  due  to  the  impedance  Zg  of  the  free  medium  as 
seen  by  the  piston  (i.e.,  the  inertia  and  radiation  properties  of  the  "piston") 
and  the  second  is  the  pressure  represented  by  the  mutual  radiation.  See 
Fig.  3.1  for  a  qualitative  view  of  the  aperture  piston  layer  forces. 


FIGURE  3.1 

FLUID  "PISTON"  ACOUSTIC  APERTURE  MODEL 
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Given  the  model  described  above  the  boundasy  conditions  can 
be  symbolically  stated  as  follows. 

p(il)  =  u(a)Z3 -u(-4)Z„  (3.3a) 

and 

p{-i)  =  -u{-jl)Z3  +  u(i)Z^.  (3.3b) 

where  the  self-impedance,  Zg,  is  defined  as  the  ratio  of  the  free  medium 
piston  surface  pressure  p^y^  to  the  aperture  particle  velocity  u(±il)  in  the 
absence  of  the  radiation  at  the  second  aperture.  The  mutual  radiation 
impedance,  Z^,  between  apertures  is  defined  as  the  ratio  of  the  pressure 
caused  at  one  aperture  due  to  the  particle  velocity  at  the  other  aperture.  The 
minus  signs  appearing  with  u(-i)  arise  from  the  convention  of  defining  a 
positive  velocity  as  being  outwardly  directed  from  the  surface  of  the  sphere. 
Relative  to  the  origin  in  the  center  of  the  duct  a  positive  velocity  at  the  -  Jl  end 
is  negative  (see  Fig.  3.2). 

An  analytical  solution  for  the  free  medium  pressure  field  as  a 
function  of  aperture  velocity  is  possible  if  the  apertures  are  represented  as 
uniformly  vibrating  circular  pistons  on  a  rigid  spherical  baffle  as  shown  in 
Fig.  3.1 .  Note  that  if  the  apertures  are  not  circular  an  equivalent  area  circular 
piston  may  be  defined  and  the  circular  piston  analysis  used  provided  the 
wavelength  is  not  small  relative  to  the  equivalent  circle  radius,  or  when  tne 
product  of  the  free  medium  wave  lumber  and  the  aperture  radius  is  much 
less  than  one  (kR3«1). 

There  is  an  analytical  solution  for  free  medium  acoustic  field 
generated  by  circular  pistons  on  a  rigid  spherical  baffle  in  terms  of  the 


>■6 


SOURCE 


VOLUMN 


spherical  harmonic  function  expansion.  Sherman'’^'^'*  has  derived  analytical 
results  for  the  self-  and  mutual  radiation  impedances  of  an  arbitrary  number 
of  circular  sources  with  arbitrary  locations  on  a  spherical  baffle  using  the 
spherical  harmonic  function  expansion.  His  results  are  specialized  here  to 
the  case  of  two  sources  (apertures)  located  180°  apart. 

Two  items  should  be  mentioned  at  this  point  concerning  the  use 
of  the  spherical  harmonic  expansion  as  a  boundary  condition.  (1)  the 
expansion  is  in  terms  of  spherical  waves,  ei*^''/r,  and  the  conditions  at  the 
aperture  are  in  terms  of  plane  waves,  ei*^^.  The  assumption  concerning  the 
matching  of  the  fields  is  that  the  radius  of  curvature  of  the  surface  of  the 
aperture  (radius  of  the  spherical  baffle,  Rg)  is  great  compared  to  the  radius  of 
the  aperture,  i.e.,  that  the  aperture  piston  layer  surface  motion  is 
approximately  planar,  and  not  radial.  (2)  The  pressure  at  an  aperture  due  to 
the  mutual  radiation  is  calculated  as  an  average  pressure  of  the  spherical 
waves  over  the  aperture  surface  area.  The  implication  of  (1)  and  (2)  is  that 
the  model  is  valid  for  R,/R^«1  andkR_«1. 

The  self-impedance,  Z^,  and  the  mutual  impedance,  which 
are  in  terms  of  the  spherical  harmonic  function  expansion  of  the  free  medium 
pressure  field,  p^j,,^  ,  are  defined  as 


Zs  {±^) 


J  P±!tfm  ('■.V)  dS 
S(±il) 


u(±jl) 


dS 

S(±.!l) 


(3.4a) 
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Zm(-A) 


(3.4b) 


p±wm  (r.M')  dS 
S(±il) 


u(+il)  I  dS 
S(±.(l) 

with  similar  definitions  for  the  aperture  at  -il. 

The  free  medium  pressure  ^  function  of  u(i),  u(-il), 

spherical  baffle  radius  Rg,  aperture  radius  R^,  and  the  free  medium  acoustic 
impedance.  Due  to  the  symmetry  of  the  sources  the  pressure  is  spatially  a 
function  of  one  angle,  \y,  and  the  range,  r.  The  source  geometry  is  shown  in 
Fig.  3.2.  See  Appendix  B  for  the  functional  form  of 


B.  JUSTIFICATIONS  OF  THE  BOUNDARY  CONDITIONS 

Before  going  on  to  the  solution  of  the  wavefield  it  is  necessary  to 
show  that  the  mutual  radiation  pressure  is  a  significant  factor.  An 
experimental  investigation  was  carried  out  with  a  transducer  having  an 
effective  spherical  baffle  radius  of  0.05  m  and  an  aperture  radius  of 
0.0210  m.  The  actual  apertures  were  square.  The  fluid  medium  in  the 
waveguide  was  a  6%  NaCI-water  solution  with  Pp=1041  kg/m^  and 
Cq=1600  m/s.  The  free  medium  was  fresh  water. 

The  self-  and  mutual  impedances  were  computed  numerically  by 
the  computer  routine  MISP,  using  Sherman's  results,  and  the  real  and 
imaginary  components  graphed  as  a  function  of  frequency  in  Figs.  3.3(a),  (b) 
and  3.4(a),  (b),  respectively.  The  use  of  these  impedances  as  boundary 
conditions  for  the  above  geometric  dimensions  is  valid  for  frequencies  less 
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IMPEDANCE  -  SAYLS  X  10  IMPEDANCE  -  RAYLS  X  10 


FREQUENCY  -  Hz  X  10^ 


(a) 

REAL  MUTUAL  IMPEDANCE 


(b) 

IMAGINARY  MUTUAL  IMPEDANCE 

FIGURE  3.4 

MUTUAL  IMPEDANCE  OF  THE  MHD  WAVEGUIDE  APERTURES 
t  =  0.10  m,  gg  =  0.10  m.  Cq  =  1550  m/s,  R  =  0.05  m,  R„=  0.02  m  arlut 
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than  10,000  Hz. 


The  mutual  impedance  components  (Fig.  3.4)  display  a 
characteristic  oscillation  which  arises  from  constructive  and  destructive 
interference  as  the  pressure  wave  length  varies  relative  to  the 
circumferential  distance  between  apertures.  For  example,  the  first  peak  in 
the  mutual  impedance  should  be  expected  at  the  frequency  at  which  a 
quarter  wavelength  is  equal  to  half  the  sphere  circumference;  for 
Cq=  1500  m/s  this  occurs  at  approximately  2390  Hz.  In  Fig.  3.4  the  first 
peak  of  the  imaginary  component  is  at  2700  Hz. 

To  show  that  the  inclusion  of  the  mutual  impedance  is  important  it 
is  necessary  to  compare  the  magnitude  of  the  mutual  impedance 
components  to  the  self-impedance  components.  Intuitively  one  would  expect 
the  mutual  impedance  to  be  most  significant  at  low  frequencies  where  the 
Rayleigh  distances  of  the  sources  are  very  short  and  therefore  radiation  is 
omnidirectional.  This  observation  explains  the  dropoff  at  high  frequencies  of 
the  envelope  of  the  mutual  impedance  plots. 

Figure  3.5  shows  the  ratio  of  mutual  to  self-impedance 
magnitudes.  The  peak  of  this  plot  occurs  at  approximately  7500  Hz  for  the 
geometry  described  above,  with  a  ratio  of  0.05.  One  might  conclude  that  the 
mutual  impedance  is  small,  however  the  radiation  of  sound  is  dependent  on 
the  real  component  of  the  acoustic  self-impedance.  Hence  it  is  more  relevant 
to  consider  a  ratio  of  the  real  impedance  parts.  Figure  3.6  shows  this  ratio, 
assuming  that  the  particle  velocities  of  the  duct  are  180°  out  of  phase  and 
are  of  equal  amplitude  (first  mode  of  the  waveguide  with  a  source 
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FREQUENCY  -  Hz  X  10 


FIGURE  3.5 

RATIO  OF  MUTUAL  TO  SELF-IMPEDANCE  MAGNITUDES 
0.10  m,  Sg  =  0.10  m,  Cg  =  1550  m/s,  R  =  0.05  m,  Rg  =  0.02  m 
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FREQUENCY  -  Hz  X  10 


FIGURE  3.6 

RATIO  OF  REAL  COMPONENTS  OF  THE  SELF- 
AND  MUTUAL  IMPEDANCE 

0.10  m,  Cg  =  0.10  m,  Cg  =  1550  m/s,  R  =  0.05  m,  Rg  =  0.02  m 


ARL:UT 

AS-86-441 


distribution  symmetric  about  the  x  origin).  Note  that  the  ratio  at  250  Hz  is 
nearly  one.  Thus  at  low  frequencies  the  real  component  of  the  mutual 
impedance  is  an  important  effect.  The  real  component  will  have  the  effect  of 
altering  the  amplitude  of  the  resonance  characteristics  of  the  wavefield  in  the 
waveguide,  and  to  a  much  lesser  extent  the  imaginary  component  will 
modify  the  periodicity  of  the  resonance  response  spectrum. 

C.  ACOUSTIC  WAVEFIELD  DUE  TO  THE  MAGNETOHYDRODYNAMIC 
SOURCE 

The  Fourier  transform  with  respect  to  time  of  the  inhomogeneous 
pressure  wave  equation  for  the  MHD  source  is  given  by 

p,,  +  k2p=(J(X)XB),.  (3.5) 

where  k  =  w/Cq  is  the  wave  number.  The  inhomogeneous  wave  equation  with 
inhomogeneous  boundary  conditions,  Eqs.  (3.3a)  and  (3.3b),  are  in  terms  of 
two  acoustic  field  variables  p  and  u.  To  make  the  solution  process  more 
convenient  a  change  of  variables  from  pressure  and  particle  velocity  to 
velocity  potential,  <}),  is  made  as  follows. 

u  =  (})^  (3.6a) 

and 

p  = -po(l),  =  icopo({)  .  (3.6b) 

Substitution  of  the  above  transform  pair  into  Eqs.  (3.5),  (3.3a), 
and  (3.3b)  yields  the  two-point  boundary  value  problem  for  the  wavefield  in 
the  waveguide  in  terms  of  one  dependent  variable,  p.  The  inhomogeneous 
wave  equation,  Eq.  (3.5),  becomes 
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{3.7a) 


-j 

“Po 

subject  to  the  boundary  conditions, 

'  -icop^iKA)  +  -  $,(-A)Z„  =  0  (3.7b) 

and 

-I  0)  5.)  -  (|>,(-  il)  2,  +  (|)^(  il)  =  0  .  (3.7c) 

The  solution  of  Eqs.  (3.7)  may  be  achieved  by  use  of  a  Green's 
function.  This  approach  was  chosen  because  the  Green’s  function  is  easy  to 
determine  and  the  solution  lends  itself  to  convenient  numerical  computation 
for  arbitrary  source  functions.  The  inhomogeneous  Green's  wave  equation 
for  the  above  problem  is 

G,x(x.Xo)  +  k2  G(x,xJ  =  -  6(x  -  xj  ,  (3.8) 

wnere  G(x,Xq)  is  the  Green's  function  for  the  wavefield  as  a  function  of  x  due 
to  a  Dirac  delta  function  distributed  source  at  x^.  Shown  in  Fig.  3.7  is  the 
delta  function  source  relative  to  the  coordinate  frame  and  other  geometric 
lengths. 

A  traveling  plane  wave  solution  in  the  following  form  is  proposed. 

[  A  +  R  oei'^’^)  •  il  <  x  <  x„  <  il,  (3.9a) 

I  B  (e*'"’'  +  Rj^e-i''")  -il  <  x^  <  x  <  il  ,  (3.9b) 

where  R.  j  and  R^^  are  the  reflection  coefficients  at  the  apertures. 

The  proposed  solution  represents  plane  waves  propagating  to 
the  left,  e  and  to  the  right,  ei'^’^.  Equation  (3.9a)  is  the  wavefield  to  the  left 
of  the  source  discontinuity  at  x^,  w  lere  the  left  traveling  wave,  e■i'<^  is 
incident  on  the  -il  aperture  and  R.j^ei*^^  represents  the  reflected  right 
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DELTA  FUNCTION  SOURCE  REPRESENTATION 
FOR  THE  PLANE  WAVE  GREEN’S  FUNCTION  SOLUTION 
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traveling  wave.  Equation  (3.9b)  is  the  wavefield  to  the  right  of  the  source, 
including  a  right  traveling  wave  and  a  left  traveling  reflected  wave.  In  general 
R  is  a  complex  quantity  representing  the  phase  shift  and  attenuation  of  the 
wave  upon  reflection  from  the  aperture. 

Equations  (3.9a.b}  must  be  compatible  with  the  conditions  at  the 
source  where  the  values  of  the  two  solutions  must  be  equal  (continuous)  and 
the  first  derivatives  must  satisfy  the  slope  discontinuity.  The  continuity 
condition  is  stated  as  follows. 

G(x=Xp,Xq)  for  -il  <  X  <  Xq  <  il  =  G(x=Xq,Xq)  for  -  A  <  x^  <  x  <  it 

(3.10) 

The  discontinuity  condition  is  found  by  integrating  Eq.  (3.8)  with  respect  to  x 
and  taking  the  limit  as  the  integration  bounds  converge  to  the  source  at  x^. 
This  yields 

lim  {  G,(Xo  +  e,x„)  -G^(Xo-e,  xJ}  =-1.  (3.11) 

e-^o 


The  result  of  applying  the  source  conditions  stated  in  Eqs.  (3.10) 
and  (3.11)  to  the  general  solution,  Eqs.  (3.9a,b),  is  the  Green's  function  for 
the  waveguide  field. 


G(x.Xo)  = 


j(eK  +  R^e-K)  (e-)'"’'  +  R.^e)'^’') 

2k(1-RfR.j)  -0_<x<Xo<5,  (3.12a) 

j(e-K  +  R.^ei'^M  (e'*'"  +  R^e-i''’^) 

2k(1-R,R.;)  -.<Xq<x<C.  (3.12b) 


The  solution  of  the  wavefield  in  terms  of  the  velocity  potential  can  therefore 
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be  stated  in  integral  form  as 


1 


<^{x)  = 


(J  X  B)^  G(x,Xj,)  dx^ 


.  A 


(3.13) 


For  any  frequency  o),  Eq.  (3.13)  can  always  be  evaluated 
numerically  because  the  ^lection  coefficents  are  constant  multipliers  and 
can  be  factored  out  of  the  integration.  A  particularly  simple  source 
distribution  function  to  evaluate,  and  the  one  investigated  experimentally,  is 
that  of  a  rectangular  electromagnetic  field  source, 

JXB  =  JoBj,(H(x -a)-H(x-b)) ,  (3.14) 


where  H  is  the  Heaviside  function.  Note  that  a  and  b  are  the  x  coordinant 
limits  of  the  rectangular  distribution  defined  in  Fig.  3.7  and  also  in  Fig.  3.2  . 
With  Eq.(3.14)  as  the  source  function  Eq.  (3.13)  becomes 


4>{x) 


-HoB, 

Wpo 


il 

(5(Xo-  a)  -  6(Xo-  b))  G(x,xJ  dx^  . 

V 


(3.15) 


To  evaluate  Eq.  (3.15)  the  integration  interval  may  be  broken  into 
three  pieces;  to  the  left  of  the  source,  the  source  region,  and  to  the  right  of 
the  source.  The  result  is  the  generation  of  three  functions  which  are  valid  to 
the  left  of  the  source  volume,  within  the  source  volume,  and  to  the  right  of 
the  source  volume. 
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(5(Xo-  a)  -  5(Xo-  b))  G(x,xJ  dx^ 


(3.16a) 


4>(x)  = 


-jJoB 


o 


®Po 


-il  il<x<a<b  <jl, 


X 

r 

(5(Xo-a)-5(Xo-b))G(x,xJdXo  + 

V 

a 


b 


J  (5(Xq-  a)  -  5(Xo-  b))  G(x.xJ  dx^ 

X  -i).<a<x<b  <il, 

il 

(5{Xj,-  a)  -  5(x„-  b))  G(x,xJ  dx^ 


(3.16b) 


(3.16c) 


The  evaluation  of  the  integrals,  Eqs.  (3.16a,b,c),  are  given  below. 
In  the  notation  used  below,  Gl  denotes  the  Green's  function  given  in 
Eq.  (3.12a),  the  left  volume  solution,  and  Gr  denotes  that  given  in  Eq. 
(3.12b),  the  right  volume  solution. 


4>(x)  = 


-jJoBo 


“Po 


'  (GL(x,a)  -  GL(x,b)) 

-il  <x<a<b<  il 

(3.17a) 

(GL(x,b)  -  GR(x,a)) 

-il  <a<. x<b<  il 

(3.17b) 

.  (Gp,(x,a)  -  GR(x,b)) 

-il<a<b<x<il 

(3.17c) 

Equations  (3.17a,b,c)  represent  the  wavefield  to  the  left  of  the  source 
volume,  within  the  source  volume,  and  to  the  right  of  the  source  volume, 
respectively. 
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At  this  point  the  general  solution  to  the  MHD  wave  equation  is 
not  complete  since  the  reflection  coefficients,  Rg  and  R.g,  have  not  been 
determined.  However,  these  reflection  coefficients  may  be  found  by 
substitution  of  Eqs.  (3.17a,b)  into  the  boundary  conditions  given  by 
Eqs.  (3.3a,c).  The  result  of  this  substitution  is  the  following  quadratic 
equations  in  terms  of  the  reflection  coefficients. 

C|  +  C2R.  +  CsR,  =0  (3.18a) 

C4+  C5Rg+  CgRg^  =0  (3.18b) 

where  the  C's  are  complex  constants  (see  Appendix  C  for  the  explicit  forms). 

The  computation  of  the  roots  of  Eqs.  (3.18a,b)  was  performed  by 
the  subprogram  unit  REEL  using  the  IMSL  routine  ZQADC.  The  magnitude  of 
the  reflection  coefficient  must  be,  by  definition,  less  than  or  equal  to  one  so 
the  proper  root  choice  was  unambiguous.  A  plot  of  the  reflection  coefficient  is 
given  in  Fig  3.8  for  the  waveguide  dimensions  given  previously.  Notice  that 
the  reflection  coefficient  declines  with  increasing  frequency.  This  is  because 
the  real  component  of  the  self-impedance  is  increasing  with  frequency  and 
therefore  allowing  more  acoustic  power  to  be  transmitted  out  of  the  wave¬ 
guide  into  the  free  medium. 

The  solution  for  the  wavefield  in  terms  of  the  velocity  potential  is 
now  completely  described  by  Eqs.  (3.17a,b,c)  with  the  solutions  of  Eqs. 
(3.18a,b)  for  the  reflection  coefficients.  The  expressions  for  the  pressure  field 
and  particle  velocity  field  are  recovered  oy  the  transform  equations, 
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Eqs.  (3.6a,b),  and  expressed  below  in  terms  of  the  Green's  functions, 
Eqs.  (3.1 2a, b),  and  the  derivatives  of  the  Green's  functions. 

(ei^^o  +  R,^e-i‘«o)  -  R.^ei'^’^) 

2(1-RoRo)  -.d  <x<x.<  il,  (3.19a) 

Gx(x.Xo)  =1 

-(e-K  +  R.^eiNCe'*^  - 

2(1-R^Rj^)  -il<XQ<x<  d,  (3.19b) 


p(x)  = 


The  complex  pressure  wavefield,  p(x),  for  the  three  regions  of  the  transducer 
waveguide  as  a  function  of  the  axial  coordinate  x  are  as  follows. 

Jq  Bg  (GL(x,a)  -  GL(x,b))  -d<x<a<b<d 
Jq  Bq  (GL(x,b)  -  Gp,(x,a))  -d<a<x<b<d 

I  Jg Bq (Gpj(x,a) - GR(x,b))  -d<a<b<x<d 

The  complex  particle  velocity  wavefield,  u(x),  is  as  follows. 

(GxL(x,a)  -  GxL(x,b))  -d<x<a<b<d 

-  (GxL(x,b) - GxR(x,a))  -d  <a<x<b<  d 

I  (GxR(x,a) - GxR(x,b))  -d<a<b<x<d 


u(x)  = 


jJoBo 


«Pc 


(3.20a) 

(3.20b) 

(3.20c) 

(3.21  d) 
(3.21b) 
(3.21c) 


Several  plots  of  the  modal  particle  velocity  and  pressure  field 
amplitudes  as  a  function  of  the  waveguide  x  coordinate  and  frequency  are 
shown  in  Figs.  3.9,  3.10,  and  3.11.  Figure  3.9  is  for  the  case  relevant  to  the 
experiment  conducted.  Shown  in  Fig.  3.9(a)  is  the  modal  particle  velocity  for 
a  waveguide  sound  speed  of  1536  m/s.  It  is  important  to  note  the  step 
change  in  particle  velocity  at  the  aperture  where  there  is  a  step  change  in 
acoustic  impedance.  Notice  that  as  the  frequency  is  increased  the  particle 
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MAGNITUDE  OF  THE  MODAL  PRESSURE 


FIGURE  3.9 

MAGNITUDE  OF  THE  MODAL  PARTICLE  VELOCITY 
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5 


(a) 

MAGNITUDE  OF  THE  MODAL  PARTICLE  VELOCITY 


(b) 

FORWARD  VIEW 


(c) 


180‘  ROTATION  OF  (b) 
MAGNITUDE  OF  THE  MODAL  PRESSURE 


FIGURE  3.10 
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velocity  in  the  waveguide  and  the  step  change  in  velocity  at  the  aperture 
decrease.  The  decrease  in  waveguide  particle  velocity  with  increasing 
frequency  is  due  to  the  inertia  effect  of  the  fluid.  The  decreasing  step  change 
at  the  aperture  is  due  to  the  increasing  aperture  impedance  with  frequency 
(see  Figs.  3.3(a)  and  (b)). 

In  Figs.  3.10(a)  and  (b)  the  sound  speed  in  the  waveguide  was 
set  at  650  m/s  to  demonstrate  the  resonance  behavior  of  the  waveguide  at 
lower  frequencies.  At  approximately  9000  Hz  there  is  a  resonant  standing 
wave.  This  first  resonance  is  approximately  a  3/2  wavelength  mode  and  is 
the  first  mode  because  at  this  frequency  the  aperture  impedance  has 
become  large  enough  to  allow  a  "hard"  reflective  waveguide  termination  and 
thus  support  of  a  reduced  aperture  particle  velocity. 

Figure  3.11(a)  and  (b)  show  the  modal  behavior  generated  by  an 
MHD  source  which  is  0.05  m  and  centered  in  a  0.10  m  waveguide.  This  plot 
is  for  a  reduced  sound  speed  of  650  m/s,  but  what  is  exemplified  is  the 
discontinuity  of  the  wavefields  due  to  the  presence  of  the  MHD  volume  force. 
The  FORTRAN  program  MHD  has  the  capability  to  calculate  the  acoustic 
response  from  any  rectangular  MHD  source  distribution  anywhere  in  the 
region  of  the  defined  waveguide  length. 

D.  ACOUSTIC  WAVEFIELD  DUE  TO  THE  THERMOACOUSTIC  SOURCE 

Considered  first  was  the  pressure  wavefield  due  to  the  MHD 
source  mechanism.  The  pressure  field  contribution  of  the  thermoacoustic 
mechanism  is  now  derived.  The  Fourier  transform  with  respect  to  time  of  the 
inhomogeneous  wave  equation  for  the  thermoacoustic  source  is  given  by 


(3.22) 


Pxx 


+  k^p  = 


jnpj2(x) 

2crCp 


where  k  =  Q/Cq,  and  Q  is  defined  as  2  co.  The  solution  procedure  is  the  same 
as  that  performed  for  the  MHD  source.  The  thermoacoustic  inhomogeneous 
wave  equation,  Eq.  (3.22),  has  the  same  inhomogeneous  boundary 
conditions  as  the  MHD  wave  equation  solved  previously.  Again  a 
transformation  to  the  velocity  potential  is  made  using  Eqs.  (3.6a,b).  The  two- 
point  boundary  value  problem  for  the  thermoacoustic  source  generated 
wavefield  in  the  waveguide  in  terms  of  one  dependent  variable,  4),  is  thus 


,,  PJ'M 

‘J'xx  =  - -  • 

2ap,Cp 

subject  to  the  boundary  conditions, 

(3.23a) 

and 

-inp„4.(-i)  -  $,(i)z„,  =0. 

(3.23b) 

(3.23c) 

The  solution  of  Eqs.  (3.23a,b,c)  may  be  achieved  by  use  of  the 
Green's  function  derived  in  Section  3.3,  and  rewritten  below. 

j  {e^^o  +  Rj^e-K)  (e-i'«  + 

2k(1  -  Rj^R.g^)  -i).<x<XQ<  il,  (3.24a) 

G(x,x-)  =  i 

j  (e-K  +  R.j^eK)  (e)'«  + 

2k(1-Rj,R.^)  -il<Xo<x<  A.  (3.24b) 


The  solution  of  the  thermoacoustic  pressure  wavefield,  in  terms  of  the 
velocity  potential  and  the  above  Green's  function,  is  stated  in  integral  form  as 


4>(x) 


r 

j2(x)o  G(x.xJ  dx„ 


(3.25) 


The  source  distribution  function  investigated  in  the  experiment, 
is  a  rectangular  electromagnetic  field.  Therefore  the  current  density  has  the 
rectangular  distribution 

J(x)  =  Jj,  (H(x  -  a)  -  H(x  -  b)) ,  (3.26) 


where  is  the  current  density  amplitude  and  H  is  the  Heaviside  function. 
With  Eq.(3.26)  as  the  source  function  Eq.  (3.25)  becomes 

il 

No  f 

(|)(x)  =  - -  (H(Xo-  a)  -  H(Xo-  b))  G{x,x^)  dx^  .  (3.27) 

2apoC  J 

^  -A 


The  iterated  integration  of  Eq.  (3.27)  for  the  three  regions  of  the 
waveguide  are  given  on  the  next  page. 

f 

b 

(H(Xq-  a)  -  H(Xo-  b))  G(x,xJ  dx^ 

-il  il<x<a<b<il,  (3.28a) 

X 

(H(Xo-  a)  -  H(Xo  -  b))  G(x,x^)  dx^  + 


J  (H(x„-  a)  -  H(x^,-  b))  G(x.Xo)  dx^ 

X  -il  <a<x<b<il,  (3.28b) 
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il 

(H(x,-  a)  -  H{x„-  b))  G{x.xJ  dx„ 

-il<a<b<x<jl.  (3.28c) 

The  results  of  evaluating  the  integrals  for  the  velocity  potential  vvavefield  are 
as  follows  for  the  step  function  source. 

(.(eika  -  ei><b)  +  R^(e-i'<a  -  e-)kb))(e-i'«  + 

2k2(l-R^Rj)  “ 

-il<x<a<b<il,  (3.29a) 
(e-)'<3- R  4ei''3)(eik>«  + R^e-ji^) 


W  = 


2opoCp 


2k2(1-R^R5^) 

(ei''*'  ■  R^e-jkb)(e-)'«  +  R.^^ei"^) 


1 


2k2(1.Rj^Rj^)  k2 

•.(l<a<x<b<  il, 

((e-ika  .  e-jkb) .  R  (eika  .eikb))(eikx  +  r  g-ikx) 


(3.29b) 


2k2(l-R^R^) 


•il  <a<b<x  <1  (3.29c) 


Equations  (3.29a, b,c)  represent  the  velocity  potential  wavefield  to  the  left, 
within,  and  to  the  right  of  the  source  volume,  respectively. 

The  reflection  coefficients  are  found  by  substitution  of 
Eqs.  (3.29a,c)  into  the  boundary  conditions  given  by  Eqs.  (3.3a, b).  The 
result  of  this  substitution  is  the  following  quadratic  equations  in  terms  of  the 
reflection  coefficients, 

Ci  +  CgR  +  C3R  j^2  ^0  ,  (3.30a) 

C4  +  CgRj^  +  CeRj^2  ^0  ,  (3.30b) 
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where  the  C's  are  complex  constants  (see  Appendix  C  for  the  explicit  forms). 
The  computation  of  the  roots  of  Eqs.  (3.30a,b)  was  performed  by  the 
subprogram  unit  REFL  of  program  THERMO  using  the  IMSL  routine  ZQADC. 

A  plot  of  the  reflection  coefficient  magnitude  is  given  in  Fig.  3.8  for  the 
waveguide  dimensions  given  previously.  Note  that  Fig.  3.8  was  computed 
using  the  MHD  solution;  however  the  magnitude  of  the  reflection  coefficient 
is  the  same  for  both  cases;  the  difference  occurs  in  the  phase. 

The  solution  for  the  thermoacoustic  wavefield  in  terms  of  the 
velocity  potential  is  now  completely  described  by  Eqs.  (3.29a, b,c)  with  the 
solutions  of  Eqs.  (3.30a, b)  for  the  reflection  coefficients.  The  expressions  for 
the  pressure  field  and  particle  velocity  field  were  recovered  by  the  transform 
equations,  Eqs.  (3.6a,b),  applied  to  Eqs.  (3.29a,b,c)  and  these  are  stated 
below. 


p(x)  = 


2«JPoCc 


(.(eika .  eikb)  +  Rj^(e-ika .  e-i''''))(e-i''’'  +  R.^^ei'"’') 


2k(1  - 


jl<x<a<b<il,  (3.31a) 


(e-ika  -  R.j^ej''3)(ei''’'  +  R^e'^) 


2k(1  -  R^R.^) 

(ei’^b  -  R4e-i'''')(e-i'«  +  R.j^ei^’') 


1 


2k(1  -  Rj^R  j^) 


jl<a<x<b<  il, 

((e-jka .  Q-jkb) .  R  ^(eika  .eikb))(eikx  R^e-i'"’') 


2k(1-RjR  j) 


(3.31b) 


-il<a<b<x<il.  (3.31c) 
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u(x)  is  the  complex  particle  velocity  field  in  the  transducer  waveguide. 


((eika  -  -  R^(e-i''3  -  e-)''^)){e-i''='  -  R.^ei'"’') 

2k{1  -  li)  ~ 

-il<x<a<b<  il, 

(e-i''^  -  R.4ei''®)(ei'<’‘  - 

2k(1  -  R^R  i) 

(e)''^  -  R^e-j''‘'){e-j''’‘  - 

2k(l  -  R,R  ,)  ”  ^  ^ 

-il<acx<b<3l, 

((e-jka .  e-jkb) .  R_^(eika  -eikb))(ei'«  -  R^e’i'''') 

2k(1  -  R^R  i)  ^ 

-Jl<a<b<x<  L 


(3.32a) 


(3.32b) 


(3.32c) 


Unfortunately  time  did  not  permit  the  author  to  write  the  computer 
code  necessary  to  generate  the  three-dimensional  plots  of  the  pressure  and 
particle  velocity  modes,  as  was  done  for  the  MHD  mechanism. 


CHAPTER  4 


THEORY;  ELECTROACOUSTIC  RELATIONS 

This  chapter  is  concerned  with  the  determination  of  the 
transduction  characteristics:  input  electrical  impedance,  real  and  reactive 
input  electrical  power,  power  efficiency,  current  and  voltage  transmitting 
sensitivity  for  the  MHD  acoustic  source  mechanism,  and  thermoacoustic 
source  mechanism. 

A.  POLARIZATION  IMPEDANCE  OF  THE  MHD  WAVEGUIDE  ELECTRODES 
In  general  when  considering  the  impedance  of  an  electrolyte 
fluid  contained  in  a  cell  with  metal  electrodes  two  independent  processes 
must  be  analyzed:  the  frequency  dependent  nature  of  the  electrolyte  and  the 
electrical  characteristics  of  the  electrolyte  interaction  with  the  electrode 
surface. 

The  conductivity  of  an  electrolyte  not  in  the  vicinity  of  an  electrode 
surface  is  "classically"  analyzed  using  the  theories  of  Debye  and  HuckeP^ 
and  Debye  and  Falkenhagen.''®  The  conductivity  for  a  1,1  valent  electrolyte, 
such  as  NaCI  in  water,  is  dependent  on  the  applied  electric  field  strength, 
concentration  of  NaCI,  and  frequency.  However,  for  NaCI  in  water  at  1  M 
concentration  with  electric  fields  less  than  1X10^  V/cm  and  frequencies  less 
than  18  MHz,  the  conductivity  is  constant  (  a  =  5.28  1/(ohm  m)). 

At  an  electrode-elecirolyta  interface  a  potential  difference  is 
present.  The  potential  arises  from  the  electrons  on  the  metal  surface  and  the 
ions  in  the  electrolyte  attracting  each  other  according  to  Coulomb's  law.  This 
charge  distribution  at  an  interface  is  known  as  the  electrical  double  layer  in 


the  electrochemical  literature.''^  The  important  concept  is  that  the  mobile  ions 
have  a  tremendous  capacity  for  charge  storage  and  thus  a  very  significant 
capacitive  impedance  when  a  time  harmonic  signal  is  applied  to  the 
electrode. 

For  the  case  of  planar  electrodes,  of  interest  in  this  investigation, 
the  Gouy-Chapman'®  model  of  the  double  layer  capacitance  is  given  as,''® 

e„  F2  z  N  (b-a)  L 

Cd=  - ^  (4.1) 

4kRT 

where 

Cq  =  double  layer  capacitance,  F, 
e^  =  dielectric  permitivity,  8.85x10*'^  F/m 
F  H  Faraday  constant,  96524  C/mole, 

N  s  concentration  of  electrolye,  moles/m^, 

R  s  ideal  gas  constant,  8.31  joules/°K  mole, 

T  s  absolute  temperature,  293  °K  (20°  C) 
z  s  charge  valence  of  the  ionized  molecule,  1 ,  and 
(b  -  a)  Ay  =  electrode  area. 

For  a  1  M  NaCl  electrolyte  solution  the  Chapman  model  predicts  the 
capacitance,  on  a  per  unit  area  basis,  to  be  52  mF/m^.  Measurements  made 
by  Vetter^®  indicate  a  range  of  values  of  100-400  mF/m^.  This  will  be 
discussed  more  in  Chapter  5  in  re^^ard  to  the  value  measured  for  the  MHD 
transducer. 

In  addition  to  ihe  double  layer  capacitance  of  the  electrode 
interface  there  is  an  additional  impedance  effect  observed  due  to  an  ion 
concentration  gradient  near  the  electrode.  This  is  due  to  electrochemical 
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reaction  processes.  The  transfer  of  charge  at  the  interface  causes 
electrochemical  reactions  which  result  in  the  depletion  of  the  orig!  lal 
electrolyte  ions  and  the  production  of  other  ions  and  molecules.  Since  the 
ion  concentration  gradient  represents  a  charge  gradient  one  can  expect  an 
electric  potential  to  result.  The  Nernst  equation^^  describes  this  potential, 

RT  rN(o.t)^ 


Vd=  —  In 


(4.2) 


nF 


V 


N, 


) 


in  terms  of  N(o,t),  the  ion  concentration  at  the  electrode  surface  (^=0)  as  a 
function  of  time,  and  N^,  the  equilibrium  ion  concentration  "far"  from  the 
electrode,  n  is  the  number  of  moles  of  electrons  required  to  form  one  mole  of 
reaction  product. 

In  order  to  apply  Eq.  (4.2)  the  electrochemical  reactions  at  the 
electrodes  must  be  determined.  The  kinetically  acceptable  electrode 
reactions  for  NaCI  and  water  are 

2Cr Cl2(gas)  +  2e‘  ,  (4.3) 


which  is  the  anode  or  oxidation  reaction.  This  reaction  represents  the 
conversion  of  the  chloride  ion  to  chlorine  gas  with  the  donation  of  electrons 
to  the  electrode.  The  other  electrode  reaction  is  the  cathode  or  reduction 
reaction, 

2H2O  +  2e-  ->  H2(gas)  +  20H-.  (4.4) 


Note  that  dissociation  of  water  in  the  cathode  reaction  is  preferred  over  the 
reduction  of  sodium  ion  to  solid  sodium  (Na"*"  +  e’  -»  Na(solid)).  The  overall 
chemical  reaction  can  be  stated  as 
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(4.5) 


2Na+  +  2C|-  +  2Hp  ->  2Na+  +  20H-  +  H2(gas)  +  Cl2(gas). 

It  can  be  seen  from  the  reaction  equation,  Eq.  (4.5),  that  the 
chloride  ion  will  experience  a  concentration  gradient.  Note  that  water  does 
not  contribute  to  the  ion  concentration  gradient  of  the  electrolyte  since  water 
is  not  ionized  and  is  the  solvent  and  is  thus  in  great  supply  in  the  vicinity  of 
the  electrodes.  Since  the  hydroxide  ion  is  produced  one  might  ask  whether  it 
is  a  more  preferable  oxidation  reaction  molecule  than  the  chloride  atom. 
Kinetically  it  is  not. 

One  last  item  needs  to  be  addressed  with  regard  to  the  formation 
of  hydrogen  and  chlorine  gas.  The  reaction, 

H2(gas)  +  Cl2(gas)  ->  2HCL(gas),  Ag^  =  -45  kcal/mole  (4.6) 

is  a  spontaneous,  very  exothermic  reaction  as  evidenced  by  the  large 
reduction  in  the  Gibbs  free  energy,  Ag^,.  However,  as  a  practical  matter  the 
gases  do  not  always  combine  but  can  nucleate  on  the  electrodes.  The  result 
of  the  gas  formation  and  the  effect  on  operation  of  the  transducer  is 
discussed  later  in  this  chapter  and  in  Chapter  5. 

The  concentration  distribution,  N(^,t),  for  the  evaluation  of  the 
Nernst  equation,  Eq.  (4.2),  is  based  on  the  concentration  of  the  chloride  ion 
as  a  function  of  the  rate  at  which  it  is  oxidized.  vSince  the  boundary  of  the 
electrolyte  is  a  plane  electrode  surface,  which  is  the  boundary  of  the 
waveguide  of  the  same  cross-section,  the  concentration  wHI  only  vary  in  the 
y  direction  of  the  waveguide  coordinate  system.  Therefore  the  distribution 
should  only  be  a  function  of  one  spatial  dimension,  The  diffusion  equation 


56 


is  stated  as 


(Pick's  second  law)  N(^,t),  =  D  N(^,t)^t ,  (4.7a) 


subject  to  the  following  boundary  conditions. 

J  V 


(Pick's  first  law) 

N(o,t)  =  e-)“- 

nFD 

(4.7b) 

and 

lim  N(^,t)=N,, 

?“>oo 

(4.7c) 

and  the  initial  condition. 

N(^,o)  =  N,. 

(4.7d) 

In  the  above  equations,  D  is  the  diffusion  constant,  J  is  the  current  density 
amplitude  defined  previously,  and  v  is  the  number  of  moles  of  ions  needed 
to  produce  cne  mole  of  product  (v  =  2  for  the  chloride  ion  reaction).  The 
solution  of  Eq.  (4.7a)  is  given  as^^ 

J  V  , _  , _ 

N(^,t)  =N  +  _  e-Vw/2D^  g.j(a)t  -  Vay2D  Jt/4)  _  (4 

°  n  F  VDw 

Equation  (4.8)  represents  an  exponentially  damped,  dispersive, 
outgoing  diffusion  wave.  However,  one  might  not  accept  the  boundary 
condition,  Eq.  (4.7c),  and  thus  the  solution,  Eq.  (4.8),  as  being  a  valid 
solution  to  the  boundary  value  problem,  which  actually  involves  not  an 
infinite  positive  domain,  but  a  finite  domain  problem  with  another  electrode 
at  ^=ily.  From  the  coefficient  e‘ ^  of  Eq.  (4.8)  it  is  seen  that  the 
concentration  distribution  is  spatially  exponentially  damped  and  that  for 
NCO/2D  ^  =  4  the  diffusion  wave  is  98%  decayed.  At  500  Hz  and  D  =  1.0X10'^ 
m\s,  the  diffusion  wave  is  damped  98%  in  0.003  mm,  which  is  a  very  small 
distance  compared  to  3.81  cm,  the  electrode  separation  used  in  the 
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experiment.  Hence  Eq.  (4.8)  is  a  good  approximation  for  this  investigation. 

The  Nernst  equation,  Eq.  (4.2),  can  now  be  evaluated  using 
Eq.  (4.8).  Note  that  in  order  to  linearize  ln( )  in  terms  of  the  current  density 
so  that  a  linear  relation  between  current  and  voltage  is  obtained,  it  must  be 
assumed  that  the  change  in  concentration  represented  by  the  second  term  of 
Eq.  (4.8)  is  small  compared  to  N^.  This  condition  is  met  if 

N^FnV^ 

J«  - - -  .  (4.9) 


At  500  Hz  the  right  side  of  Eq.  (4.9)  equals  1.71  X  10^  A/m^  .  The 
experiments  were  performed  at  a  current  density  of  1.5  X  lO'^  A/m^,  which  is 
a  factor  of  10  less  than  the  constraint  and  represents  approximately  a  5% 
linearization  error  of  ln( ).  Substituting  Eq.  (4.8)  into  Eq.  (4.2)  and  linearizing, 
the  Nernst  equation  becomes 


JRTv2 

V .  = _ e  . 

n^F^N^V^ 


(4.10) 


From  Eq.  (4.10)  the  diffusion  impedance  for  the  waveguide  electrode  can  be 
defined  as 


RTv2(1  -j) 

ily(b  -  a)  n2  f2  uj2W(x> 


(4.11) 


where  the  current  is  defined  as  J  Ay  (b  -  a). 

The  electrical  circuit  interpretation  of  Eq.  (4.11)  is  that  of  a 


resistor  and  capacitor  in  series,  with  the  resistive  and  capacitive  impedance 
having  the  equal  magnitudes  and  creating  a  current  lag  of  -45°.  Notice  that 
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the  resistive  component,  R^,  is  proportional  to  1/\w.  The  capacitive  nature  of 
Eq.  (4.11)  can  be  separated  from  the  resistive  by  assuming  the  impedance 
form  1/j  (oC^,,  which  yields  the  expression  for  C^,.  the  diffusion  capacitance, 
that  is  also  a  function  of  1/Vw, 


fly  (b  -  a)  n2  f2 
RTv^V^ 


(4.12) 


Finally  the  complete  electrolyte-electrode  circuit  topology^^  can 
be  assembled,  given  the  expressions,  Eqs.  (4.1)  and  (4.11),  for  the  values  of 
the  elements,  as  shown  in  Fig.  4.1.  The  diffusion  impedance  in  parallel  with 
the  double  layer  capacitance  is  termed  the  polarization  impedance,  Zp,  of 
the  electrode.  Note  that  in  Fig.  4.1  the  impedance  of  the  electrolyte  is 
represented  by  the  complex  impedance  element  Z.  This  is  done  so  that  the 
impedance  due  to  the  acoustic  coupling  can  be  included  in  a  complete 
circuit  model  of  the  input  electrical  impedance.  The  analysis  of  Z  is 
discussed  in  the  next  section  of  this  chapter. 

B.  ELECTRICAL  IMPEDANCE  OF  THE  MHD  TRANSDUCTION  PROCESS 
The  electrical  impedance  of  the  transducer  due  to  the  interaction 
of  the  electrolyte  with  the  acoustic  field  is  defined  as  the  voltage  across 
divided  by  input  current  through  the  impedance  element  Z  of  Fig  4.1.  The 
input  current  I  is  found  by  integrating  the  current  density  over  a  surface  near 
the  electrodes  (four  diffusion  layers  away,  4/Vco/2D). 

I* 

I  =  J  •  dS  (4.13) 

S 
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FIGURE  4.1 

EQUIVALENT  CIRCUIT  OF  THE  MHD  INPUT  IMPEDANCE 
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The  voltage  is  defined  as  the  path  integral  of  the  electric  field  between  the 
electrodes,  but  not  near  the  surface, 


m 


V  = 


E-  d! 


path  between 
electrodes 


(4.14) 


Using  the  above  integral  definitions  the  electrolyte  impedance  is  defined  as 
follows. 


E-  dl 


path  between 
electrodes 


z  = 


J  •  dS 


(4.15) 


.'vV 


Evaluation  of  the  numerator  of  Eq.  (4.15)  follows  from  assuming 
that  the  surface  integration  near  the  electrodes  is  at  a  constant  electric  field 
value. 

Substituting  Eq.  (2.28)  into  Eq.  (4.15),  to  take  account  of  the 
induced  electric  field,  and  substituting  Eq.  (3.21c)  for  u,  the  impedance  is 


dy 


Z  = 


j 

0 


XF 


(4.16) 


a 

^0 


(1  +  j  aB^2  (GxL(x,b)  -  GxR(x,a)))  dx 


0  a 


61 


-'-'C''-' 


HI  )*1  \  H  , 


Performing  the  integration, 


Z  = 


'  a  il^  (b-a)  j  g2  6^2 

- 5 - '  r.n  T  (^L(b.b)  -  GJa.b)  -  Gplb.a)  +  GR(a.a)) 

^  Xy  0)  Pq  34y 


(4.17) 


Shown  in  Figs.  4.2(a)  and  (b)  are  the  real  and  imaginary 
components  of  the  electrolyte  impedance  for  the  waveguide  filled  with  a  6% 
NaCI-water  solution.  The  transducer  boundary  conditions  of  the  waveguide 
are  those  of  the  transducer  used  in  the  experiment  discussed  in  Chapter  5. 
For  the  6%  NaCI  solution  at  500  Hz  the  impedance  is  1.581  -  j  0.2X10’®. 
The  important  result  is  that  the  real  component  for  6%  NaCI  is  relatively 
constant  over  the  frequency  band  from  500  to  13000  Hz  and  the  imaginary 
component  is  negligible. 

From  the  preceding  discussions  it  can  be  concluded  that  the 
electrolyte  Impedance  can  be  assumed  to  be  real  and  constant  for 
NaCI-water,  except  at  "low"  frequencies  and/or  high  magnetic  fields  (B»1T) 
for  a  high  conductivity  fluid  such  as  mercury.  It  is  justifiable  therefore  to  state 
Ohm's  law  as  Eq.  (2.2)  for  the  electrolyte  impedance. 
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IMPEDANCE  -  ohm 


1.580768310 
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(a) 

REAL  ADMITTANCE 


(b) 

IMAGINARY  ADMITTANCE 


FIGURE  4.2 

MHD  MECHANISM  IMPEDANCE  FOR 
6%  NaCI-H20  SOLUTION 
eg  =  0.10  m,  C  =  0.10  m,  Cq  -  1550  m/s,  R  = 


C.  MHD  TRANSDUCTION  PROCESS  POWER  EFFICIENCY 

The  transduction  power  efficiency  is  defined  as  the  real 

component  of  the  radiated  acoustic  power  divided  by  the  real  component  of 
the  input  electrical  power  plus  the  real  component  of  the  radiated  power. 
The  complex  input  power  expression  is 

We=  l^z,.  (4.18) 

where  Z,  is  the  transducer  terminal  impedance  defined  as 

Z,  =  2Zp  +  Z.  (4.19) 

To  calculate  the  total  radiated  acoustic  power  the  real  component 
of  the  acoustic  self-impedance  Re(Z5)  and  the  acoustic  particle  velocity 
evaluated  at  the  apertuie  are  integrated  over  the  surface  of  the  aperture. 

Wa=  j|u(il)|2  Re(Z3)  ds  +  f|u(-il)|2  Re(Z3)  ds  (4.20) 

aperture  apenure 

at  5.  at-t 

To  calculate  the  power  efficiency,  the  ratio  of  Eq.  (4.20)  to  the 
real  component  of  Eq.  (4.18)  is  taken, 

Wa 

T1  =  -  •  (4.21) 

Re(W3) 

Shown  in  Fig.  4.3(a)  and  (b)  are  plots  for  mercury  and  the  6%  NaCI  solution 
used  in  the  experiment,  respectively.  At  low  frequencies  the  efficiency  is 
higher  because  the  particle  velocity  is  higher  since  there  is  low  acoustic 
aperture  impedance  at  low  frequencies.  As  the  frequency  increases,  up  to 
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the  first  resonance,  the  particle  velocity  decreases  faster  than  the  real 
component  of  the  radiation  impedance  increases,  and  thus  there  is  a 
decreasing  trend  in  the  efficiency. 

D.  MHD  TRANSMITTING  CURRENT  SENSITIVITY 

The  transmitting  current  sensitivity  (TCS)  is  defined  as  the  ratio  of 

output  acoustic  pressure  at  1  m  at  an  angle  y  to  1  A  electrical  input  current. 
The  TCS  is  also  defined  for  the  special  case  of  the  transducer  waveguide 
axis  (x  axis  response,  \j/=  0°,  and  range  of  1  m).  The  expression  for  the 
farfield  pressure,  Eq.  (B.4),  is  analytically  and  computationally  simpler  than 
the  full  spherical  harmonic  function  expansion  used  to  express  the  mutual 
acoustic  impedances.  The  farfield  equation  has  the  Hankel  function 
components  replaced  by  the  large  argument,  kr»1,  asymptotic  expression 
el'^^  Although  the  farfield  approximation,  in  general,  is  not  always  valid  at  a 
distance  of  1  m,  this  treatment  is  standard  for  determining  transmitting 
sensitivities.  The  advantage  of  this  definition  is  that  the  pressure  field  may  be 
extrapolated  to  the  farfield  region  by  assuming  ordinary  spherical  spreading 
(which  is  generally  the  region  of  interest).  The  TCS  expressions  for  the  MHD 
source  mechanism  in  the  farfield  are  given  as 

TCS  =  p„(r  =  1  meter,  xy)/! 

where,  from  Eqs.(B.4)  and  (3.21),  the  TCS  can  be  stated  as 


2(op„kilj(b-a) 

(P„.,  (cose)  -  P„^,  (cose))  P„  (-cosy) 

dx^n(^)  x=kR 


(G^„(i,a)-Gx„(!l,b))X 


(P„.,  (cose)  -  P„^,  (COS0))  P„  (cos  V) 


d  .  1 ,  , 
zr  h„ (x) 

dx  n'  '  x=kR 


(4.22a) 


The  X  axis  transmitting  sensitivity  is  given  below  as  a  special  case  of 


Eq.  (4.22a). 


TCSq  s  pj,(r  =  1  meter,  \j/  =  0°  )/l 

jpq  _  ^0  P<m^(m 

2cop  kil  fb-a) 


-(Gx,_(-Jl,a)-GxL(-Ji,b))X 


(P„.,  (cosG)  -  P„^,  (COS0))  (-i)"e-P""*'>'2 


-f  h'(x) 
dx  n  '  x-kR 


V  (Pn  1  (COS©)  -  P„  ,  (cos©)) 

(G,,(i,a)  ■  G,,  (il,b))£  - Vf-i - 


;;r  h^(x) 

dx  n'  /  x»kR 


(4.22b) 


In  general  the  transmitting  sensitivity  is  the  most  important 
frequency  dependent  relation  for  a  transducer.  The  expressions  stated 
above  have  many  geometric  parameters  which  affect  the  transmitting 
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sensitivity  in  very  pronounced  ways.  It  is  not  the  intent  of  this  thesis  tc  probe 
into  a  detailed  parametric  investigation  of  the  sensitivity  functions  but  to  point 
out  a  few  important  points.  In  Figs.  4.4(a),  (b),  4.5(a),  (b),  and  4.6(a),  (b) 
plotted  in  three  dimensions  is  the  current  sensitivity  on  the  waveguide  axis 
(TCSq)  versus  frequency  versus  a  geometric  parameter  or  fluid  property  for  a 
specified  set  of  constant  parameters. 

The  variation  of  the  sensitivity  with  waveguide  sound  speed 
proved  to  be  of  practical  utility  when  analyzing  the  experimental  data.  When 
current  is  passed  through  the  NaCI  solution  in  the  transducer  it  is  inevitable 
that  some  hydrogen  and  chlorine  gas  bubbles  will  be  generated.  The  effect 
of  the  gases  is  to  lower  the  acoustic  phase  speed  due  to  the  large  increase 
in  the  small  signal  compressiblity  of  the  medium. 

The  transmitting  current  sensitivity  surface  shown  in  Fig,  4.4(a)  is 
for  the  case  of  a  0.03  m  source  centered  in  a  0.10  m  long  waveguide  with  a 
0.09  m  baffle  radius  (an  aperture  radius  of  0.0196  m  is  used  in  all  plots). 
Compared  to  Fig.  4.4(b),  which  is  for  a  source  of  0.10  m  with  the  same  baffle, 
the  shorter  source  seems  to  excite  a  sharper  resonance  due  to  the  source 
volume  being  spatially  a  better  fit  to  the  excited  "3/2  wave  length  pressure 
mode".  (See  Figs.  3.8,  3.9,  and  3.10  for  diagrams  of  the  particle  velocity  and 
pressure  modes.) 

At  higher  sound  speeds  a  null  forms  in  the  sensitivity  surface.  The 
null  corresponds  to  the  destructive  interference  of  the  acoustic  radiation  in 
the  farfield  and  is  primarily  dependent  on  the  radius  of  the  baffle.  Notice  that 
the  null  occurs  at  the  same  frequency,  approximately  6400  Hz.  in  both 


68 


FIGURE  4.5 

TRANSMITTING  CURRENT  SENSITIVITY  LEVEL 
AS  A  FUNCTION  OF  BAFFLE  RADIUS  AND  FREQUENCY 


(a) 

Cq  =  650  m/s 
R  =  0.09  m 


2  700  (b) 


Cq  I  1536  m/s 
R  =  0.09  m 
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FIGURE  4.6  scs-GA 

TRANSMITTING  CURRENT  SENSITIVITY  LEVEL  ^  ^7.85 

AS  A  FUNCTION  OF  SOURCE  LENGTH  AND  FREQUENCY 


figures,  which  indicates  it  is  not  a  feature  which  is  dependent  on  the  sound 
upeed  in  the  waveguide.  The  ha'^  circumference  of  the  baffle,  which  is  the 
separation  distance  of  the  apertures,  corresponds  to  3/4  of  a  wavelength  at 
6400  Hz  with  a  free  mcoium  sound  speed  of  1500  m/s.  Since  the  apertures 
radiate  180°  out  of  phase  for  MHD  sources  centered  symmetrically  about  the 
waveguide  length,  signal  cancellation  occurs  at  this  wavelength. 

The  transmitting  sensitivity  as  a  function  of  baffle  radius  is  shown 
in  Figs.  4.5(a)  and  (b),  where  the  periodicity  of  the  nulls  is  seen  to  be  a 
function  of  baffle  radius  for  the  reason  given  above. 

Another  important  feature  of  the  plots  of  Figs.  4.5(a)  and  (b)  is  that 
as  the  baffle  radius  decreases  the  overall  sensitivity  increases.  The  reason 
for  this  behavior  is  that  the  smaller  baffle  provides  a  smaller  acoustic 
impedance,  and  thus  a  larger  particle  velocity  at  the  aperture  results.  As  the 
size  of  the  baffle  increases,  the  acoustic  impedance  increases  and  the 
particle  velocity  is  diminished.  Note  that  the  increase  in  radiating  efficiency 
resulting  from  increasing  the  baffle  surface  area  does  not  increase  the 
radiated  pressure  amplitude  (sensitivity).  When  the  baffle  surface  area  is 
increased  the  aperture  impedance  increases  and  has  the  effect  of 
decreasing  the  particle  velocity  faster  than  the  improvement  in  the  baffle 
reflecting  property. 

The  higher  peak  amplitude  in  Fig.  4.4(a)  relative  to  Fig.  4.4(b)  is 
due  to  the  achievement  of  the  "3/2  wavelength"  resonance  at  a  lower 
frequency  due  to  the  slower  sound  speed  in  the  waveguide. 
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Shown  in  Figs.  4.6(a)  and  (b)  is  the  transmitting  sensitivity  as  a 
function  of  the  source  length,  (b  -  a)  or  Jig.  Again  it  is  lower  sound  speeds 
that  cause  the  greater  amplitudes  at  low  frequencies.  The  effect  of  source 
length  is  not  as  strong  for  the  higher  sound  speed  plot  in  Fig.  4.5(b)  since  a 
resonance  condition  has  not  been  reached.  The  primary  effect  in  this  case  is 
the  fixed  radius  of  t.ne  baffle  causing  the  resulting  destructive  interference 
null  structure  of  the  plot  at  the  same  frequency. 

E.  ELECTRICAL  IMPEDANCE  OF  THE  THERMOACOUSTIC 
TRANSDUCTION  PROCESS 

The  electrolyte  impedance  characteristics  and  transmitting 
pressure  sensitivity  of  the  ohmic  source  are  not  linear  relations,  that  is,  the 
electrical  current  is  not  linearly  related  to  the  voltage  and  the  transmitted 
pressure  is  not  linearly  related  to  the  input  voltage  (or  current).  The 
nonlinearity  of  the  ohmic  mechanism  is  evident  from  the  fact  that  an  electrical 
input  signal  of  a  single  frequency  co  produces  an  acoustic  output  of 
frequency  2  co,  as  demonstrated  in  Chapter  3.  The  nonlinear  behavior  of  the 
ohmic  mechanism  is  evident  in  the  source  term  of  the  inhomogeneous  wave 
equation,  -p  (J-J),/aCp  ,  where  the  dot  product  of  the  current  density  field 
with  itself  is  a  squaring  or  quadratic  nonlinearity.  One  obvious  consequence 
of  this  nonlinear  behavior  is  that  for  muhple  frequency  inputs  superposition 
of  the  single  frequency  Fourier  transform  solution  does  not  hold.  Therefore 
what  was  derived  in  Chapter  3  is  valid  for  the  single  frequency  input  co  and 
not  for  linear  combinations  of  this  single  frequency  solution,  it  is  important  to 
note  that  the  wave  equation  is  linear,  that  is,  the  spectral  content  of  the  wave 


fields  produced  by  the  nonlinear  source  do  not  undergo  frequency  domain 
corruption  due  to  propagation.  An  acoustic  signal  at  frequency  co  generated 
by  the  MHD  source  mechanism  will  be  propagated  with  frequency  co. 
Similarly  a  single  frequency  electrical  input  signal  of  frequency  co,  through 
the  ohmic  mechanism,  will  generate  a  wavefield  of  2co,  but  the  medium  will 
transmit  the  signal  uncorrupted  at  frequency  2(o.  However,  the  boundary 
conditions  are  frequency  dependent  and  therefore  the  wavefields  are 
frequency  "dependent",  and  thus  the  spectral  content  of  the  waves  emitted 
by  the  source  must  be  known  in  order  to  determine  the  wavefield. 

Fortunately  the  source  nonlinearity  can  be  Fourier  transformed  by 
using  convolution  and  thus  the  spectral  output  of  the  thermoacoustic  source 
can  be  determined.  If  the  current  density  has  the  functional  form, 

J{x.t)  =  J(x)  f(t)  .  (4.23) 

then  the  nonhomogeneous  ohmic  source  term  takes  the  form, 

-li 

—  (J(x)-J(x))  f(t)  f(t) ,  (4.24) 

oCp 

where  the  dot  product  operation  of  the  spatial  component  of  J  was  noted  in 
Chapter  2.  The  Fourier  transform  can  now  be  performed  on  f(t)  f(t)  by 
convolution  in  the  frequency  domain  as  foiiows. 

Q(a))  =  f{f(t)  f(t)}  =  J  F(a)')F(o)  -  co')  do)’,  (4.25) 

-oo 

where  Q  is  the  spectrum  of  the  convolution,  F  is  defined  as 

F(a))sf(f(t)},  (4  26c.; 
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and  -f  is  the  exponential  Fourier  transform  operator  defined  as 

f{ }  -=  f{ }  ej“'’dt.  (4.26b) 

-oo 

As  a  note,  it  should  be  kept  in  mind  that  when  the  multiplication  of  time 
harmonic  functions  is  performed,  f(t)  f(t),  the  e  i®’’  "notation"  must  be  replaced 
by  the  real-time  functions  such  as  cos((ot).  The  reason  is  simply  that 
multiplication  is  a  nonlinear  operation  and  when  it  is  performed  on  the 
complex  time  functions  cor  mcniy  used  for  notational  ease  incorrect  or 
incomplete  harmonic  components  are  produced  by  the  convolution 
operation.  As  an  example  of  the  convolution  operation  consider  the  following 
time  dependence  for  J,  assuming  co^  >  cog, 

f(t)  =  A^  cos(a)^t)  +  Aj  cos(co2t) .  (4.27) 

The  Fourier  transform  of  Eq.  (4.27)  is 

F(a))  =  A,(8{a)+o),)  +  5(a)-co^))  +  A2  •  (4.28) 

The  convolution  results  in  the  following  source  spectrum,  Q, 

Q(co)  =  A,2(5(co+2a)^)  +  5(a)-2a3^))  +  Ag^  (5(o)+2o)2)  5(a)-2a)2)) 

+  2A,  A2  (5(o)+a)^+o)2)  +  5(o)-(o^-o)2)  +  5((o+(o^-o)2)  +  5((o-o)^+o)2)) , 

(4.29) 

which  indicates  that  the  following  wavefield  harmonics  are  produced;  2(Op 
0),+  0)2,  20)2,  and  0)2- 

Notice  that  the  output  of  the  simple  two-frequency  input  is  the 
generation  of  four  different  harmonics.  For  this  two-frequency  input  the  single 
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frequency  solutions  of  Chapter  3  can  now  be  written  in  terms  of  the  above 
fiequency  components  and  the  resultant  pressure  fields  can  be  summed 
(since  the  acoustic  equation  is  linear)  to  obtain  the  total  acoustic  field. 

Due  to  the  nonlinearity  of  the  source  mechanism,  the  linear 
impedance  is  not  defined.  What  can  be  defined  instead  is  the  quadratic 
current  transfer  function  (QCTF).  The  QCTF  is  defined  as  the  ratio  of  the 
output  current  to  the  square  of  tne  input  voltage  for  the  input  of  two  frequency 
components  as  follows.-'^ 


QCTFh 


VaK) 

The  expression  for  the  ohmic  source  is  as  follows. 

II  oBQu(co^4-tog,  x)  dx  dz 


(4.30) 


QCTF 


(4.31) 


I  -f 

■p  B  (R. +et2ka .  )+  R  +  2e’i'^<®+‘^))) , 

8.) 


(4.32) 


where  k  =  an  example  of  the  eva'uation  of  Eq.  (4.32),  the 

QCTF  for  difference  frequency,  o.-cog,  would  be  found  by  substituting  (0y^=  co^ 
and  (Og  =  -(Og  for  the  case  of  a  single  frequency  input  (o,  co^  =  co,  and  Wg  =  co. 

The  QCTF  should  be  regarded  as  a  term  of  second  order  for  NaCI 
solutions,  due  to  the  combination  of  low  thermal  exparsion,  conductivity,  and 
magnetic  field.  The  QCTF  for  a  6%  NaCI  solution  is  shown  in  Fig.  4.7(a)  and 
(b)  for  a  single  frequency  input  signal.  The  amplitudes  of  the  plots  are 
insignificant  at  any  reaiistic  voltage  levels. 
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For  liquids  with  a  very  high  conductivity  (such  as  mercury)  the 
QCTF  will  not  be  of  second  order.  Plots  for  mercury  in  the  same  waveguide 
configuration  as  the  NaCI  show  much  higher  QCTF  magnitudes  (see 
Fig.  4.8(a)  and  (b)). 

Since  the  QCTF  is  the  current  relationship  the  total  complex 
current,  ij,  (MHD  and  thermoacoustic  effects)  is 

i,  AV  +  QCTFV2  ,  (4.33) 

where  A  is  the  admittance  due  to  the  MHD  mechanism  and  is  equal  to  1/Z. 
The  equivalent  circurit  of  the  combined  MHD  and  thermoacoustic  impedance 
including  the  polarization  impedance  is  shown  in  Fig.  4.9. 

F.  THERMOACOUSTIC  TRANSDUCTION  PROCESS  POWER  EFFICIENCY 
The  thermoacoustic  t.'ansduction  power  efficiency  is  defined  as 

the  radiated  acoustic  power  due  to  the  thermoacoustic  mechanism  divided 
by  the  real  input  electrical  power  to  the  MHD  transducer.  The  complex  input 
power  is  defined  in  Eq.  (4.18),  but  if  the  QCTF  is  significant  the  effect  needs 
to  be  accounted  for  in  the  electrical  power  expression.  The  total  complex 
electrical  power,  Wj^,  is  then  given  by 

WeTH  =  ViT.  (4.34) 

The  radiated  acoustic  power  is  defined  using  Eq.  (4.20)  with  the 
particle  velocity  expressions  of  Eqs.  (3.32a,c)  being  substituted,  and  the 
same  expressions  for  the  real  component  of  the  self-impedance  being  valid. 
With  Eq.  (4.20)  and  Eq.  (4.34)  the  thermoacouctic  process  efficiency  can  be 
defined  as 
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FIGURE  4.9 

EQUIVALENT  CIRCUIT  OF  THE  COMBINED  MHD 
AND  THERMOACOUSTIC  INPUT  IMPEDANCE 
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Plots  of  the  thermoacoustic  transduction  process  efficiency  are 
given  in  Fig.  4.10(a)  and  (b)  for  a  6%  NaCl-water  solution  and  mercury, 
respectively.  Comparing  the  MHD  efficencies  given  in  Figs.  4.3(a)  and  (b)  to 
the  thermoacoustic  efficiencies  in  Fig.  4.10(a)  and  (b)  the  MHD  process  is 
much  more  efficient  for  both  fluids. 

A  rough  relationship  between  the  fluid  properties  of  the  medium 
and  the  thermoacoustic  process  efficiency  is  obtained  from  the  leading 
coefficients  of  the  ratio  of  acoustic  power  expression  Eq.  (4.20)  to  the 
expression  for  the  real  dc  ohmic  power  dissipation  (the  primary  mechani.sm 
of  power  dissipation  in  the  experimental  transducer), 


JoV  C,3 


coCj  apH. 


(4.36) 


The  same  type  of  expression  can  be  written  for  the  MHD  process  and  is 
stated  as 


^MHD  = 


(4.37) 


Some  interesting  insights  are  gained  by  these  two  equations. 
First,  the  thermoacoustic  efficiency  is  a  function  of  the  current  density 
squared,  or  power,  and  thus  makes  comparison  of  efficiencies  a  matter  of 
first  stating  the  current  density.  Second  a  low  conductivity,  mass  density, 
specific  heat,  and  high  coefficient  of  thermal  expansion  are  desirable  for  an 
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efficient  thermoacoustic  source.  For  the  MHD  source  the  opposite  is  true  with 
regard  to  \he  conductivity;  a  high  value  is  desired,  so  that  little  ohmic  losses 
are  generated.  As  the  plots  in  Fig.  4,3(a)  and  (b)  show,  mercury  is  a  more 
efficient  fluid  for  the  MHD  process  than  the  NaCI  solution  by  a  factor  of 
60  dB,  due  to  mercury's  1,0  X  10®  1/(ohm  m)  conductivity  versus  the  5.26 
1/(ohm  m)  conductivity  of  the  NaCI  solution.  The  opposite  is  of  course  true  for 
the  thermoacoustic  case  where  the  NaCI  solution  is  approximately  70  dB 
more  efficient.  However,  mercury  does  have  a  value  of  CpOf  140  J/kg/C 
versus  4100  J/kg/C  for  the  NaCI  solution,  but  the  mass  density  of  mercury  is 
13600  kg/m®  versus  1041  for  the  NaCI  water  solution. 

G.  TRANSMITTING  SENSITIVITY  OF  THE  THERMCACOUSTIC  PROCESS 
The  transmitting  sensitivity  of  a  transducer  is  a  transfer  function, 

like  the  admittance  or  impedance,  and  is  generally  defined  in  the  frequency 
domain  for  a  single  input  frequency  assuming  that  the  same  output 
frequency  is  produced.  The  transmitting  sensitivities  of  the  thermoacoustic 
mechanism  cannot  be  defined  using  the  single  frequency  definition  since  the 
thermal  source  mechanism  is  nonlinear  and  produces  spectral  pressure 
signal  components  different  from  the  input  electrical  signal.  A  quadratic 
transfer  function  will  be  defined  for  two  input  frequencies  as  the  quadratic 
voltage  transmitting  sensitivity  (QVTS), 

P«(o)a+cOr,  r  =  1  meter,  \j;  =  main  beam  axis) 

QVTS  =  .  -  ,  (4.38) 

VAK)VB(t%) 

where  V^(a);^)  and  Vq(cob)  are  the  Fourier  transformed  expressions  for  the 
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voltage  amplitudes  at  co;^  and  ©g.  Substituting  Eq.  (8.4)  for  Pjj,  Eq.  (3.32b)  for 
the  particle  velocity  distribution,  and  Eq.  (4.1)  for  the  potential  V  into 
Eq.  (4.38),  the  analytical  expression  for  the  QTVS  results, 


QVTS  = 


BpgPfmCfm 


2k{1  -  R,R  ,) 


jkfi. 


) 


(P„.^  (COS0)  -  (COS0))  P„(-cos\i/) 


n=o 


’“''lin'’ 


((e+> .  e+'>) .  R  (gik^ .  (6'"“ •  R,  el^) 


2k(l-R,R,) 


(P„.,  (COS0)  -  P„^,  (COS0))  P„  (cosy)  e'f 


n=o 


x-kf^R 


(4.39) 


where  kj^  =  (w^  +  co0)/c,f„  and  k  =  (©^^  +  W0)/Co:  for  a  single  frequency  input  ©, 

'<fm  =  2o)/c,^andk  =  2©/Co. 

Time  did  not  permit  a  parameter  investigation  such  as  that 
performed  for  the  the  MHD  sensitivity.  K  wever,  two  plots  of  the  single 
frequency  input  thermoacoustic  transmitting  sensitivities  (QVTS)  for  fixed 
parameters  are  given  in  Fig.  4.11(a)  and  (b)  for  the  6%  NaCI  solution  and 
mercury.  Notice  that  the  mercury  transducer  shows  a  higher  sensitivity  than 
the  NaCl  solution.  However,  the  sensitivity  is  in  terms  of  the  input  voltage. 
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With  mercury's  very  high  conductivity  a  very  small  voltage  is  always  going  to 
result,  at  the  terminals,  and  thus  a  low  source  level.  Thi.'>  is  because 
generating  large  voltages  on  a  practical  thermoacoustic  merely*  transducer 
is  tantamount  to  huge  electrical  currents  due  to  the  real  input 
impedance. 

More  discussion  of  the  QVTS  is  included  in  C;  apter  5  with  the 
experimental  results. 
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CHAPTER  5 

TRANSMITTING  EXPERIMENT 


To  verify  the  theoretical  predictions  presented  in  Chapters  2,  3, 
and  4,  an  MHD  transmitter  was  constructed  and  sensitivity  and  beam  pattern 
measurements  were  taken.  The  transducer  used  in  the  experiments  is 
depicted  in  the  photograph  of  Fig.  5.1. 

The  planned  operating  regime  of  the  the  MHD  transmitter  with 
respect  to  the  frequency  bandwidth  and  output  level  was  primarily  dictated 
by  the  characteristics  of  the  power  amplifier  that  was  available  and  the 
spatial  extent  of  the  magnetic  field  that  could  be  generated.  Fortunately  a 
high  power,  relatively  broadband  power  amplifier  was  available  at  ARL:UT. 
The  magnetic  field  was  generated  using  samarium-cobalt  permanent 
magnets.  Permanent  magnets  were  employed  primarily  because  of  the 
convenience  of  not  needing  a  power  supply,  but  due  to  the  expense  of  the 
magnets  a  limited  volume  of  low  permeable  space  could  be  energized 
uniformly.  An  important  consideration  was  that  the  permanent  magnet  be 
relatively  compact  and  of  simple  design  since  the  size  of  the  transducer  has 
an  important  effect  on  the  acoustic  radiating  characteristics. 

A.  OPTIMIZATION  OF  THE  MHD  TRANSMITTING  CURRENT  SENSITIVITY 

The  prime  objective  in  designing  the  MHD  transducer  was  to 

maximize  the  transmitting  sensitivity  across  a  broad  frequency  bandwidth. 
Tantamount  to  this  objective  is  determining  the  waveguide  dimensions,  the 
permanent  magnet  spatial  distribution,  and  the  electrode  surface  distribution. 
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FIGURE  5.1 

MAGNETOHYDRODYNAMIC  TRANSDUCER 
APERTURE  VIEW 
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From  the  transmitting  current  sensitivity  expression,  Eq.  (4.22b),  it  is  obvious 
that  a  large  magnetic  induction  is  desirable,  but  maximizing  the  magnetic 
induction  does  not  optimize  the  design  goal  using  permanent  magnets. 
Since  the  radiating  aperture  area  is  an  important  parameter,  it  would  seem 
that  a  large  aperture  together  with  a  high  magnetic  field  would  contribute  to 
maximizing  the  sensitivity.  However,  the  size  of  the  aperture  is  directly 
related  to  the  pole  face  separation  of  the  magnets,  and  therefore  increasing 
the  aperture  area  increases  the  reluctance  of  the  magnetic  circuit.  There  is  a 
definite  optimization  of  the  magnetic  induction  by  the  aperture  size  which  will 
maximize  the  transmitting  sensitivity  for  a  given  quantity  and  spatial 
arrangement  of  permanent  magnets. 

The  functional  dependence  of  the  sensitivity  on  the  aperture 
radius  and  baffle  radius  is  not  explicit  in  Eq.  (4.14b):  the  uepsndence  arises 
in  the  computation  of  the  reflection  coefficients.  The  relationship  between 
sensitivity  and  aperture  size  is  displayed  in  Fig.  5.2  for  a  spherical  baffle 
radius  of  0.05  m  and  frequency  of  1000  Hz.  The  sensitivity  function,  for  the 
purpose  of  the  optimization,  is  well  fit  by  a  linear  relation  in  the  region  from  0 
to  29  mm  (the  general  trend  is  for  the  sensitivity  to  increase  over  the  entire 
frequency  band  as  seen  in  Fig.  5.3).  At  1000  Hz  and  for  a  baffle  of  radius 
0.05  m  the  relation  is 

TCSoc  727.3  Rg,  (5.1) 

where  R3  has  the  units  of  meters. 

It  is  also  necessary  to  relate  the  magnetic  induction  to  the 

aperture  radius  by  the  appropriate  magnetic  circuit  model.  The  analysis  of 
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the  magnetic  circuit  follows  from  Cousins  and  Nash^^  (they  constructed  a 
cloud  chamber  using  permanent  magnets).  The  important  technical  aspect  of 
their  work  is  the  development  of  an  empirical  relation  tor  II'.c  magnetic  flux 
leakage  (fringing  losses)  associated  with  square  or  rectangular  magnetic 
pole  faces  with  various  geometric  restrictions.  Due  to  the  manufactured 
shape  of  the  permanent  magnets  a  limited  number  of  practical  geometric 
configurations  were  possible.  Since  the  research  budget  allowed  the 
purchase  of  40  magnets  of  dimensions  2"  x  2"  x  0.5"  an  arrangement  of  four 
stacks  of  five  magnetics  above  and  below  the  field  gap  was  intuitively 
obvious.  This  was  not,  however,  an  ad  hoc  choice;  the  circuit  model 
predicited  that  larger  but  fewer  stacks  would  result  in  large  fringing  losses. 
Therefore,  given  the  four-stack  arrangement,  the  relationship  between  the 
magnetic  induction  and  aperture  radius  was  found  to  be 

B  = -1 5350  R3  + 8250,  (5.2) 

where  the  units  of  B  are  tesla.  The  com|.uter  generated  curve  is  shown  in 
Fig.  5.4. 

Stated  previously,  the  sensitivity  is  directly  proprortional  to  B,  so 
Eq.  (5.1)  and  Eq.  (5.2)  can  be  combined  to  yield  an  empirical  relation  for  the 
TCS  as  a  function  of  R^, 

TCS-  (-15.35  Ra  + 0.8250)  1.106  Rg.  (5.3) 

To  maximize  the  relationship  in  Eq.  (5.3)  the  expression  is  differentiated  with 
respect  to  Rg  and  the  result  set  equal  to  zero.  The  value  for  Rg  is  found  to  be 
0.027  m,  which  is  a  field  gap,  Jl^,  of  0.0479  m.  Due  to  an  error  in  the 
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MAGNETIC  INDUCTION,  B  -  GAUSS 


calculations  made  at  the  time  the  MHD  transducer  was  designed  a  value  of 
0.0381  was  used  for  the  field  gap.  The  construction  details  of  the  waveguide 
are  shown  in  Fig.  5.5. 

The  estimate  of  the  magnetic  induction  for  the  0.0479  m  field  gap 
is  4500  G  and  maximum  measured  value  for  the  0.0381  m  gap,  displayed  in 
the  distribution  shown  in  Fig.  5.6,  is  4327  G.  It  was  concluded  that  either  the 
assumptions  concerning  the  geometricVi  restrictions  on  the  validity  of  the 
fringing  theory  are  in  slight  viola.icn,  or  small  gaps  in  the  transducer  magnet 
stack  raised  the  circuit  reluctance  higher  than  expected.  However,  the 
magnetic  induction  is  very  uniform  as  shown  in  Fig  5.6  for  the  distribution 
measured  at  the  midpoint  plane  of  the  waveguide. 

B.  ACOUSTIC  SPEED  MEASUREMENT  OF  THE  WAVEGUIDE 

The  phase  or  plane  wave  acoustic  speed  for  the  6%  NaCI-water 

solution  in  the  MHD  waveguide  was  measured  to  confirm  the  boundary 
condition  assumption  that  the  waveguide  walls  were  acoustically  rigid.  Two 
methods  were  used  to  measure  the  wave  speed,  (1)  a  pulse  and  (2)  a 
resonance  technique. 

The  pulse  technique  used  a  transmit  frequency  of  150  kHz  and  a 
time  duration  of  three  cycles  or  10  ps.  This  method  yielded  a  wave  speed  of 
1550  m/s. 

The  rcoonance  method  employed  a  conti.ious  wave  white  noise 
signal  transmitted  into  the  waveguide  with  the  apertures  terminated  by  "rigid" 
bra.ss  blocks.  The  boundary  conditions  were  thus  taken  to  be  rigid  walled 
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cavity  (see  Fig.  5.7).  By  measuring  the  amplitude  spectrum  of  the  cavity 
response  to  the  white  noise  signal  the  frequencies  of  the  standing  wave 
modes  were  detet  mined  and  used  to  calculate  the  wave  speed.  From  the 
1 ,0,0  pressure  mode  of  the  cavity  the  wave  speed  was  calculated  to  be 
1600  m/s.  For  most  of  the  theoretical  calculations  1600  m/s  was  used 
because  the  sound  speed  was  determined  by  the  resonance  method  at 
7500  Hz,  which  is  within  the  frequency  band  of  the  experiment.  The 
theoretical  sound  speed  for  the  6%  NaCI-water  solution  at  20°  C  was 
calculated  to  be  1536  m/s  by  using  of  the  equation  given  by  Coppens^®  for 
the  speed  of  sound  in  seawater  as  a  function  of  salinity  and  temperature. 

C.  MEASUREMENT  OF  THE  ELECTRICAL  IMPEDANCE  OF  THE 
TRANSDUCER 

The  input  impedance  of  the  MHD  transducer  was  measured 
using  the  system  of  Fig.  5.8  with  the  transducer  in  the  water.  The  input  signal 
to  the  transducer  was  white  noise  at  an  rms  amplitude  of  20  V  and  the 
current  was  sensed  using  a  Pearson  current  transformer.  The  transfer 
function  mode  of  the  signal  analyzer  was  used  to  determine  the  impedance 
using  the  voltage  as  the  response  signal  to  the  input  current  signal.  The 
transfer  function  is  calculated  as  the  crosscorrelation  between  the  input  and 
output  complex  signal  spectra  divided  by  the  autocorrelation  of  the  complex 
input  signal  spectrum.  A  plot  of  the  magnitude  and  phase  of  the  impedance 
measurement  is  shown  in  Fig.  5.9. 

The  measured  impedance  indicates,  by  the  lagging  phase,  a 
significant  capacitive  impedance  at  frequencies  less  than  4000  Hz.  This 
capacitive  effect  is  due  to  the  polarization  impedance  of  the  electrodes. 
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which  was  investigated  theoretically  in  Section  A  of  Chapter  4.  Plotted  in 
Fig.  5.10  are  the  real  and  imaginary  components  of  the  theoretical 
impedance  represented  by  the  network  shown  in  Fig  4.1  with  data  points 
taken  from  Figs.  5.9.  Notice  the  ordinate  axis  is  the  reciprocal  of  the  square 
root  of  the  angular  frequency  so  that  the  linear  dependence  of  the  diffusional 
impedance  on  1/Vco  can  be  exhibited  (see  Eqs.  (4.11)  and  (4.12)).  One  can 
clearly  see  the  linear  trend  in  the  data  at  low  frequencies.  It  is  also  noted 
that  the  best  fit  of  the  theoretical  curves  occurred  for  a  double  layer 
capacitance  value  of  762  |is  which  is  200  mF/m^,  and  an  electrolye 
conductivity  of  5.56  1/ohm  m,  which  is  a  dc  resistance  of  1.5  ohms,  was 
included  in  the  calculation  of  Z  (Eq.  (4.17)).  The  imaginary  component  of  the 
MHD  impedance  was  found  to  have  a  negligible  effect  on  the  input 
impedance  model  and  can  therefore  be  discarded  in  most  preliminary 
calculations  using  a  low  conductivity  fluid. 

It  is  concluded  that  the  polarization  impedance  model  fits  the  data 
well;  however,  the  model  would  become  extremely  difficult  to  apply  if 
electrolytes  other  than  1 ,1  valent  salts  are  used.  This  might  be  the  case  with 
seawater.  Vetter  gives  a  very  complete  discussion  on  the  impedance 
properties  of  other  common  electrolyte  chemistries^  and  some  of  the 
difficulties  of  complex  chemistries. 

D.  TRANSMITTING  CIRCUIT  AND  ACOUSTIC  SIGNAL  MEASUREMENT 
SYSTEM 

The  efficiency  of  both  the  MHD  and  thermoacoustic  mechanism  is 
extremely  low  in  the  case  of  low  conductivity  fluids  as  evidenced  by 
Fig.  4.3b.  As  stated  previously  the  efficiency  is  better  at  low  frequency 
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since  the  particle  velocity  is  larger.  However  an  efficieny  of  -108  dB  implies 
that  an  input  of  5  kW  of  electrical  power  will  produce  only  20  mW  of  radiated 
acoustic  power.  The  problem  of  performing  the  experiment  is  whether  the 
signal-to-noise  ratio  (S/N)  is  sufficient  to  produce  meaningful  data. 

The  power  amplifier  used  to  power  the  trancducer  was  a  CML 
Corp.  model  B5K  which  can  deliver  5000  V-A  loaded  with  power  factors 
between  1  and  0.1  over  the  band  from,  100  to  10000  Hz,  ±0.5  dB.  Since  the 
maximum  lag  in  the  phase  of  the  impedance  is  -25°  (Fig.  5.11),  which 
corresponds  to  a  power  factor  of  0.91,  the  compensation  of  the  transducer 
was  not  deemed  necessary.  The  block  diagram  of  the  transducer  circuit  is 
shown  in  Fig.  5.11. 

An  impedance  matching  transformer  was  therefore  designed  to 
match  the  dc  resistance  of  1 .5  ohms  of  the  transducer  to  the  amplifier.  The 
transformer  turns  ratio  is  1.43,  with  38.6  turns  on  the  primary  and  30.6  turns 
on  the  secondary.  The  matching  transformer  has  a  2  mil  silicon  steel 
laminate  core  capable  of  broadband  power  transfer  with  a  maximum 
hystersis  loss  of  10  W  at  400  Hz.  The  total  loss  of  the  transmission  lines  and 
transformer  was  measured  at  300  W.  For  all  the  transmitting  experiments  the 
transducer  current  was  typically  60  A  peak  with  a  terminal  voltage  of  75  V 
peak. 

The  experiment  was  performed  at  the  large  outdoor  wooden  tank 
facility,  as  depicted  in  Fig.  5.1 1.  A  broadband  noise  survey  was  performed  to 
access  an  S/N  for  the  experiment.  Both  the  rms  level  and  the  coherency 
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between  two  hydrophones  1  m  apart  were  measured  and  averaged  over  a 
period  of  500  ms.  Over  the  band  from  500  to  3500  Hz  the  rms  noise  SPL  was 
approximately  97  dB  (re  1  jis)  and  then  proceeded  to  drop  off  at 
20  dB/decade.  The  coherency  was  measured  to  be  iess  than  0.4  over  much 
of  the  band  from  500  to  10000  Hz,  which  suggests  that  coherent  ping  to  ping 
averaging  of  the  received  acoustic  signal  should  be  a  successful  method  for 
eliminating  the  narrowband  noise  after  filtering. 

The  calculated  SPL,  on  axis,  at  15  ft  range  for  5000  W  of  input 
power  at  500  Hz  is  89  dB  which  gives  an  S/N  of  -8  dB,  this  represents  the 
worst  S/N  figure.  At  5000  Hz  the  SPL  is  123  dB  and  the  S/N  is  12  dB.  Signal 
averaging  thus  became  a  mandatory  process.  Typically  ten  signal  averages 
were  required  to  reduce  the  noise  by  20  dB.  And  typically  a  center  frequency 
to  filter  band  ratio  of  10  was  maintained  with  the  bandpass  filter  to  eliminate 
broadband  signal  noise. 

The  acoustic  signal  measurement  system  consists  of  a  broad¬ 
band  hydrophone  (NRL  USRD  H56  reference  standard  hydrophone),  signal 
amplifier,  bandpass  filter,  and  a  Nicolet  4094  digital  oscilloscope  (the  Nicolet 
has  a  signal  averaging  mode).  The  signal  amplifier  provided  40  dB  of  gain 
necessary  for  lifting  the  signal  above  the  elect,  ical  noise  produced  by  the 
transmission  cable's  150  ft  run  from  the  tank  to  the  building. 

The  MHD  transducer  was  mounted  on  a  rotator  column,  depicted 
in  Fig.  5.12,  and  lowered  to  a  depth  of  15  ft.  The  MHD  transducer  and 
hydrophone  were  separated  by  a  distance  of  18  ft,  which  placed  the 
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hydrophone  well  into  the  farfield  over  the  entire  trequency  band,  and 
allowed  a  maximum  gated  transmit  signal  of  5  ms;  the  gate  time  was  limited 
by  the  surface  reflection  path  length  of  the  transmitted  acoustic  signal.  The 
5  ms  gate  time  set  a  lower  limit  on  the  frequency  band  of  the  transmitted 
signal  to  500  Hz,  which  is  2.5  cycles. 

E.  MEASUREMENT  OF  THE  MHD  TRANSMITTING  SENSITIVITY  AND 
DIRECTIVITY 

The  TCS  on  the  waveguide  axis  was  measured  for  the 
frequency  band  from  500  Hz  to  13000  Hz.  A  typical  time  domain  pressure 
signal  received  while  making  the  sensitivity  measurement,  along  with  the 
current  waveform,  are  shown  in  Fig.  5.13.  The  acoustic  signal  shown  is  the 
result  of  ten  ping-to-ping  averages  and  bandpass  filtering. 

The  transmitting  current  sensitivity  is  defined  in  this  thesis  as  20 
log  10  0^  pressure  divided  by  the  peak  current.  Figure  5.14  shows 

the  comparison  of  measured  current  sensitivity  and  the  theoretical  prediction 
from  Eq.  (4.22b).  Good  agreement  was  found  between  theory  and 
experiment  for  the  500  Hz  to  4000  Hz  band  with  most  measured  values 
within  2  dB  of  the  theory.  Erratic  behavior  of  the  measured  data  for 
frequencies  greater  than  4000  Hz  can  probably  be  attributed  to  two  causes. 
First,  two  restrictions  were  placed  on  the  theory  on  Chapter  3,  R3/Rg«1  and 
kR3«1.  The  first  constaint  is  set  by  the  geometry  to  be  effectively  0.392, 
which  therefore  might  raise  some  doubt  about  the  transducer  design. 
However,  the  good  fit  of  the  data  at  less  than  4000  Hz  seems  to  relax  the 
importance  of  this  constaint.  At  4000  Hz,  kR3=  0.33  (which  again  is  not  much 
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less  than  1),  may  indicate  that  the  aperture  radius  to  wavelength  restriction  is 
more  severe  than  the  constraint  of  planar  motion  of  the  particle  velocity  at  the 
aperture. 

The  second  cause  is  the  non-spherical  construction  of  the 
transducer  (see  Fig  5.1).  Since  the  measured  data  do  not  diverge  smoothly 
from  the  theorv  one  would  expect  that  diffraction  caused  by  the  irregular 
geometry  of  the  magnetic  r’.uuit  structure  is  at  fault.  At  4000  Hz  the  wave 
length  is  37.5  cm  and  the  diameter  of  the  ring  structure  is  29.2  cm,  so 
diffraction  should  be  expected.  Given  the  geometrical  restrictions  of  the 
theory  and  good  fit  of  data  well  beyond  those  restrictions  one  can  conclude 
that  the  theoretical  approach  is  appropriate. 

Shown  in  Figs.  5.15(a),  (b),  (c),  and  (d)  are  the  measured  and 
predicted  directivity  patterns  for  1000,  4000,  6000,  and  10000  Hz  (the 
theoretical  directivity  patterns  were  computed  by  the  FORTRAN  program 
MHD).  The  th'  oretical  directivities  are  shown  for  two  different  baffle  radii  to 
emphasize  the  importance  of  the  baffle  size  in  determining  the  directivity 
pattern.  In  Fig.  5. 16  tfie  equivalent  baffle  radius  is  compared  to  the  actual 
transducer  geometry. 

The  directivity  patterns  were  measured  at  10°  intervals  for  a 
complete  360°  rotation  of  the  MHD  transducer.  At  1000  Hz  there  is  excellent 
data  agreement  with  theory.  At  4000  Hz  there  is  some  significant  divergence, 
and  at  6000  and  10000  Hz  the  data  and  theory  do  not  correlate,  except  near 
the  90°  and  270°  null.  The  directivity  pattern  at  4C00  Hz,  and  higher 
frequencies,  indicates  that  diffraction  by  other  parts  of  the  transducer  is 
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significant  and  is  thus  the  cause  of  the  divergence  of  the  on-axis  sensitivity 
data  above  4000  Hz. 

F.  MEASUREMENT  OF  THE  THERMOACOUSTIC  TRANSMITTING 
SENSITIVITY 

Measurement  of  the  QVTS  was  made  on  the  waveguide  axis  for 
the  single  input  frequency,  in  the  band  between  500  and  4000  Hz,  which 
corresponds  to  the  transmitted  acoustic  frequency  doubled  band  of  1000  to 
8000  Hz.  In  order  to  measure  the  thermoacoustic  signal  it  was  necessary  to 
remove  the  waveguide  from  the  magnetic  ring  structure,  thus  eliminating  the 
MHD  signal.  Accurate  measurement  of  the  thermoacoustic  signal  in  the 
presence  of  the  MHD  signal  proved  to  be  difficult  because  of  MHD  signal 
interference.  Even  though  the  bandpass  filter  rejected  40  dB  of  the  MHD 
signal  amplitude,  spectral  leakage  was  a  problem. 

A  typical  thermoacoustic  time  waveform  is  shown  in  Fig.  5.17. 
Notice  the  frequency  doubled  pressure  signal  relative  to  the  input  current 
signal.  Noise  was  much  more  of  a  problem  in  measuring  the  thermoacoustic 
signal  since  the  pressure  signal  is  typically  10  dB  below  the  levels  measured 
for  the  MHD  process. 

The  measured  QVTS  data  and  the  theory  are  plotted  in  Fig.  5.18 
as  a  function  of  the  input  signal  frequency.  The  QVTS  in  this  plot  is  defined 
as  the  peak  pressure  divided  by  the  square  of  the  peak  terminal  voltage.  In 
this  case  the  quadratic  current  transmitting  sensitivity  is  related  to  the  QVTS 
by  substacting  7  dB  from  the  latter. 
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The  number  of  data  points  shown  in  Fig.  5.18  are  limited  due  to 
difficulties  with  the  transducer  fluid  chemistry.  The  thermoacoustic  data  were 
gathered  several  months  after  the  experiment  was  set  up,  and  during  this 
time  the  electrode  surfaces  had  begun  to  corroJe.  The  electrodes  are  copper 
with  nickel  and  gold  electroplated  on  the  surface.  At  the  time  of  the 
thermoacoustic  experiment  the  copper  had  apparently  reacted  with  the 
nickel  and  blistered  the  gold  plate  off  the  surface,  exposing  the  now  oxidized 
copper  surface.  As  a  result  large  quantities  of  undissolved  gases  were 
generated  due  to  a  copper  oxide-chlorine  reaction.  After  virtually  each  ping 
the  transducer  had  to  be  raised  from  the  tank  and  checked  for  visible  gas 
Dubbles,  and  the  sodium  chloride  and  water  solution  changed.  Thus  few 
reliable  data  points  were  generated.  However,  the  few  points  that  were 
obtained  generally  fit  within  2  dB. 

As  a  note  to  the  problem  of  hydrogen  and  chlorine  gas  evolution 
within  the  transducer,  the  gases  seemed  to  stay  dissolved  until  the  saturation 
concentration  was  reached  for  CIg,  which  is  0.0004  moles/liter.  In  order  to 
produce  this  concentration  in  the  waveguide  a  total  charge  transfer  of  93°C 
is  required.  At  60  A-peak  current  the  average  current  is  38.2  A;  thus  1 .54  s  of 
current  is  required  to  produce  the  saturation  concentration  of  chlorine  gas. 
This  is  equivalent  to  309  5  ms  pings,  which  is  enough  data  to  generate  the 
data  in  Fig.  5.14.  However,  the  transducer  must  be  periodically  shaken  to 
dissolve  small  nucleated  gas  bubbles  which  can  sometimes  form  on  the 
surface  of  the  electrode.  The  curious  inconsistency  in  the  production  of 
hydrogen  and  chlorine  is  that  it  does  not  seem  to  react  to  form  hydrochloric 
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acid  according  to  Eq.  (4.6)  and  thus  neutralize  the  buildup  of  sodium 
hydroxide  that  results  according  to  the  overall  electrode  reaction,  Eq.  (4.5). 
The  dissolved  chlorine  gas  seems  to  reach  saturation  when  bubbles  are 
noticed  in  the  waveguide  after  approximately  250  pings. 


CHAPTER  6 


CONCLUSIONS 

This  thesis  has  examined  the  behavior  of  an  underwater  acoustic 
transmitter  based  on  the  magnetohydrodynamic  principle.  The  transducer 
consists  of  a  waveguide  filled  with  salt  water  and  exposed  to  orthogonal 
magnetic  and  electric  fields.  The  magnetic  field  was  generated  by 
permanent  magnets,  and  the  electric  field  was  time  harmonic  in  order  to 
produce  oscillatory  force  on  the  salt  water  causing  an  acoustic  signal.  This 
electric  field  results  in  electric  current  flow  which  causes  heating  and 
thermoacoustic  radiation  of  sound  as  well  as  producing  sound  through  the 
MHD  mechanism.  Both  sound  generating  mechanisms  were  examined. 

The  linear  inhomogeneous  pressure  wave  equation  which 
characterizes  the  MHD  and  thermoacoustic  sound  mechanisms  was  derived. 
The  inhomogeneous  equation  was  solved  assuming  a  plane  wave  pressure 
field  within  the  waveguide.  The  plane  wave  field  was  matched  at  the 
waveguide  apertures  to  a  spherical  wave  pressure  field  in  the  free  medium 
through  impedance  conditions  at  the  waveguide  apertures  om  the  wave 
field  solutions  the  MHD  and  thermoacoustic  mechanisn  -ansmitting 
sensitivities  were  derived  to  predict  the  farfield  radiation  .rom  an 
"hypothetical"  spherically  baffled  MHD  transducer. 

There  were  three  phases  of  the  experimental  investigation: 
measurement  of  the  electrical  input  impedance,  measurement  of  the  on-axis 
transmitting  sensitivities,  and  measurement  of  the  acoustic  radiation 
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directivity.  Due  to  the  weight  limitations  the  transducer  which  was  built  did 
not  resemble  a  sphere  with  a  channel  through  the  center,  as  was  assumed 
in  the  theoretical  development,  but  the  measurements  made  correlated  very 
well  with  the  theoretical  predictions. 

A  model  of  the  electrode  polarization  impedance  was  used, 
which  produced  accurate  theoretical  predictions.  The  polarization 
impedance  took  into  account  the  chemical  diffusional  impedance  and 
included  a  relatively  simple  model  for  the  double  layer  impedance.  It  can  be 
concluded  from  the  impedance  measurements  that  a  model  of  the 
transducer  which  places  the  polarization  impedance  in  series  with  the 
electrolyte  impedance  is  appropriate  for  NaCI-water  electrolytes  and  noble 
metal  planar  electrodes. 

Measurements  of  the  MHD  transmitting  sensitivity  were,  in 
general,  in  good  qualitative  and  quantitative  agreement  with  the  theory.  The 
data  comparison  shows  that  up  to  4000  Hz  the  agreement  is  within  2.5  dB. 
Above  this  frequency  the  predicted  and  measured  MHD  directivities  diverge 
significantly,  due  to  the  difference  between  the  assumed  and  actual 
transducer  geometries. 

The  measured  MHD  directivity  pattern  at  1000  Hz  showed 
excellent  agreement  with  the  computed  directivity.  At  4000  Hz  the  directivity, 
while  corrupted  due  to  the  non-spherical  nature  of  the  transducer  does 
agree  fairly  well  with  predictions  assuming  a  spherical  baffle  of  5  cm  radius. 
In  general,  at  low  frequencies,  the  transducer  makes  an  excellent  dipole 
radiation  transducer. 
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The  measurements  of  the  thermoacoustic  radiation  were  very 
difficult  to  make,  but  given  the  limited  amount  of  data  one  can  see  a 
reasonably  good  theoretical  agreement.  Since  time  was  not  available  to 
measure  reliable  directivity  patterns  for  the  thermoacoustic  source 
mechanism,  absolute  validity  of  the  sensitivity  data  is  lacking.  The  directivity 
of  the  thermoacoustic  source  is  predicted  by  the  computer  code  to  be 
omnidirectional  within  0.05  dB  for  1000  Hz.  Measurements  on  axis  and  90° 
off  axis  were  made  and  the  values  were  within  2  dB;  however,  this  is  not 
proof  that  there  was  actually  a  monopole  directivity. 

Considering  the  geometrical  and  wave  number  constraints  of  the 
theory,  all  of  the  measurements  are  in  good  agreement  and  use  of  the  theory 
for  even  a  "moderately"  exceeded  operating  frequency  range  seems 
justifiable. 

The  experiments  reported  here  were  limited  to  a  frequency  band 
500  to  13000  Hz.  At  the  low  end  of  this  regime,  the  size  of  the  testing  tank 
was  the  limiting  factor.  Frequencies  lower  that  500  Hz  would  have  been 
interesting  to  examine  since  the  lower  frequencies  necessarily  imply  larger 
particle  velocity  amplitudes.  Studies  at  higher  velocity  amplitudes  might 
have  given  some  indication  of  the  range  of  validity  of  both  the  anerture 
acoustic  impeda-  model  and  electrical  impedance  model.  However, 
facilities  for  a  larger  free  field  environment  and  a  means  to  increase  the 
source  level  were  not  available. 

With  regard  to  the  input  power  available  for  the  experiment, 
sufficient  current  was  not  available  to  test  the  limit  of  the  approximation 
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imposed  on  the  use  of  the  linearized  Nernst  equation  for  derivation  of  the 
diffusional  impedance.  Thus,  nonlinear  behavior  of  the  polarization 
impedance  at  large  current  densities  would  be  an  interesting  phenomenon 
to  study  more  closely  since  any  source  using  the  MHD  mechanism  to 
produce  "practical"  pressure  levels  (using  an  electrolyte)  would  be  operating 
in  this  nonlinear  impedance  regime. 

Since  the  constraint  on  the  small  temperature  rise  assumption  of 
the  thermoacoustic  assumption  allows  (effectively)  for  a  very  large  current 
density,  testing  the  effects  of  exceeding  this  limit  might  require  the  use  of  a 
pulsed  power  genen'.tor  or  capacitor  bank. 
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APPENDIX  A.1 


DERIVATION  OF  THE  MAGNETIC  DISPERSION  RELATION 

The  derivation  of  the  magnetic  dispersion  relation  to  follow  is 
included  to  supplement  the  discussion  of  Chapter  2.  The  small  signal 
momentum,  continuity,  and  state  equations  of  Chapter  2  are  restated  below 
with  the  exclusion  of  the  thermoacoustic  term  present  in  the  continuity 
equation,  Eq.  (2.24),  because  the  term  simply  generates  terms  of  second 
order  compared  with  the  dispersion  term  of  the  MHD  source. 

The  equations  derived  in  Chapter  2  are: 
the  small  signal  momentum  equation, 

p„u,  +  Vp=  JXB,  (2.25) 

small  signal  continuity  equation, 

Pj  +  PoV-u  =  0,  (A.1.1) 

and  small  signal  state  equation, 

p  =  c„2  6p.  (2.23) 

The  expression  for  the  current  density  J  from  Chapter  2  is  (neglecting  the 
displacement  current), 

J  =  a  (E  +  uxB)  .  (2.1) 

Inserting  the  expression  for  J  into  Eq.  (2.31), 

PqU,  +  Vp  =  a  (E  +  uxB)  X  B  (A.1 .2) 

a  momentum  equation  in  terms  of  the  new  source  field  E  and  the  damping 
term  ouxBXB  is  found.  The  divergence  of  Eq.  (A.1 .2)  is  now  taken, 

V2p  +  p^V'U,  =  aV-(E  +  uxP)XB  ,  (A.1 .3) 
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and  the  continuity  equation,  Eq.  (A.1.1),  is  differentiated  with  respect  to  time 
with  the  state  equation  substituted  for  p, 


Substituting  Eq.  (A.1.4)  into  Eq.  (A.1.3)  results  in  the  inhomogeneous  wave 
equation, 

V2p  - 1  =  aV-  (E  +  uxB)XB  .  (A.  1 .5) 

Since  the  goal  here  is  the  derivation  of  the  plane  wave  dispersion 
relation,  the  inhomogeneous  term,  aV-  E,  can  be  dropped  and  the 
homogeneous  wave  equation  can  be  written  in  one-dimensional  form, 

Pxx  '  =  0  .  {A.1.6) 

From  the  continuity  equation,  Eq.  {A.1.1)  and  Eq.  (2.29),  u,^  is  related  to  p  as 

Ux  =  -Pt/poCo^-  (A.1.7) 

Substituting  Eq.  (A.1.7)  into  Eq.  (A.1.6)  yields  the  final  form  of  the 
homogeneous  wave  equation, 

Pxx  -  I  Ptt+  Pt  =  0.  (A.  1.8) 

The  dispersion  relation  is  found  by  substitution  of  the  assumed 
harmonic  progressive  form  of  the  plane  wave  solution  into  Eq.  (A.1.8), 

p(x,t)  =  pQ  e  ,  (A.  1.9) 

which  produces 

(-k2  +  {1/Cq2)  0)^  -  j(aBg2/p^Co2)  co)  pQ  e  iCo'-®*)  =  0  .  (A.1 .1 0) 
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Hence,  the  magnetic  field  dispersion  relationship  is 


k  = 


VV^o  J 


j  (ooB, 
PoC<“ 


2^ 


1/2 


0^0  / 


(A.1.11) 


APPENDIX  A.2 


THERMOVISCOUS  PROCESSES  WITHIN  THE  WAVEGUIDE 

This  appendix  will  threat  the  justification  for  neglecting  the 
thermoviscous  processes  in  the  theoretical  analysis  of  Chapters  2,  3,  and  4. 
Primarily  the  effects  of  viscosity  are  being  neglected  so  as  to  simplify  the 
analysis  and  allov/  an  unobcured  treatment  of  the  magnetohydrodynamic 
and  thermoacoustic  acoustic  source  mechanism.  The  assumptions  made  in 
the  modeling  process  of  the  MHD  transducer  are  justifiable  in  the  framework 
of  the  theory  applied  to  the  experimental  transducer  study  and  are  not  in 
general  always  applicable. 

Losses  in  a  thermoviscous  medium  occur  due  to  two  primary 
mechanisms,  heat  conduction  and  viscous  molecular  interaction.  The  losses 
from  these  two  mechanisms  will  be  considered  in  terms  of  the  acoustic 
resonance  quality  factor,  of  the  nth  resonance  mode.  The  quality  factor 
will  be  used  to  compute  the  half  power  bandwidth  of  the  nth  harmonic  and 
compared  to  estimates  from  the  theory  cf  Chapter  3,  which  take  into  account 
only  radiation  losses. 

For  plane  wave  propagation  in  seawater  the  attenuation  due  to 
thermoviscous  absorption  at  10000  Hz  is  a-p^  of  8.1  X  10'^  Nep/m.^^  From 
Kinsler  and  Frey^®  the  expression  for  the  Q  is 

Qn  =—  (A.2.i) 

2a 

The  first  resonance  mode  is  approximately  the  full  wavelength 
pressure  mode;  for  a  waveguide  0.10  m  long  the  frequency  is 
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approximately  15000  Hz.  Taking  the  wave  speed  to  be  1500  m/s,  Q-]  is 
260000.  The  half  pov^er  bandwidth  is  therefore  0.06  Hz,  a  very  "sharp" 
resonance. 


The  attenuation  coefficient  due  to  the  boundary  layer  losses  for 
an  isothermal  waveguide  wall  is  given  as^® 


®eL~ 


1 


f 


cop. 


2poC, 


1  + 


V 


p-1 

Pr 


(A.2.2) 


with  the  constraint  on  Eq.  (A.1.1)  that 


kSvisc 


1 

«  kRg 

«  j-~ — 

500  Hz  ^  °visc 

10000  Hz 

500  Hz 


(A.2.3) 


where  5^,3^.  =  2p/copQ  is  the  viscous  boundary  layer  thickness.  For  seawater 
the  following  fluid  property  values  are  used  for  computing 


p  =  0.001  N  s/m2 

Y=1.01 

p^=  1026  kg/m3 
Cq=  1 500  m/s 


coefficient  of  viscosity 
ratio  of  specific  heats 
mass  density 
adiabatic  sound  speed 


For  an  aperture  radius  Rg  of  0.02  m  the  theory  is  valid  for  the  entire 
frequency  spectrum  from  500  to  10000  Hz,  with  bounds  of  0.000234  « 
0.0419  «  955.0.  The  OqlIS  0.0035  Nep/m  which  yields  a  of  5984  and  a 
half  power  bandwidth  of  2.5  Hz. 

Figure  A.2.1  shows  a  graph  of  the  maximum  acoustic  pressure 
predicted  by  the  theory  presented  in  Chapter  3,  which  takes  account  of 
radiation  losses  only.  It  is  seen  from  this  graph  that  the  half  power  bandwidth 
of  9400  Hz  is  very  large  compared  to  either  of  the  numbers  calculated 
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above.  Also  plotted  in  the  graph  is  the  maximum  pressure  for  a  sound  speed 
of  650  m/s;  again  a  large  bandwidth  is  noted.  The  conclusion  to  be  drawn  is 
that  thermoviscous  mechanisms  play  an  insignificant  role  ii'  the  dynamics  of 
the  acoustic  processes  occurring  within  the  waveguide  of  the  specified 
geometry  and  frequencies  of  interest  noted  here. 
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FIGURE  A.2.1 

PEAK  MODAL  PRESSURE  IN  THE  MHD  WAVEGUIDE  AS 
A  FUNCTION  OF  FREQUENCY  Q  DETERMINATION  PLOT 
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APPENDIX  B 


PRESSURE  FIELD  SPHERICAL  HARMOMIC  FUNCTION  EXPANSION, 
RADIATION  IMPEDANCES,  AND  FARFIELD  PRESSURE 


The  calculation  of  the  mutual  radiation  impedance  requires  that 
the  acoustic  field  be  specified.  Therefore  the  first  task  is  to  specify  the  source 
and  boundary  conditions.  In  this  case  the  sources  are  two  circular  pistons 
on  the  surface  of  a  rigid  spherical  baffle  radiating  into  an  acoustic  free  space. 
The  pistons  are  positioned  180°  apart  on  the  sphere  which  gives  the 
radiation  an  axis  of  symmetry  through  the  piston  centers.  Thus  there  is  only 
one  angle  y,  on  which  the  field  will  be  dependent  (see  Fig.  B.1).  The  results 
that  are  presented  follow  directly  from  Sherman's^^  results  with  the  above 
conditions. 

The  pressure  at  any  point  in  free  space  outside  the  sphere  can  be 
expressed  as  the  sum  of  the  pressure  fields  generated  by  the  individual 
pistons, 

=  Pitafr.f)  +  P-lJW)  .  (B.1) 

The  spherical  harmonic  function  expansion  is  stated  below  for  the 
pressure  field  resulting  from  one  piston  on  the  sphere. 


P.t(m(''’'1')=±iPfmP|mP<i^) 


oo 

I 

n=o 


(Pn-I  (COS0)  -  (cos©))  (cos\i/)  h^(kr) 


x-kR 


(B.2) 
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FIGURE  B.1 

APERTURE  GEOMETRY 
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In  Eq.  (6.2).  and  Pp  are  the  ordinary  Legendre  polynomials  of  order  n±1 
and  n,  cos0,  and  cosy  are  the  arguments  of  the  Legendre  polynomials 
where  0  is  the  half  angle  of  the  source  radius,  and  \j/  is  the  coordinate  angle 
referenced  to  the  source  symmetry  axis.  h^„(kr)  is  the  Hankel  function  of  the 
first  kind  of  order  n  and  is  a  function  of  kr,  where  r  is  the  radial  distance  to  a 
point  in  free  space  and  R  is  the  sphere  radius. 

The  calculation  of  the  mutual  radiation  impedance  follows  directly 
from  Eq.  (B.2).  The  pressure  is  evaluated  on  the  surface  of  the  sphere,  r=  R, 
substituted  into  Eqs.  (3.4a,b)  and  integrated  over  the  angular  sector  of  the 
source  aperture  from  \jr  =  0  to  V  =  ©•  t^or  the  case  of  the  waveguide  geometry 
used  in  the  experiment  the  mutual  impedance  is  as  follows 


2  = 


CX3 

I 

n=:0 


2.  1 


^  (P„.i  (COS0)  -  P^^^  (COS0))  h^(kR) 


2n+1 


-f  hVx) 

dx  ' 


x=kR 


(B.3) 


where  p  is  the  relative  temporal  phase  between  the  sources. 

Equation  (B.3)  is  valid  for  both  apertures  provided  they  are  of 

equal  radius  or  as  stated  above,  angular  sector  0.  The  expression  is 
independent  of  the  magnitude  of  the  particle  velocity  since  the  aperture 
velocity  distribution  was  assumed  to  be  purely  planar  motion. 

Calculation  of  the  farfield  pressure  is  accomplished  using  the 
farfield  approximation  to  Eq.  (B.1),  which  is  to  say  the  approximations  are 
made  in  Eq.  (B.2).  Note  that  Eq.  (B.2)  is  an  exact  solution  to  the  pressure 
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field  for  the  stated  boundary  conditions.  However,  when  in  the  farfield,  the 
large  argument  form  of  the  Hankel  functions  save  much  computational 
overhead.  Without  proof  the  farfield  analog  of  Eq.  (6.2)  is  stated  below.  See, 
again,  Sherman  for  a  more  complete  discussion. 


Pff  = 


JPlmCfn,® 


■UW)  X 


{COS0)  -  (COS0))  P„(-cos\k) 


d  .  1 ,  . 


+  uW  X 


(P„.^  (COS0)  -  P^^^  (COS0))  P„(cosv) 


(B.4) 
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APPENDIX  C 


REFLECTION  COEFFICIENT  EQUATIONS 

Due  to  the  complexity  of  the  coefficients  a  hierarchy  of  dummy 


variables  are  defined  below  to  construct  more  simplistic  forms. 

A  =  e-j'^^  (C.1a) 

B  =  ej'^^  (C.1b) 

C=  ei'^^-  ei'^ 

D=  8‘i'^^  (C.1d) 

The  dummy  variables  for  the  MHD  reflection  coefficient  equations  are  as 
follows. 

A  =  BD(1  -  Z3/p„c,)  -  BC  Z>,c,  {C.2a) 

i  =  AD(1+Z3/p,c,)-BDZ>,c,  (C.2b) 

C=BC(1-Z3/p,c,)+AC7^/p^c,  (C.2C) 

Q  =  AC(1  +  Z3/P0C0)  +  AD  Z^/p^Co  (C.2d) 

E  =  BC(1-Z3/p,c,)-BDZ>,c,  (C.2e) 

F  =  BD(1-Z3/p,c,)  +  ADZ>„c,  {C.2f) 

Q  =  AC(1  +  Z3/p,cJ  -  BC  Z Jp„c,  {C.2g) 

H  =  AD(1  +  Z3/p„c,)  +  AC  Z^/p„c,  (C.2h) 

The  coefficients  for  the  MHD  reflection  coefficient  equations,  Eqs.  (3.1 8a, b), 
are  as  follows  in  term  of  the  above  underlined  dummy  variables. 

Ci  =  AF-S£  (C.3a) 

C2  =  AH-  B£  +  C£-0E  (C.3b) 
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Cg  =  QH-QQ. 

(C.3c) 

C4  =  ££-AS 

(C.3d) 

C5=  D.£+  £F-fi£-Aid 

{n.3e) 

Cg  =  DF-SM 

(C.3f) 

The  dummy  variables  for  the  thermoacoustic  reflection  coefficient  equations 

are  as  follows. 

A  =  BD(1  -  Z3/p,c,)  +  BC  ZJp,c, 

(C.4a) 

B  =  AD(1+Z3/p,cJ-BDZ>,c, 

(C.4b) 

£  =  BC(-1  +  Z3/P0CQ)  -  AC  Z^/PqCq 

(C.4c) 

D  =  -AC{1  +  Z3/PQCJ  +  AD  Z  JPqCq 

{C.4d) 

E=BC(1.Z3/p„c,)  +  BDZ>„c„ 

(C.4e) 

E  =  BD(-1  +  Z3/P0CJ  -  AD 

(C.4f) 

e  =  AC(1  +Z3/poCo)-BCZ>gCo 

{C.4g) 

H  =  -AD(1  +  Z3/p,c,)  +  AC  Z  Jp„c, 

(C.4h) 

The  coefficients  for  the  thermoacoustic  reflection 

coefficient  equations, 

Eqs.  (3.30a,b),  are  as  follows  in  term  of  the  above  underlined  dummy 

variables. 

Cl  =  AF-BE 

(C.5a) 

C2=  AH-  Ba  +  CF-DE 

(C.Sb) 

C3  =  BH-Da 

(C.5c) 

C4  =  a£-AB 

(C.5d) 

C5=  a£+aF-Ea-AH 

(C.5e) 

Ce  =  DE-EM 

(C.5f) 

APPENDIX  D 


COMPUTER  PROGRAMS 


The  numerical  computation  of  the  theoretical  relations  presented 
in  the  text  of  this  thesis  were  computed  by  the  FORTRAN  programs  MHD  and 
THERMO.  MHD  is  the  analysis  code  for  the  magnetohydrodynamic  source 
mechanism  and  is  capable  of  computing  the  following  relations  as  a  function 
of  frequency. 

1)  Aperture  Reflection  Coefficient 

2)  Complex  Aperture  Pressure  and  Particle  Velocity 

3)  Complex  Aperture  Real  and  Self-Impedance 

4)  Complex  Electrical  Terminal  Impedance 

5)  MHD  Transduction  Process  Power  Efficiency 

6)  Transmitting  Current  and  Voltage  Sensitivity 

7)  Farfield  Pressure  Directivity  (single  frequency) 

The  input  parameters  for  "MHD"  are 

1)  Frequency  Band  and  Frequency  Step  Size 

2)  Waveguide  Dimensions 

3)  Electrode  Dimensions  and  Position  within  the  Waveguide 

4)  Spherical  Baffle  Radius 

5)  Fluid  Properties  Inside  and  Outside  the  Waveguide 

6)  Conductivity  of  the  Fluid  Inside  the  Waveguide 

7)  Magnetic  Induction  and  Current  Density  Amplitude  in 
the  Waveguide 

THERMO  is  the  analysis  code  for  the  thermoacoustic  source 
mechanism.  The  following  relations,  as  a  function  of  frequency,  are 
computed  by  THERMO. 


1)  Aperture  Reflection  Coefficient 

2)  Complex  Aperture  Pressure  and  Particle  Velocity 

3)  Complex  Aperture  Real  and  Self-Impedance 

4)  Quadratic  Current  Transfer  Function  (QCTF) 

5)  Thermoacoustic  Transduction  Process  Power  Efficiency 

6)  Quadratic  Voltage  Transmitting  Sensitivity  (QVTS) 

7)  Farfield  Pressure  Directivity  (single  frequency) 

THERMO  has  the  same  input  parameter  list  that  is  given  above 
for  MHD  except  the  following  properties  for  the  fluid  within  the  waveguide 
are  needed;  specfic  heat  capacity  at  constant  pressure  and  coefficient  of 
thermal  expansion. 
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PROCRAM  MHDUNPUT.  OUTPUT.  DATA.  TAPE2- INPUT.  TAPE4-DATA 
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PROCRAn  MHD  PERFORMS  A  FPEOUENCV  DOMAIN  ANALYSIS  OF  THE  ELECTROACOUSTIC 
TRANSMITTINC  CHARACTERISTICS  OF  A  PLANE  WAVE  MODE  MHO  TRANSDUCER 
MHO  GENERATES  GRAPHIC  OUTPUT  FOR  THE  FOLLOWING  TRANSDUCER 

characteristics 


MAGNITUDE  OF  THE  REFLECTION  COEFF  (MRNL)  VERSUS  FREQUENCY  <FREQ> 
MAGNITUDE  OF  APERATURE  PRESSURE  (RPNL>  VERSUS  FREQ 
08  LEVEL  OF  THE  APERTURE  PRESSURE  IRPNL)  VERSUS  LOGIO(FREO> 

OF  APERTURE  PART  VEL  (RVNL)  VERSUS  LOOIO(FREQ) 

IMAGINARY  <IADM>  VERSUS  REAL  (RADM)  INPUT  ELECTRICAL  ADMITTANCE 
(ZELI)  •  “  (2ELR>  -  “  IMPEDANCE 

06  LEVEL  OF  RADIATED  ACOUSTIC  POWER  (PR)  VERSUS  LOGIO(FREO) 

REAL  COMPONENT  OF  INPU'  ELECTRICAL  POWER  (PER)  VERSUS  FREQ 
-  -  -  -  (PERECT) 

10)  DB  LEVEL  or  TRANSDUCTION  POWER  EFFICIENCY  (EFFIC)  VERSUS  LOGIO(FREO) 
ANSMITTING  VOLTAGE  SENSITIVITY  (TVSL)  VERSUS  LOGFREO 
-  CURRENT  "  (TCSL) 

13)  APERTURE  ACOUSTIC  SELF  IMPEDANCE  REAL  COMP  (2FMR)  VERSUS  Ff.-.Q 
"  "  IMAGINARY  -  (ZFMI) 

MUTUAL  "  Rf-AL  *  (2MR)  - 

••  -  IMAGINARY  •  (ZMI  ) 

lAONlTUDES  OF  ZM  /  ZFM  VERSUS  FREO 


1) 

MAGNITUDE 

2) 

MAGNITUDE 

3) 

OB  LEVEL  ( 

4) 

“  “  ( 

9) 

IMAGINARY 

6) 

" 

7) 

DB  LEVEL  ( 

0) 

REAL  COMP( 

9) 

REACTIVE 

10) 

DB  LEVEL  ( 

11) 

DD  LEVEL  ( 

12) 

OB 

13) 

APERTURE  i 

14) 

" 

19) 

■ 

16) 

• 

17) 

RATIO  OF  ' 

18) 

PRESSURE  ( 

19) 

ANGLE  OF  t 

DATA  input  is  THROUGH  THE  •.  DATA  STATEMENT  AT  THE  END  OF  THE 

MAIN  PROGRAM 


COMMON  /ZIP/  A,D.L.LY,LZ.SRAD.R0D.  COD.  ROM.  C0M.P1.6H 
I.  CO.  CON.  FO.  FF.FO 

0 1 MENS ION  FREQ ( 1 000 ) . FLOG ( 1 000 ) . PR ( 1 000 ) .  PSPt ( 1 000 ) 

1. RPNL( 1000).  RVNL(1000>.RADM( 1000). I ADM( 1000 ) . MRNL ( lOOO) 

2.  TC^c ( 1 000 ) .  TVSL ( 1000 ) . ANGHAX ( 1 000 ) .  R A ( 1 000 ) 

3.  ZMR(IOOO).  ZMI (1000). ZFMR< 1000). ZFMZ( iOOO).  RATIO! 1000) 

4.  ZCLRdOOO).  ZELI  ( 1000).  RVL(  1000) 

9. P£R( 1000).  PERECT! 1000).  EFFIC( 1000) 
complex  RNL(IOOO).RL<lOOO).PNL(IOOO).PL(1000) 

1.  VNL<1000).VL(1000> 

2.  CLNLA. GLNL8.  GRLA. ORLD. GLBD. GLAB. GRAA.  GRBA 

3.  GXLNLA.  GXLNL8.  CXRLA.  GXRLB.  PE 

4. ADMl(100O).ZEL<I00O> 

PEAL  IAOM.K.L.LY.LZ.MRNL 
INTEGER  I.N 

N  •  (FF-FO)/FO  ♦  1 
C 

C  PRIMARY  CALCULATION  LOOP 
C 


DO  1  I-l.N 

C  CALCULATE  THE  WAVE  NUMBER  AND  FREQUENCY 

K  •  <2  •PI*FO  ♦  FL0AT(I-l)*2  •PI*FO)/COD 
FREOd)  •  MC0D/<2  •PI) 

FLOG(I)  •  AL0G10(FRE0(I)) 

C  REFL  EVALUTES  THE  APERTURE  REFLECTION  COEFFICirNT 
CALL  REFL(K.  RNL(I).RL<I).MRNL(I)>  ZMR  ( 1 ) .  ZMI  ( 1 ) 

1.  ZFMR(I).  ZFMI(I>> 

C  GREEN  EVALUATES  THE  PLANE  WAVE  GREENS  FUNCTION 

CALL  OREENIK.  A.  8.  L.  RNL<  I ).  RL<  1 ).  GLNLA.  OLNLB.  GRLA.  GRLB 
1. GXLNLA. GXLNL8.  OXRLA. GXRLB. GL8B.  GLAD.  GRAA.  0R8A) 

C  GO  TO  90 

C  ZM  AND  ZS  ARE  THE  MAGNITUDES  OF  THE  APERTURE  SELF  AND  MUTUAL  IMPEDANCES 
ZM  •  (ZHR(I)*«2  ♦ZMI<I)««2  )••  9 
ZS  -  (ZFMR(I)**2  ♦ZFM1(I)**2  >•*  9 
RATIO! 1)  •  ABS(ZM/ZS> 

C  PRES  EVALUATES  THE  APERTURE  PRESSURE  IN  THE  WAVE  GUIDE 
CALL  PRES (GLNLA.  CLNLB.  GRLA.  GRLB. PNL ( 1 ) .  PL < 1 ) ) 

RPNLd)  •  CAB8(PNLd)) 

PSPLd)  *  20  •ALOOIOCRPNLd)/  000001) 

C  VEL  EVALUATES  THE  APERTURE  PARTICLE  VELOCITY 
90  CALL  VcL (K.  GXLNLA. OXLNLD.  OXRLA.  OXRLO.  VNL (I ) .  VL (Z ) .  UJ .  U2. P 1 .  P2) 

C  GO  TO  100 

RVNLd)  •  CA88<VLd>> 

C  ADM  EVALUTES  THE  INPUT  ELECTRICAL  ADMITTANCE 
CALL  ADM(H.  GLBB.  GLAB.  GRBA.  GRAA.  ADMl  (I ) ) 

RADMd)  •  REALCCMPLXd  .0  )/ADN|d)) 
lADMd)  •  AlMAOtCMPLXd  .0  )/ADMld)) 

C  ZCL  18  THE  COMPLEX  ELECTRICAL  IMPEDANCE 
ZEL(')  •  CMPLXd  0  Ol/ADMld) 

ZELRd)-  REAL(ZCLd)> 

ZELld)-  AIMAOiZELd)) 

C  RADPR  EVALUATES  THE  RADIATED  ACOUSTIC  POWER 
CALL  RADPRILY.LZ.  ZFMRd  ) .  VL  d  ) .  PRR  J 
PRd)  -  20  •ALOGlOlPRR) 

C  PE  IS  THE  AVERAGE  COMPLEX  INPUT  ELECTRICAL  POWER 
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PE  -  CMPLX(CO*CO*CB-A)*Ly»L2/<2  •CON).0  0)  ♦ 

1  CrtPLX(0  .  CD*CO«8n*Drt*LY*LZ/<K*ROD*COO*S  M* 

2  <CL8e*CLAS'CRBA«CRAAI 
PER(I)  •  P£AL<Pe> 

PERECTd)  •  AinAC(P£> 

C  EFFIC  IS  THE  TRAMSDUCTION  PROCESS  EFFICIENCy 
EFF|C<n  -  20  •AL0C10<PRR/PER<in 
C  CO  TO  1 
C 

C  SEN  EVALUATES  THE  PRESSURE  RADIATION  PATTERN.  TRANSMITTING  CURRENT  AND 
C  VOLTAGE  SENSITIVITIES  ANO  THE  ANCLE  OF  THE  MAIN  LODE 
C 

C  IF  A  DIRECTlViry  PLOT  IS  DESIRED  SET  FLAG  •  1  0  < DO  SCALE) 

C  ANO  /  OR  SET  FLAG2  «  1  0  (LINEAR  SCALE) 

C  SET  FLAG  AND  FLAC2  *00  OTHERUISE 
C  SET  THE  FLAGS  IN  THE  IF  STATEMENTS  DELOW.  ONLV 
IOC  FLAG  "00 

FLAG2  *00 
IF<I  EO  1>  FLAG-1  O 

IF<I  EO  1)  FLAC2*  O 

C  IF  THE  SCNSITIVITy  ON  THE  X  AXIS  IS  DESIRED  SET  FLAGS  -  0  0 
C  IF  THE  MAXIMUM  SENSITIVITY  IS  DESIRED  SET  FLAGS  -  1  0 
C  FLAGS  MUST  EQUAL  1  0  WHEN  DIRECTIVITY  PLOTS  ARE  REQUESTED 
C 

FLAGS  •  t  0 
C 

CALL  SEN(FReO( I ) . U1 . U2. PI . P2. CXLNLA. GXLNLD.  CXRLA.  GXRLB 
1.GLB8.CLAQ.  CRDA.  CRAA.  7CSL(I>.  TVSLI I  ).  ANCMAXd  ) .  FLAG.  FLAG2.  FLAGS) 

C  URITE(4.*)  FREOd).FLOG(I>.TCSLd) 

1  CONTINUE 

C 

C  02PLT  PRODUCES  A  NEUTRAL  PLOT  FILE  FOR  PLOTTING  ON  THE  TEKTRONICS 
C 

CALL  02PLT<N. FREO  MRNL.  1 ) 

CALL  02PLT(N.  FREQ.RPNL.2^ 

CALL  02PLT(N.FLOG.PSPL.S) 

CALL  02PLT(N.  FREQ.  RVNL.  4) 

CALL  02PLT(N.  FREO.  IADM.  S> 

CALL  D2PLT(N.FRE0.RADM.  6) 

CALL  02PLT<N.FLOG.PR.7> 
call  02PLT<N. FREQ,  per,  B) 

CALL  02PLT(N.FREO.PERECr,9) 

CALL  02PLT(N. FLOG.  EFFIC.  10) 

CALL  02PLr(N.FL0G.  TVSL.  ID 
CALL  02PLT<N.FLOG.  TCSL.  12) 

CALL  02PLT(N.FR£Q.  ANGMAX.  13) 

CALL  02PLT(N.FRE0.  ZMR. 14) 

CALL  02PLT(N,FR£0.  2MI.  IS) 

CALL  D2PLT(N.FREQ.  ZFMR.  16) 

CALL  02PLT(N.FRE0.  ZFHI,  17) 

CALL  D2PLT(N.FREQ.RAT10.  18) 

STOP 

END 

C 

c 

BLOCK  DATA 
C 

c  var;-^ble  inputs  to  program  mho 
C 

C  A 

C  0 
C  L 
C  LY 
C  LZ 
C  SRAO 
C  ROD 
C  COO 
C  ROM 
C  COM 
C  6M 
C  CO 
C  CON 
C  PO 
C  FF 
C  PD 
C 

COMMON  /IIP /A.  0.L.  LV.  LZ.  SRAD.  ROD.  COD.  ROM.  COM.  PI.  BM 
1.  CD.  CON.  FO.  FF.FD 
real  L  LY.tZ 

DATA  A.OL.LY.LZ  SRAO.  ROD.  COO.  ROM.  COM.  P I .  DM,  CD.  CON,  FO.  FF.  FD/ 

1-0  05,0  05.  10.0  0310.0  0301.  05.1041  .1550.990.1401 
2.3  14159.0  40,  MA36  5  20.250  .10000  .50  / 

C 

C  PHYSICAL  PROPERTIES 
C 


C  6X  NACL-H20  SOLUTION 
(, 

ROD  • 

1041 

COD  •  1536 

1600  (MEASURED) 

C  MERCURY 

ROD  • 

13600 

COD  -  1450 

C 


LEFT  BOUNDARY  OF  SOURCE  VOLUMF  (METERS) 

RIGHT  «  V  - 

total  WAVEGUIDE  LENGTH  (METERS) 

TRANSVERSE  WAVEGUIDE  LENGTH  Y-DIR  (METERS) 

Z-DIR 

RADIUS  OF  SPHERICAL  BAFFLE  (METERS) 

FLUID  DENSITY  IN  WAVE  GUIDE  (K0/HETER**3) 

•  SOUND  SPEED  "  “  (METERS/SEC) 

fluid  density  in  free  medium  (KC/MeTER*«3) 

“  '3VN0  SPEED  '  "  (METERS/SEC) 

MAGNETIC  INDUCTION  FIELD  AMPLITUDE  (WEDERS/METER**2) 

CURRENT  DENSITY  AMPLITUDE  (AMPERE5/METER**2> 

ELECTRICAL  CONDUCTIVITY  OF  FLUID  IN  WAVEGUIDE  d/OHMS  METER) 
initial  calculation  FREOUCNCY  (HERTZ) 

FINAL  "  -  - 

FREQUENCY  STEP  SIZE  (HERTZ) 
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rj*  ‘J*  'J  •j*  V  ✓  V  /  •*-  A 

‘fh  ■  *  "  A  Kr  *.'  •.*'  1^,  I*-  **•  S 


CND 

C 

C»«*«*«  •«•••*<  . . . . . . . 

c 

suoRour t^e  rcflck, rsnl. rsl. RSNLn, ZKR.  zni*  zsr.  zsi ) 
c 

C  SUBROUTINE  RCFL  CALCULATES  THE  COERPICENTS  AND  ROOTS  OP  THE 
C  REFLECTION  COEFFICIENT  COMPLEX  QUADRATIC  EQUATION 
C 

COMMON  /ZIF/A.  fl.  L^LV.LZ.  SRAD.  ROD.  COD.  ROM.  COM.  PI.  0M,  CD.  CON.  FO.  FF.  FD 
complex  ha.  hD.  KL?.  ZS.  zn.  A1 . 0 1 .  C  I .  DI  .  a?.  B2.  C7.  D2.  E2.  F?.  C3.  H2 
1 .  C 1 1.  C?2.  C33.  C44.  CSS.  C66.  RSNL.  RSL.  RLL.  RLNL 
REAL  L  LY.LZ.K.PC 
PC  •  ROO*COD 
HA  •  CMPLXIO  O.K*A) 

KB  •  CKPLXIO  O.K«B) 

HL3  -  CMPLXIO  O. K*L/r) 

C  CALCULATE  the  aperture  IMPEDANCES 

CALL  ntSPIK*COD.  ZTtR.  ZMI.  ZSR.  ZSI  > 

ZS  •  CMPLXIZSR.  ZSn/CMPLX<PC.O  0) 

ZM  •  CMPLXIZMR. Zni>/CMPLX<PC.  0  0> 

C  CALCULATE  THE  QUADRATIC  EQUATION  COEFFICIENTS 
Cl  •  CEXPCKA)  •  CEXP<k0) 

01  «  CexP<*HA>  -  CCXP<-HB> 

Ai  -  CCXP<-HL?> 

01  -  CExPIKtS) 

A?  •  Ol*Dl«<CnPLX<l  .0  0)-ZS>-8I*CI»ZM 

02  •  AI*01*<CI1PLX<I  .0  0)»ZS>-BI*DI*ZM 

cr  •  OMCi*<cnPLx<i  .o  o>-zs>*ai*ci*zm 

02  -  AI*C1*(CMPLX(1  .0  0>*ZS)*Ai*Dl«ZM 

E2  •  Dt«Cl«CCnPLX(l  .0  0>-ZS)*01*Dl*ZM 

F2  ■  OI*01*<CMPLX< 1  .0  0)-ZS)*Al»Dl»ZM 

C2  •  Al*Cl*<CnPLX< 1  .0  0)*ZS)*DI*CI*ZM 

H2  -  A1*01*(CMPLX<1  ,0  0)»ZS>*A1*CMZM 

Cl  I  •  A2«F2  -  02«£2 

C22  •  A2*H2  -  B2*C2  ♦  C2*F2  -  02«E2 

C33  •  C2«H2  -  02*C2 

C44  ■  C2*E2  •  A2*C2 

CSS  •  02*£r  ♦  C2«F2  -  B2*C2  -  Ar*H2 

C66  •  02*F2  -  B2*H2 

C  CALCULATE  the  COMPLEX  ROOTS  USINC  THE  IMSL  ROUTINE  ZOADC 
C  the  ROOT  WITH  THE  MACNITUDE  LESS  THAN  ONE  IS  PHYSICALLY 

c  meanincful  the  other  is  discarded 

CALL  ZQADC(C33.C22.CU.RSNL.RLNL.  lER) 

CALL  Z0A0C<C66.CSS.C44.RSL.RLL.  lERI 
RSNLM  •  CAaS(RSNL> 

RETURN 

END 

C 

. . . 

c 

subroutine  MISPIU.RM.CXM.R.CX) 

c 

c  MISP  EVALUATES  THE  SELF  AND  MUTUAL  ACOUSTIC  IMPEDANCES  IF  THE  HAVECUIOE 
C  APERTURES  THE  RESULTS  ARE  INPUT  TO  REFL  TO  EVALUATE  TmE  REFLECTION 
C  COEFFICIENTS 

C  MISP  USES  the  equations  DEVELOPED  BY  C  H  SHERMAN  TO  EVALUATE  TME  MUTUAL 
C  AND  SELF  impedances  OF  TWO  CIRCULAR  SOURCES  IDO  DECREES  A.’»ART  ON  THE 
C  SURFACE  OF  A  SPHERICAL  BAFFLE 
C 

DIMENSION  P<  1000 >.  WH( 2000).  Ji  ( 1000).  YI ;  1000). 

IPMIOOO) 

COMMON/  Z  IP/A.  B.  L.  LY.  LZ .  SRAD.  ROD.  COD.  ROM,  C  3M  P  I ,  OM.  CD.  CON.  FO.  FF.  FO 
REAl  Ul.  YI.L.LV.LZ 
INTEGER  I.n.N 

C  N  IS  THE  NUMBER  OF  TERMS  IN  THE  SVM  REQUIRE'/  FOR  CONVERGENCE  AT 
C  LARGE  VALUES  OF  H«RO 
N  ■  30 

C  RD  IS  THE  EFFECTIVE  RADIUS  OF  THE  CROSS  SECUON  OF  A  RECTANGULAR 
C  WAVEGUIDE 

RD  -  S0RT(LY«LZ/PI> 

T  -  ATAN<RO/SRAD) 

0  •  ROM*COn«SRAD*SRAD/<RD*RD) 

C  INITIALIZATION  OF  LEGENGRE  POLVNOMIALI 
P«I)  •  I  0 

p<r)  •  I  o 

p(3)  •  COSCT) 

PUD  •  I  O 
PU2)  --10 
M  •  N*2 

C  CALCULATION  OF  LEGENDRE  POLYNOMIALS  BY  RECURSION 
00  1  H«3>  ri 
RI  •  FLOAT<H) 

R2  •  RI  -  I  0 

P(H*D  •  P<3)«P(K)*<2  cRDI  '-P<H-D»RI/(RD1  > 

PUH)  •  PU2)*PUH>1>»>2  •R2*l  )/<R2  1  > -P  U  H-2 )  •P2/ <R2*  I  ) 

1  CONTINUE 

C  THE  BESSEL  FUNCTION  EVALUATIONS  ARE  PERFORMED  BY 
C  MHOOUR  and  MMOSVN  (IMSL  ROUTINES) 

X  «  W»SRAD/COM 


V  V-  A  O'  r.  *  .*-F  V  . 


CALL 

MMBSJR(X. 

D.NO.  J1«WK.  lER) 

CALL 

MMOSYNIX. 

S.N-3.  Yl.  ICR) 

ZRMl 

•  0  0 

2CMI 

-  0  0 

ZRSl 

•  0  0 

ZCSl 

•  0  0 

N2  • 

N-l 

C  TKC  SERIES  ARE  SimnEO  IN  THE  1  LOOP 
C  70  CALCULATE  THE  SELF  AND  nUTUAL  inPEDANCES 
DO  3  t-l<N2 
SI  •  FLOATTt) 

D  •  (((SI-1  )/X>*Jl(ll-Jl<|«lM*»3  ♦ 
l<((St-l  >/X>*Yl(I>-yi<I«l>)*«2 
D!  •  <(SI-1  )/X)*( J1 (t)*Jl( I )*yi ( I MVI ( 1 > >  - 
Mui(i»*ji(t*i)«yt(t >*yi(i*i}> 

D2  -  -Jt  (I>«yi(Ull«Jl(l«l)«Yl  (  1  » 

RS  ■  02/0 
CPU  •  -01/D 

C2n  •  <1  /<2  *(51-1  >»i  )>*<(p(n-p(i»2M**2>«Pi(n 

C2S  -  (1  /(2  *<51-1  »♦!  I >*(P( I )-P( 1»2) »#*2 

2Rni  -  C2n«RS«>Zftni 
icni  »  C2n*CPL*2cni 
2RSt  •  C2S*RS*ZRS1 
ICSl  •  C2S*CPL*2CS1 
3  CONTINUE 
C  RESULT 

Rn  •  C*7Rni 

CxM  ■  C*2Cttl 

R  •  C*2RM 

CX  •  C*2CS1 

RETURN 
END 
C 

c*«*«***»««**.  >«**««**«t««»««****«»***»****»«**«**«*t**»«**»«***i^ 

c 

SUBROUT INE  CNE£N(K.  A. B. L.  RNL.  RL.  CLNLA.  CLNLB.  CRLA. CRLD 
1.  CXLNLA.  CXLNLO.  CxRLA.  CXRLD  clob.  clad»  craa.  croa> 

c 

C  CREEN  evaluates  ThE  PLANE  WAVE  vREEN'S  FUNCTION  FOR  THE  WAVE  FIELD  IN 
C  THE  WAVECUIOe  AT  VARIOUS  LOCATIONS 
C 

^  COnPLEX  CC.KA.KO.HLS.RNL.RL’CLNLA,  CLNLO.  CRLA.  GRLB.CLDD.  CLAD 

1 .  CRAA.  CROA.  LA.  LB.  HA.  RB 

2.  CXLNLA.  CXLNL9.  CXRLA.  CXRlD 
REAL  K.L 

KA  •  CnPLXTO  .K*A) 

KB  •  CnPLX<0  .H«B> 

KL2  ■  CHPLXIO  .K*L/2  > 

CC  -  CnPLXCO  *1  )/(Cr(PLX<2  *K.O  )»(CHPLX(1  .0  >-RNL*RL>> 

LA  •  CC*(CEXP<KA>  ♦  RL«CEXP(-KA) > 

RA  •  CC»<CEXP(-KA)*flNL*CEXP<KA>) 

LO  •  CC»<C£XP(KS>»RL«C£XP(-HB)> 

RB  •  CC*<CEXP(-KB»*RNL»C£XP<K0> I 
C  the  CREEN 'S  functions 

CLAB«  L0*<CEXP<-KA>»RNL»CEXP(KA>) 

CRAA«  RA«{CEXP<KA)*RL*CEXP(-KA) > 

CLBB-  L0*<CCXP<-KD)»RNL*CEXP(KB) > 

CROA-  RA*(CEXP(KB)*RL*CEXP(-KD>) 

CLNLA  •  LA»(CEXP<KL2»*RNL*CEXP(-KL2M 
CLNL0  •  L8*(CEXP<KU2»*RNL*CEXP<-KL2) > 

CRLA  ■  RA»(CEXP(KL2>*RL*CEXP(-KL2) ) 

CRLO  •  R0*(CEXP(KL2»*RL»CEXP(-KU2> > 

C  DERIVATIVES  OF  THE  CREEN'S  FUNCTIONS  (WHT  X) 

CXLNLA  •  CMPLX(0  0. -K>*LA*<CEXP<KL2>  -  RNL^CEXPI -KL2) I 
CXLNLB  •  CftPLX<0  0. -Kl  •LB« (CEXP  ( KL2 >  -  RNL«CEXP  ( -KL2 )  i 
CXRLA  •  C»PLX(0  0.K)«RA*(CEXP(KL2)  -  RL*CEXP ( -KL2> ) 

CXRL8  •  CMPLX(0  O.K)cR0*(CEXPCHL2I  -  RL*CEXP I -KL2) ) 

RETURN 

END 

C 

C*********************!  . . . . . . 

c 

subroutine  PRES(CLNLA. CLNLD. CRLA, CRLB.  PNL- PL) 

C 

c  PRES  Evaluates  the  aperture  presure 
c  PNL  IS  the  complex  pressure  at  The  left  aperture 
C  PL  "  ^  -  n  -  RJQHT 

c 

COMMON  /IIP/  A.  D.  l.LY.  LI.  SRAD.  ROD.  COD.  ROM  COM.  PI.  DM.  CD 
I.CON.  FO.FF.FO 

complex  PNL.  PL.  CLNLA.  CLNLB.  CRLA.  CRLB 
PNL  •  CMPLX(0M*CD  O  )*(CLNLA  -  CLNLB) 

PL  ■  CMPLX<8M*CD.  0  )*<CRLA  *  CRLB) 

RETURN 

END 

C 

c«*«  . 

C 

subroutine  VEL(K  CXLNLA.  CXLNLB.  CXRLA.  CXRLB.  VNL.  VL.UI.  U2.  pi.  P2) 

C 
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j»  V*  ^ 


-r-V' 


A*" 

i  -  ^  - 


C  V£L  evaluates  THE  PAATtCLC  VELOCITY  AT  THE  APERTURE 
C  VNL  IS  The  complex  PARTICLE  VELOCITY  AT  THE  LEFT  APERTURE 
CVL---  - 

C 

CmnON  /ZIP/  A«8.L.LY.LI«SRAD.R0D.C0D»R0M.C0M.  PI.Bn«CO.  CON.FO.  FF*FD 
COMPLEX  VL. VNL.  OXLNLA. CXLNLO. CXRLA.  OXRLB 
REAL  K 

VNL  •  CMPLXIO  O. -6n*CO/(H«ROO*COO)}*(CXLM.A  -  6XLNL8» 

VL  •  CRPLXIO  0.>Bn*CO/<K*RuO*COD})*<CXRLA  «>  OXRLB) 

U1  -  CASS<VL) 

U?  •  CABS<VNL) 

UM  •  RCAL<VL> 

U?S  •  -REAL (VNL) 

C  THE  CONDITIONALS  BELOW  EVALUATE  THE  RELATIVE  PHASES  OF  THE  APERTURE 
C  PARTICLE  VELOCITIES  REFERENCED  TO  THE  SURFACE  OF  THE  SPHERICAL  BAFFLE 
IF(Un  CE  0  0  AND  U??  LT  O  0>  THEN 
PI  -  3  |4l439r6D4 
P2  -  0  0 


ELSElFlUt 1 
PI  -  0  0 

GC 

0 

0 

AND 

U22 

CE 

0 

0) 

THEN 

P2  •  0  0 
ELSEIFIUI 1 

LT 

0 

0 

AND 

U22 

LT 

0 

0) 

THEN 

PI  •  3  141392634 
P2  -  3  141392634 
ELSE 

Pi  •  3  141392634 

P3  •  0  0 
ENOIF 
RETURN 
END 
C 

. . . 

c 

SUBROUTINE  ADM(K. CLBO. CLAD. CRDA. CRAA. ADMl ) 

C 

C.  AON  evaluates  THE  ELECTRICAL  INPUT  ADMITTANCE  OF  THE  TRANSDUCER 
C  ADNl  IS  THf  COMPLEX  INPUT  AONfTTANCC  (OHNS> 

C 

COnnON  /IIP/  A.  B.L^LV.LZ.SRAD.ROD.COD.  ROH.  CON.  PI.  BN.  CD.  CON.  FO.  FF.  FD 
CONPLEX  AONt.CLBB.  CLAB.CROA.CRAA 

real  L.LY.LI.K 

ADNl  •  CNPLX<LZ*CON»(8-A)/LY.O  )  ♦ 

I  CNPLX(0  0.  -8N*OH*LI*CON*CON/<MCOD*ROO*LY)>*(CLBB-OLAB-CRBA*ORAA) 

RETURN 

END 

C 

C**##***##**** #•*•«•••« ••*•««««««•*«•••*••*••»#••»••#•*•••••••••••» 

C 

SUBROUTINE  RAOPRILY.LZ.ZR.VL.PR) 

C 

C  RAOPR  EVALUATES  THE  RADIATED  ACOUSTIC  POWER  AT  THE  APERTURES 
C 

CONPLEX  VL 

REAL  LY.L2.PR  v 

PR  •  <CA0S<VL)**r  )«ZR«LY*LZ  ' 

RETURN  \ 

END 
C 

c 

SUBROUTINE  SEN(FREQ.  Ul.  U2.  NUI.  HU2.  OXLNLA.  OXLNLB.  CXRLA.  CXRLB 
C 

C  SEN  EVALUATES  DTME  FARFIELD  PRESSURE  DIRECTIVITY,  GIVEN  THE  AMPLITUDES 
C  OF  THE  APERTURE  PARTICLE  VELOCITIES  AND  PHCSES 

C  2>TME  TRANSMITTINO  SENSITIVITIES  AT  ONE  NETER  PER  UNIT  AMP  AND  VCLT 
C  THE  FARFIELD  PRESSURE  IS  CALCULATED  USING  C  H  SHERMANS  RESULTS 
C 

l.CLBB.  CLAB.CRDA.  CRAA,  TCSL.  TVSL.  ANCMAX.FLAC.  FLAC2.  FLAG3) 

COMMON  /ZIP/  A.  D.L.LY.LZ.SRAD.  ROD.  COD.  ROM.  COM.  PI.BM.  CD.  CON 
l.FO.  FF.FO 

dimension  LO(I03>.L1(100.  100).  L3<  100.  100).  Jt(lOO) 

I. VI< 100).  WK(410>.P3(103> 

S.PSLdOS).  AN0<103) 

3.PN(I03) 

CONPLEX  CXLNLA.  GXLNLB.  OXRLA.  CXRLB.  GLBB.  CLAB.  GRBA.  GRAA 
1.  Pi  1  ( 103).  OH.  Cl  l'C22.  C3.  SUMI  SUN2.  SUNl  1.  SUN22.  P  1 .  P2.  Al.  A2 
INTEGER  N. NO. 1. U 

real  LO.  LI.L2.  MUI.NU2.  KA.  KR.L  LY.LZ.J1 
C 

C  N  IS  THE  NUMBER  OF  TERMS  REQUIRED  FOR  CONVERGENCE  OF  THE  SERIES  EXPRESSION 
C  FOR  THE  FARFIELD  PRESSURE  NO  IS  THE  NUMBER  OF  POINTS  CALULATED  ON  THE 
C  DIRECTIVITY  CIRCLE  R  IS  THE  RADIAL  DISTANCElMETERS)  IN  THE  FARFIELD  AT  WHICH 
C  PRESSURE  CALCULATION  IS  MADE 
DATA  N.  NO.  R/33,  73.  I  / 

NPTS  -  ND 

RO  •  SQRT(LV«LZ/PX) 

T  ■  ATAN(RD/SRA0) 

KR  -  2  •P|«FRtQ/CON*R 
KA  -  2  •PI*FRtO/CCM»SRAO 

Al  •  CMPLXtO  . (ROM«CON*UI/R>*SORT(SRAD*SRAD/(2  cPUKA))) 
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A?  •  CMPLX(0  <fiOn*CCf1*Ur/H»*SOBT<S«AD*SRAD/<2  •PI*kA»)» 

C  INjTlALItC  TMC  LECewDRE  POLYNO-II ALS 
LO(l>  *  0  0 
LO(?)  •  1  O 
LOO)  •  COSd) 

C  CALCULATE  The  PQLYNOfllALS  OV  RECURSION 
DO  I  I«4.N*S 
Rl  •  PL0AT<l-0» 

LOU*  -tr  ♦RI*S  >/<RX*l  )*L0<3)»L0<  I-l  »-RS ''IRS*!  >»LO<l-r) 

I  CONTINUE 

c  DA  IS  the  ancle  step  size 

DA  •  ?  •Pl/FLOAT(ND-l ) 

C  CALCULATE  TM£  ANCLE  DEPENDENT  LECENDRE  POLYNOMIALS 
DO  :  J>1.ND»: 
ancle  ■  FLOAT  (J-n*OA 
LH  I.  U)  •  1  0 

LI cr. j>  •  cosiancle) 

Lrd.s/)  »  1  0 
L2<r.U»  •  C0S<PI-ANCLE> 

DO  3  K>3.N 
RK  *  FlOATCK-;*) 

L1<K.U»*(2  •RK*I  >/(RK^l  )*LI<2.  J)*L1<K-|.  J>-RK/ <RK*1  )*Ll(K-r,J> 

LS<K.J)*<2  •RK*!  >/<RK*l  )  •LTC  2.  J)  •LTT  K- 1 ,  J> -RK/ <  RK*  1  )*L2(K-2.J) 

3  CONTINUE 

IF<FLAC3  EO  0  0)  CO  TO  20 

2  continue  ^  ^ 

C  MMDSUR  AND  MMDSYN  CALCULATE  THE  OESSEL  FUNCTIONS  OF  TmE  FIRST  AND  SECOND 
C  KINDS  RESPECTIVILY 
20  CALL  MMDSJR<KA.  5.  N*  1 .  J1 ,  WK.  I  ER  » 

CALL  MMDSYN<KA.  5,N*1.  Y1,  IER) 

Cll  -  CEXPICMPLX<0  0, -MUI)> 

C22  -  CEXP<CMPLX<0  0. -nu2>) 

C  CALCULATE  THE  N  DEPENDENT  TERMS 
DO  4  1-3. N*2 

C3  •  <FL0ATU)-3  ’/KA^Jl  <  |-2)-JI  ( 1-1  > 

C4  -  <FL0AT<I>-3  ) /KA*Y 1 < 1 -2) -Y I < I - I > 

C5  «  CEXP(CMPLX(0  0.-PI/2  •FL0AT<I-2>>» 

DH  •  CnPLX<C3.  C4) 

PlHl-2>  •  CMPLX(L0<  l-2)-L0(  I  >.  0  0.*C»/DM 
4  CONTINUE 
C  SUM  THE  N  TERMS 
DO  9  J"t,ND 
SUMll  ■  CMPLXTO  0.0  0) 

SUM22  ■  CMPLXTO  0.0  0) 

DO  6  !*t.N 

sum  •  PlHI)*CnPLX<Ll(I.U).0  O)  *  SUMII 
SUM2  -  Pini>*CMPLX!L2(l.  Jl.O  0)  *  SUM22 
SUMll  -  SUMl 

SUM22  >  SUn2 

6  CONTINUE 

c  calculate  the  angle 

ANCIUI  ■  (FL0AT(U>  -I  0)«DA»ie0  /PI 
PI  -  AMCIUSUMI 
P2  -  A2*cr2«SUM2 

C  SUM  the  pressure  fields  of  the  SOURCES  COHERENTLY  AND  FIND 
c  the  magnitude 

PS(U)  ■  CADS(P1  ♦  P2> 

PSL(U)  -  20  *AL0CI0<PS(U)/1  Oe-06) 

IF(PSL(J>  LE  90  0)  PSL(J>  -  90  0 
WR|TE<4. •)PSL<J) 

IF(FLAC3  EG  0  0>  CO  TO  30 
5  continue 

F(FLAC  EO  I  0)  CALL  D2PLT(NPTS. ANC.  PSL.  19) 

C  IF(FLAC2  EO  1  0)  CALL  D2PLT<NPTS. ANC. PS. 20> 

C  LOCATION  OF  THE  MAIN  DEAM  AXIS 
30  YMAXl  •  PS< I ) 

lF(rLAC3  £0  0  0)  CO  TO  10 

ANCMAX  •  1  0 

DO  7  I-2.NPTS 

YMAX  •  AMAXl (P5< I ).  YMAXl ) 

|F<YMAX  CT  YMAXl >  THEN 
ANCMAX  •  ANG<I> 

ELSE 

ENDIF 

YMAXl  -  YMAX 

7  CONTINUE 

C  USE  the  following  conditional  only  for  the  SYMMETRIC  SOURCE 
irCANCMAX  CT  90  J>  ANCMAX  -  160  -  ANCMAX 
C  CALCULATE  The  sensitivities 

10  TCS-  YMAX1/<CD*<0-A)«LY«1  4|4) 

TVS  •  YMAX1/<CD*LZ/C0N*1  414) 

TCSL  •  20  *4100. 0( TCS) 

TVSL  •  20  •ALOClO(rvS) 

RETURN 

END 

C 

. . . . . . . . 

C 

subroutine  D2PL T<NPT&. XX  YY  FLAG) 
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C  02^LT  IS  The  TUO  OinCNSICNAU  PLUTTINC  ROUTINC  USINC  AftLLlB 

c 

DlnENSICN  XX<tOOO». YY(IOOO) 

INTEGCP  PLAC.NPTS 

xniNi  •  xxcn 

XMAXt  •  XXCll 
YrtiNi  •  vvu> 
ynAxt  •  YY(i> 

DO  100  !»2.NP1S 

xniN  •  AniNKxxiD.xniNt) 

xniNi  •  xniN 

XRAX  •  ArtAXKXXdl.XriAXl) 

IKAXI  •  XriAX 

YfllN  •  AMfNKYYCD.VntNlI 
YniNt  •  vniN 

>nAx  •  AnAXHYVdi.YriAxn 
VHAXI  •  VriAX 
100  CONTINUE 

iFivnAx  CO  vniNi  ynin  *00 
ox  «  <xnAX  -  XHlNl/rLOAT<NPTS> 

XORO«i  0 
YORG-1  O 

xoCL«*x(iAx>xniN>/ro  o 
YOCL>(YnAX-VnlN>/10  0 
AYLEN-3  o 
AXLCN«0  O 
XTITL-AXLCN/2  0 
VTITL-AVLEN^O  ^ 

OX«<XrtAX>XniN)/AXLCN 
DY-(YnAX*>YnlNl/AYL£N 
CALI  PLTLFN<L“PL0T-) 

CALL  PLTOln<ll  .6  9.  1>4  ) 

CALL  PLTOftG(XOPC. YOPC) 

CO  TO  <1. 2.0.4. 9.6.  7.0.  10.  11>  12.  13.  14.  19.  16.  17.  10 

1  19. 20>  PLAC 

1  CALL  PLTAXIS<0  .0  .A¥LEN.90  .  YM|N.  VriAX.  YDEL 
l.l3MrtAC  PEFL  COEF.  13.  2.  -  1.-  1) 

CO  TO  200 

2  CALL  PLTAXISTO  .0  .AYLEN.  90  .  YNlN.  YHAX.  YDEL 
1.22NAP£RATUft£  PRESSURE < PA).  22,  2.-  1.  -  1) 

GO  TO  200 

3  CALL  PlTAXISTO  .0  .AYLEN,  90  .  YMlN,  YMAX.  VDEL 
1.29MAPER  PRESTOD  RE  IX10-6PA>,  29.  2.  -  1,-  1> 

GO  TO  190 

4  call  PLTAXI$(0  .0  .AYLEN,  90  .  YHIN.  YHAX.  VDEL 
I.14MAPER  PART  VEL.14.2.-  1,-  1) 

GO  TO  200 

9  CALL  PLTAXiSiO  .0  .AYLEN. 90  . YNtN. YHAX.  YDEL 
1.20H|nAG  ELECT  ADrKNHOS).  20. 2.  -  1.-  1) 

GO  TO  200 

6  CALL  PLTAXISTO  .0  . AXLEN, 90  . YHIN. YMAX. YDEL 
1.20HREAL  ELECT  AON < OHMS).  20.  2.  -  1.-  I> 

GO  TO  200 

7  CALL  PLTAXISIO  .0  .  AYLEN.  90  .  YMIN.  YMAX,  YDEL 
1.23HAC0US  PWRfDB  RE  1  WATT). 23.2.-  }.-  1) 

GO  TQ  I9G 

6  CALL  PLTAXISTO  .0  . AYLEN. 90  . YMIN.  YMAX.  YDEL 
1.23HREAL  INPUT  POWER  (WATTS) .  23.  2.  -  I.-  1) 

GO  TO  200 

9  CALL  PLTAXISTO  .O  . AYLEN,  90  . YMIN.  YMAX,  YDEL 

1. 27HREACT1VE  INPUT  POUER<UATTS).  27.  2.  -  1.-  1) 

GO  TO  200 

10  CALL  PLTAXIS<0  .0  . AYLEN. 90  , YMIN. YMAX.  YDEL 
1.29HPOW6R  EFFICIENCVTDD  RE  1).29. 2,  -  1.-  1) 

GO  TO  190 

11  CALL  PLTAXISCO  .0  .AYLEN.  90  ,  Yr’N,  YMAX.  YDEL 
1.24MP/V  SENlOe  RE  1  PA.  VOLT).  24  *•.  -  I.-  1) 

CO  TO  19P 

12  CAu  .tAXI8<0  ,0  .AYLEN. 90  .  YMIN.  YMAX.  YDEL 
l.rOH'-.^I  SEN\00  RE  I  PA/AnP).  23, 2,  -  l.-  l) 

GO  10  190 

13  CALI  PLTAXISCO  .0  .AYLEN.  90  ,  YMIN,  YMAX.  YDEL 
1  lOHANG  MAX  RCSP  (DEO).  10.2.-  1.-  I) 

GO  TO  200 

14  CALL  PLTAXIS(0  .0  .AYLEN, 90  . YMIN,  YMAX.  YDEL 
1.20HREAL  MUT  IMP  (RAVLS)  20.4.-  1,-  1) 

CO  TO  200 

19  CALL  PtTAXIS<0  .0  .AYLEN.  90  .  YMIN.  YMAX.  YDEL 
l.rOMlMAO  MUT  imp  (RAYLS).20  4.-  I,  -  I) 

CO  TO  200 

19  CALL  PL»AXIS(0  .0  .AYLEN. 90  .YMIN.  YMAX  YDEL 
1  19HREAL  FM  IMP  (RAVLS) .  1 9.  4 .  -  1.-  1) 

GO  TO  200 

17  CALL  PLTAXIS«0  0  . AYLEN, 90  YMIN. YMAX  YDEL 
1  19HtMAG  FM  imp  <RAYLS).19. 4  -  1.-  1) 

GO  TO  200 

10  CALL  PLTAXIS(0  0  AYLEN. 90  YMIN. YMAX.  YDEL 
1. 21HRATI0  MA6<Zn)/nAG(Z3).  2!  4.-  l.-  l) 

GO  TO  200 

19  CALL  PlTAXISCO  0  .AYLEN.  90  .  YMIN.  YMAX.  YDEL 
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20 

ISO 

CISO 

C 

200 

250 

300 


14.4,  -  1,  . 


1 ) 


AYLcN, 90  . YMIN,  YMAX,  YDEL 
<PA),  19,  4,  -  1,  -  1 ) 

AXLEN,  0  ,  XMIN,  XMAX,  XDEL,  UHLOGIO(FREQ), 


0  .AXLEN,  0  .  XMIN,  XMAX,  XDEL.  4HRA/R, 


0  .AXLEN, 0  , XMIN,  XMAX, XDEL,  14HFREQUENCY  (HZ) 


1, 14HPRES  LEVEL (OD) 

GO  TO  250 
CALL  PLTAXIS(0  ,0 
1, 19HPRES  AMPLITUDE 
GO  TO  250 
CALL  PLTAXIS(0  ,0 
1-11,4,  -  1,  -  1 ) 

CALL  PLTAX1S(0 
1-4.4,-  1,-  1) 

CO  TO  300 
CALL  PLTAXIS(0 
1-14.  4,  -  1,  -  I  > 

GO  TO  300 
CALL  PLTAXIS(0 
1.11 MANGEL  (DEG) 

CO  TO  300 

CALL  PLTAXIS(0  ,  AYLEN,  AXLEN,  0  ,  XMIN,  XMAX,  XDEL,  LABX,  0,  0, - 
CALL  PLTAX1S(AXLEN,  0  .  AYLEN,  90  ,  YMIN,  YMAX,  YDEL,  LADY,  0,  0, 
CALL  PLTDATA(XX.  YY,  NPTS,  0,  0,  XMIN,  DX,  YMIN,  DY,  0  08) 

CALL  PLTLINE(XTITL,  YTITL,  -  1) 

CALL  PL  TEND! 11  0,  8  5) 

RETURN 

end 


0  ,  AXLEN,  0  , 
-11,4,-  1,  - 


XMIN,  XMAX,  XDEL 
1 ) 


1,  - 
1. 


1) 
1 ) 


PROGRAM  THERMO! INPUT,  OUTPUT,  DATA,  TAPE2=INPUT,  TAPE4=DATA 
I,  PLOT,  TAPE5«=PL0T) 

C 

C**«*  *«*•***«  »#««»*««*«**«*«««**«««*«««**««««»»* 

c 

c  PROGRAM  THERMO  PERFORMS  A  FREQUENCY  DOMAIN  ANALYSIS  OF  THE 
C  THERMOACOUSTIC  TRANSMITTING  CHARACTERISTICS  OF  A  PLANE  WAVE  MODE  MHD 
C  TRANSDUCER  THERMO  COMPUTES  THE  ACOUSTIC  FIELD  AS  A  FUNCTION  OF  THE 
C  THE  TRANSMITTED  ACOUSTIC  FREQUENCY,  WHICH  IS  TWICE  THE  ELECTRICAL  DRIVE 
C  FREQUENCY  THE  POWER  EFFICIENCY  AND  TRANSMITTING  SENSITIVITIES  ARE 
C  COMPUTED  AS  A  FUNCTION  OF  THE  ELECTRICAL  DRIVE  FREQUENCY 
C 

C  THERMO  GENERATES  GRAPHIC  OUTPUT  FOR  THE  FOLLOWING  TRANSDUCER 
C  CHARACTERISTICS 
C 

L  I)  MAGNITUDE  OF  THE  REFLECTION  COEFF  (MRNL)  VERSUS  ACOUSTIC  FREQUENCY  (AFREQ) 

C  2)  MAGNITUDE  OF  APERTURE  PRESSURE  (RPNL)  VERSUS  AFREQ 
C  3)  DB  LEVEL  OF  RPNL  VERSUS  LOGIO(AFREQ) 

C  4)  "  "  OF  APERTURE  PART  VEL  (RVNL)  VERSUS  LOG (AFREQ) 

C  5)  "  "  OF  radiated  ACOUSTIC  POWER  (PR)  VERSUS  LOG( AFREQ) 

C  6)  “  "  OF  TRANSDUCTION  POWER  EFFICIENCY  (EFFIC)  VERSUS 

C  L0G10(ELECTRICAL  SIGNAL  FREQUENCY)  (FLOG) 

C  7)  QUADRATIC  CURRENT  TRANSFER  »'UNCTiON(QCTF)  VERSUS  LOGFREQ 

C  8)  DB  LEVEL  OF  QUADRATIC  VOLTAGE  TRANSMITTING  SENSITIVITY(QVTS)  VERSUS  LOGFREQ 
C  9)  PRESSURE  LEVEL  DIRECTIVITY 
C  lOANGLE  OF  PRIMARY  BEAM  VERSUS  AFREQ 
C 

C  DATA  INPUT  IS  THROUGH  THE  BLOCK  DATA  STATEMENT  AT  THE  END  OF  THE 
C  MAIN  PROGRAM 
C 

C******»*«****»*»****««*#«**»*«#*#*-»***###**#»**i,#*»*##»#»*##**«## 

C 

COMMON  /ZIP/  A,  B,  L.  LY.  LZ.  SRAD,  ROD,  COD,  ROM,  COM,  PI,  BETA.  CP,  CD 
1  CON,  BM.  FO,  FF.  FD 

DIMENSION  AFREQ( 1000),  EFREO( 1000),  FLOC( 1000),  PR( 1000) 

1.  P5PL( 1000), RPNL( 1000),  RVNL( 1000),  MRNL( 1000),  FLOGA( 1000) 

2.  TCSL( 1000),  TVSL ( 1000).  ANOMAX ( 1000) ,  EFFIC ( 1000) 

3.  QTSL( 1000).  QCTFR( 1000).  QCTFI ( 1000) 

COMPLEX  RNL( 1000),  RL ( 1000), PNL ( 1000),  PL( 1000) 

1.  VNL( 1000),  VL( 1000),  QCTF( 1000) 

2.  GLNLA,  GLNLB,  GRLA,  ORLB 

3.  GXLNLA,  CXLNLB,  CXRLA,  GXRLB 
REAL  K,  L,  LY.  LZ,  MRNL 
INTEGER  !.N 

N  s  (FF-PO)/FIi  ♦  1 
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i 

1 

} 


c 

C  PRIriARV  CALCytATlON  LOOP 
C 

DO  1  l-l.N 

C  CALCULATE  the  wave  NUMOER  AND  FREQUENCY 

K  «  (T  •PI*F'0  ♦  FLOATCI-l)*?  •PI«FD)/COO 
AFREOd)  •  n*C00/<2  •Pl> 

EFREO(l)  •  AFREQd)/? 

FLOCA(I)  •  ALOCIO(AFRCO( 1 ) ) 

FLOOCI)  •  ALOC10<AFREO<t)/S  ) 

C  REFL  EVALU7ES  TmE  APERTURE  REFLECTION  COEFFICIENT 
CALL  R£FL(K.  RNL(  I  >.  RL<  I  >*  nRNL<  1  )•  2FnA> 

C  GREEN  EVALUATES  SEVERAL  COMPLEX  COEFFICIENTS 

CALL  CR£EN(K. A. 0. L. RNL < I > . RL < I ) . GLNLA. CLNLD. CRLA. CRLO 
1 . CXLNLA. GXLNLB.  CXRLA. CXRLO ) 

C  PRES  evaluates  THE  APERTURE  PRESSURE  IN  THE  WAVE  GUIDE 
CALL  PRESIGLNLA.  CLNLB. CRLA. GRLD. PNL < I > . PL ( I ) > 

RPNL< I )  •  CA0S(PNL( I ) ) 

PSPL<I»  •  ZO  •AL0GI0<RPNL< I )/  000001) 

C  VEL  evaluates  TmE  APERTURE  PARTICLE  VELOCITY 

CALL  VEL<K. CXLNLA. CXLNLD.  CXRLA,  CXRLO.  VNL( I ) .  VL < I ) • U1 .  Ur.  P 1 .  PD) 
RVNL(I)  •  DO  *ALOCIO<UD) 

C  WRITECA  *>01  P)  UD.PD 

c  OTF  IS  THE  Quadratic  current  transfer  function  routine 

CALL  OTFiK.  RNLd  >.  RLd  ).  0CTF<  I  )  ) 

OCTFRC I )  •  real  »OCrF( I ) ) 

OCTFKI)  •  AI.*1AC<0CTFM  ) ) 

C  RADPR  evaluates  the  RADIATED  ACOUSTIC  POWER 
CALL  radprily  l2  :fmr. VL( I ),PRR) 

PR(I)  •  DO  •AlOCtO(PRR) 

c  PE  IS  the  average  real  input  electrical  power 
PE  •  CD*C0*(B-A)*LY*LZ/«r  *C0N) 

C  EFFIC  IS  The  transduction  process  power  efficiency 
EFFIC<n  •  DO  •ALOCIOCPRR/PE) 

IF<£FFIC<I)  LE  -000  )EFFIC<I)  -  -300 

c 

C  OVTS  EVALUATES  THE  PRESSURE  RADIATION  PATTERN.  QUADRATIC  TRANSMITTING 
C  VOLTAGE  SENSITIVITY  AND  THE  ANCLE  OF  THE  MAIN  LODE 
C 

C  IF  A  DIRECTIVITY  PLOT  IS  DESIRED  SET  FLAG  •  I  0  <00  SCALE) 

C  IN  IF  STATEMENTS  OElOW 

C  AND  /  OR  SET  FLACD  ■  I  0  ILINEAR  SCALE) 

C  SET  FLAG  AND  FLACD  *00  OTHERWISE 
C 

C  CO  TO  I 

FLAG  -  0  0 

flags  >00 

IF(I  EO  1)  FLAG-0  0 
IF( I  EG  1 >  FLACS-0  0 
C 

C  IF  THE  X  AXIS  SENSITIVITY  IS  DESIRED  SET  FLACD  -  0  0 
C  IF  THE  SENSITIVITY  ON  THE  MAIN  DEAM  AXIS  IS  DESIRED  SET  FLAGS  •  1  0 
C  FLAC3  must  OE  I  0  WHEN  CALCULATING  DIRECTIVITY  PATTERN'S 
C 

FLAC3  -00 
C 

CALL  QVTS(AFREO< I ). AOS (U1 >.  A0S(UD ).  P 1 .  PD 
1 .  OTSL  (  1  > .  ANGMAX  (If.  FLAG.  FLACD.  FLAGS ) 

IF(0TSL<I)  LE  -105  )OTSL<I)— 105 
WR1TE(4,*)  EFREQ(  I  ).OTSL(  1  ) 

1  CONTINUE 
C 

C  DDPLT  PRODUCES  A  NEUTRAL  PLOT  FILE  FOR  PLOTTING  ON  THE  TEKTRONICS 
C 

CO  TO  D 

CALL  DDPLT (N. AFRCO.MRNL. 1 ) 

CALL  ODPLTIN.  AFREO.RPNL.D) 

CALL  ODPLTIN.  FLOCA,  PSPL.  3) 

CALL  DDPlT  (N.  FLOCA.  RVNL.  4  ) 

CA.L  ODPLKN.  FLOCA.  PR.  5) 

CALL  ODPLT(N.FLOG.EFriC.6J 
D  CALL  DDPLT<N.FlOG  OTSL.  7) 

GO  TO  3 

CALL  DDPlTin.  EFREQ.  OCTFR.  8) 

CALL  DDPLTCN.  EFREO. OfTFI,  ID) 

CALL  DDPLT  (N.  AF  REG.  ANGMAX.  9) 

3  STOP 
END 
C 

. . . . . . . . . . 

c 

OLOCK  data 
C 

c  variadle  inputs  to  program  Thermo 
c 

C  A  left  DOVNOARY  op  SOURCE  VOLUME  (MEIERS) 

C  0  RIGHT 

C  L  total  WAVEGUIDE  LENGTH  (METERS) 

C  LV  transverse  WAVEGUIDE  LENGTH  v-DIR  (METERS) 

C  ll  2-niR 
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C  SRAO  RADIUS  OF  SPHERICAL  DAPPLE  fhCTCRSt 
C  ROD  fluid  density  IN  UAVE  CUIDC  CKC/n£TeR*«3) 

C  COO  -  SOUND  SPEED  “  *  (METERS/SEC) 

c  ROM  Fluid  density  in  free  medium  <kc/m£ter»*3) 

C  COM  -  SOUND  SPEED  "  "  <METCRS/SEC> 

C  CD  CURRENT  DENSITY  AMPLITUDE  < AMPERES/METER**?) 

C  OETA  thermal  COEFFICIENT  OF  EXPANS10r4  OF  FLUID  (M*«3/M**3/HELVIN) 
C  CP  SPECIFIC  HEAT  CAPACITY  OF  FLUID  < JOULCS/KC/KELVIN) 

C  CON  electrical  conductivity  of  FLUID  IN  WAVEGUIDE  <I/OHMS  METER) 
C  FO  INITIAL  CALCULATION  FREQUENCY  (HERTZ) 

C  IF  final  ■  “  “ 

C  FD  FREQUENCY  STEP  SIZE  (HERTZ) 

C 

COMMON  ^Z IP/A.  B.  L.  LY.  LZ.  SRAD.  ROD.  COD.  ROM.  COM.  PI .  OETA.  CP.  CO 
l.COM.  DM.  FO.  FF.FO 

real  L  ly.lz 

DATA  A.D  L  LY.LZ.SRAD.ROD.COD.ROM.COM.  PI.  BETA.  CP.  CO.  CON 
1.  DM.  FO.  FF.  FD/>0  OS.  0  OD.  lO-  0316.  0301.0  032.1041  .1600  .998 
2.1461  3  141S9.0  00021.3966  .14436  .3  28.  40 

3. 500  . 13000  230  / 

C 

C  PHYSICAL  PROPERTIES  OF  NACL  6X  SOLUTION  AT  30  C 
C 


c 

CON 

5  26 

(1 /OHM/M)  MEASURED 

c 

ROD 

1041 

(KC/M**3) 

c 

COO 

1500 

(M/SEC) 

c 

CP 

3068 

( joules/kg/kelvinj 

c 

beta 

2  1C 

04  (M**3/M**3/'K£LVIN) 

C 

C  PMVSCIAL  PROPERTIES  OF  MERCURY  AT  20  C 
C 

C  CON  1020000  (l/OMMM) 

C  ROD  13600  (KG/M»*3) 

C  COO  1450  (M/SEC) 

C  CP  140  (J/KG/K) 

C  BETA  1  83E-04  (1/K) 

C 

END 

C 

. . . . . . * . . . . . . 

C 

SUOROUT INE  REFL ( K. RSNL.  RSL. RSNLM. ZSR ) 

C 

C  SUBROUTINE  REFL  CALCULATES  THE  COEFPICENTS  AND  EVALUATES  THE 
C  QUADRATIC  EQUATION  FOR  THE  REFUECTION  COEFFICIENTS 
C 

COMMON  /  Z IP/A.  D.  L.  LV.  LZ.  SRAD.  ROD.  COD.  ROM.  COM.  PI .  BETA.  CP.  CD 
l.CON,  DM,  FO.FF.FD 

COMPLEX  KA.  KO.  KL2.  ZS.  ZM.  Al.  Bl.  Cl.  Dl.  A2.  BS.  C?.  02.  E2.F2.G3.K2 
1  C 1 1 .  C22.  C33.  C44.  C33.  C66.  RSNL.  R8L.  RLL.  RL.NL 
REAL  L.  LY.LZ  K.PC 
PC  ■  ROO*CCD 
KA  -  CMPLXIO  0. K*A) 

KD  ••  CMPlX(0  0.  K*D) 

KL2  «  CMPLXIO  0.K*L/2) 

C  CALCULATE  the  APERTURE  IMPEDANCES 

CALL  M1SP(K*C0D. 2MR. ZMI.  ZSR. ZSI) 

ZS  -  CMPLX(ZSR. ZSI>/CMPLX(PC.O  0) 

ZM  -  CMPLX(ZMR. ZMI)/CrtPLX<PC.O  0) 

C  calculate  the  QUADRATIC  EQUATION  COEFFICIENTS 
Cl  -  CCXP(KA)  >  CEXP(KB) 

Dl  -  CEXP(-KA>  -  C£XP(-KB) 

Al  •  CEXP(-KL;') 

Bl  *  CEXP(KL2) 

A?  •  B1*D1*(CMPLX( 1  0  0>-ZS)»Dl*Cl*ZM 

B2  •  A1«D1*(CMPLX( 1  0  0 > ♦ZS > 'B 1 *01 vZM 

C2  ■  01*C1*(CMPLX(-1  .0  0)*ZS>-A1*C1*ZM 

Dr  ■  -Al*Cl*(CMPLXI 1  .0  0>*2S»*A1#D1*ZM 

E2  •  01*Cl*(CMPLX(t  .0  0>>ZSMOl*Dl*Zn 

F2  ■  0I*01*<CMPLX<-1  0  0)*ZS»-A1*D1*ZM 

02  •  A1»C1»(CMPLX(1  .0  0>*ZS)*D1*C1*ZM 

M2  -  ‘AWDMICMPlXI  1  .0  0)*Z8>*A1*C1»ZM 

Cll  •  A2*F*2  -  D2*E2 

C22  •  A2*M2  -  D2*Gr  ♦  C2*F2  -  D2*E2 

C33  ■  C2*M2  -  02*C2 

C44  •  Cr*£2  -  A2*C2 

C35  ■  D2*£2  ♦  C2*F2  -  32*02  -  A2*H2 

C66  •  D2*F2  -  D2*H2 

C  CALCULATE  the  COMPLEX  ROOTS  USING  THE  FOLLOWING  IMSL  ROUTINES 
C  THE  ROOT  with  THE  MAGNITUDE  LESS  THAN  ONE  IS  PHYSICALLY 
C  MEANINGFUL  The  OTHER  IS  DISCARDED 

CALL  ZQADC(C33.C22  C  11 .  RSNL ,  RLNL.  IER  ) 

CALL  Z0ADC(C66.  C55.  C44.  RSL.  RLL.  IER) 

RSNLM  •  CABS(RSNL> 

RETURN 

END 

C 


c 

c  «;sp  evaluates  the  sect  and  mutual  acoustic  impedances  ir  the 
C  WAVEGUIDE  APEPATUflES  THE  RESULTS  ARE  INPUT  TO  REEL  TO  EVALUATE  THE 
C  REPLECTION  COEFFICIENTS 

C  MISP  USES  THE  EQUATIONS  DEVELOPED  OY  C  H  SHERMAN  TO  EVALUATE  THE  MUTUAL 
C  AND  SELF  impedances  OF  TWO  CIRCULAR  SOURCES  180  DECREES  APART  ON  THE 
C  SURFACE  OF  A  spherical  BAFFLE 
C 

DIMENSION  P(  1000)UK(S000>.  JtC  1000).  Yt  UOOO). 

IPHIOOO) 

COMMON/ZIP/A.  B.  L  LV.  LZ.  5RAD.  ROD.  COD.  ROM.  COM,  PI.  BETA.  CP.  CD 
I.  CON.  Bn.  FO.  FF.FO 
real  Jl.  Yl.L.LY  LZ 
INTEGER  I.H.  N 

C  N  IS  THE  NUMBER  OF  TERMS  IN  THE  SUM  REQUIRED  FOR  CONVERGENCE  AT 
C  LARGE  VALUES  OF  R«RD 
N  -  30 

C  RD  IS  THE  EFFECTIVE  RADIUS  OF  THE  CROSS  SECTION  OF  A  RECTANGULAR 
C  WAVEGUIDE 

RD  -  S0RT(LY«1  Z/PP 
T  ■  ATAN(RD/SRAD> 

C  ■  ROn*COn*SRAD«SRAD/(RO«RD; 

C  INITIALIZATION  OF  LtCtNCRE  POLYNOMIALS 
PU )  •  1  0 
p<?)  •  1  o 
P<3)  ■  COS(T) 

PI < 1 )  ■  1  0 
Pl(2)  •  -1  0 
M  • 

C  CALCULATION  OF  LEGENDRE  POL>N0MIALS  BY  RECURSION 
DO  1  K-O.M 
Ri  •  Float <K» 

RD  •  R|  -  ■  0 

P(K*1)  ■  P«3»*P<K>#(D  •RPl  >/<RP!  )-P(K-l  )»R|/<RP1  ) 

P|(K)  ■  PP2)*PI  (^-1  )*(2  •RD*  )/<R2*l  >-PPK-D>»R2/lRD*I  > 

1  CONTINUE 

C  The  BESSEL  FUNCTION  EVALUATIONS  ARE  PERFORMED  BY 
C  MMOSJR  AND  MMBSYN  (IMSL  ROUTINES) 

I  ■  W*SRAD/COM 

CALL  MMBSJRIY.  3> N«3. JI .  UK.  lEP > 

CALL  MMSSYN(X.  3.  N*3,  YI .  IER  ) 

ZRNI  *00 
ZCMl  -  0  0 
ZRSI  ■00 
ZCSI  *00 
ND  a  N»1 

C  THE  SERIES  ARE  SUMMED  IN  THE  |  LOOP 
C  TO  CALCULATE  THE  SELF  AND  MUTUAL  IMPEDANCES 
DO  3  lai.NO 
SI  a  FLOAT! P 

Da  ((<SI-I  )/X)«JMtPJt(PI>>*«2  ♦ 

{(((Sl-i  )/X)«Yl(l)-Yl<Pl))*»D 
D1  -  ((Sl-1  )/X>*(Jl(l)*Jl<n«YI(I)»Yl(I))  - 
PJl  (  I)»U1  (  PP«Yt  ( 1  )«YI(  PI  ) ) 

DD  -  -UI  <  I  >*V1  !PP*J1  U*1  >*VP  I ) 

RS  -  D2/D 
CPL  “  -Dl/n 

CDM  -  <1  /<2  *<51-1  )♦!  ) )*< <P< I )-P< I ♦?) )««2)*P1 < I ) 

CDS  -  (1  /«2  *<SI-1  >♦!  >)*<P(I>-P<1*2))*»2 

ZRMl  -  GrM*RS»ZRMl 
ZCMl  -  CrM*CPL*ZCMI 
ZRSI  -  CSS*RS^ZRS1 
ZCSI  -  CDS*CPL*ZCS1 
3  CONTINUE 
c  result 

RM  -  C*ZRM1 

CXM  •  C«ZCM1 

R  ■  C*ZRS1 

cx  •  C»ZCSI 

RFTURN 

END 

C 

c**«*«*  . . . . . . . . 

c 

Subroutine  creenjk  a. b.l  rnl  rl. clnla. clnlb. crla.  crlb 

1>  GXLNLA.  GXLNLB.  GXRLA.  CXRLB) 

complex  CC.KA.  KD  KLr.RNL  RL  CUNLA.  CLNLB.  CRLA  CRLB 

1  LA  lb  RA  RB 

2  GXLNLA,  GXLNLB  GXRLA,  CXRLB 
real  k  L 

KA  -  CMPL  X  *0  ,  K*A/ 
aB  •  CMPLX(0  K*D' 

KLD  •  CMPLX<0  K*l  /:  > 

CC  •  CMPLX<1  0  0'/<CMPLXir  •K  0  >*<CMPLX<1  0  )-RNL»RL>> 

LA  •  CC*<CCXP<ktA>  <RU*CeXP(-KA)  ) 

RA  -  CC*<CEXP( -KA) <RNL*CEXP4KA> i 
lb  •  CC»<CEXP(K0  )-RL*C£xP(-kO) » 

SB  -  LC*<CEXP( -KB  > -*RT4LACEXP(K0  >  ) 

Clnla  •  CMPLXJO  O  -l  0)»LA*KEXP(KL2)*RNL»CEXP<-KLr)  > 

GLNlO  •  CMPLXtO  0  -I  0>*LO*<CCXP<KLD)*RNL*CEXP<-KLD) ) 
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CRLA  «  CHPLXCO  0.  t  0>«RA*<CCXP  (KL2)  *RL»CEXP  ( -'HL2 )  ) 

CRLD  •  CnPtX<0  0  t  0)*RO*(CexP<KL2l*«L*C£XP(-KL?> > 

CXLNLA  •  C«PLX<0  0*1  0)«LA*(C£XP<AUr>  -  RNL*CEXP( -KL2) > 

CXLNLO  ■  CHPLXfO  0.1  0)»LO*«CEXP«KU2*  -  RNL*CEXP< -KL?) > 

CXRLA  >  CMPLXtO  0.1  0>*RA* CCEXP f KL2)  **  RL.*CEXP ( -KL2)  ) 

CXRLO  •  CnPLXlO  0.1  0)*RO«<C£XP<HL2)  -  RL *CrxP ( -KL2 M 

RETURN 
END 
C 

. . . . . . . 

c 

SUOROUT  INE  PRES<CLN1.A.  CLNLD.  CRLA.  €RL0>  PNL-  PL  > 

C 

C  PRES  EVALUATCS  TmE  APERTURE  PRESURE 
C  PNL  IS  THE  complex  PRESSURE  AT  TmE  LEFT  APERTURE 
CPL--  -  --  -  RICHT 

C 

COMMON  /ZIP/  A.  o.L  LT.LZ.SRAO.  fiOD.CCD.  ROM.  con.  PI.  PETA.  CP.CO 

1.  CON.  DM.  PO.  ff.FO 

COMOLE  X  PNL .  PL.  Cl  NLA.  CLNLO.  CRLA.  CRufi 

PNL  ■  CMPLXi -QETA*CD*CO*COO/<2  •C0N*CP>.0  )»(CLNLA  ►  CLNLO > 

PL  ■  CMPLX<-OETA*CO»CO*COD/<2  *C0n*CP>.0  >*(CRLA  -  CRLB) 

RETURN 

End 

c 

. . .  . . . . . 

c 

SuOROUTInE  v£L(K,  CxLNLA.  CXLNLO.  CXRLA.  CXRLD.  VNL.  VL  UI.U2.PI.P2> 

c 

c  vel  Evaluates  the  particle  velocity  at  the  aperture 
C  VNL  IS  the  complex  particle  velocity  at  the  left  aperture 

C  VL  -  -  -  •  -  •  RICH! 

c 

COMMON  /ZIP'  A.p.t.LY.LZ.SRAO  ROD.  COD.  ROM.  COM.  PI .  DCtA.  CP.  CD.  CON 
I  BM.FO.FF.PD 

complex  VL  VNL.CXLNLA,  CXLNLO.  CXRLA,  CXRLO 
real  k 

VNL  •  CMPLX<-DeTA*CD*C0/<2  •CON«CP*HOD ) .  0  0)»<CXLNi.A  -  CXLNLO) 

VL  •  CMPLX<-OETA*CO*CO/<2  ♦C0N«CP*R0D>,0  0)*(CXRUA  -  CXRLB) 

Ut  •  CAQS(VL> 

U2  ■  CAaS<VNL) 

UU  •  REALIVL) 

022  •  ReAL<-VNl > 

C  the  conditionals  OElOW  evaluate  THE  RELATIVE  PHASES  OF  THE  APERTURE 
c  particle  velocities  referenced  to  the  Surface  of  the  spherical  daffle 

IF<un  CE  0  0  and  U22  LT  0  0)  THEN 
PI  •  0  0 

P2  •  3  1414992654 

ELSElF<un  CE  0  0  AND  U22  CE  0  0>  THEN 
Pi  -  0  0 
P2  •  0  0 

ELSEIF<un  CT  0  0  AND  U22  LT  0  0)  THEN 
PI  •  3  141992694 

P2  •  3  141992694 
ELSE 

PI  •  3  141992694 
P2  •  0  0 
ENOIF 
RETURN 

End 

c 

«•••«•*««•«»••*•«**«#••**•*«««•«•»«»•••»••••»••••• •#«•»•••* 

c 

subroutine  QTF<k,RNL.RL.OCTFj 

COMMON  /ZIP  A,  0  L.LV  LZ.SRAD.  rod.  cod.  ROM,  com.  pi,  ofta.  cp.cd 
1.  con.  dm.  FO-  f*“  fd 
real  K  LY  I z 

COMPLEX  X1.X2.  X3  X4.X9.X6.RNL  RL  OCTF 

XI  ■  CMPlX< -BETA#Dfl«C0N*C0N*LZ/<4  •CP*ROD«LY»L  <«K»K ) .  0  0>/ 
t(CMPLX<|  0.0  0>-RNL«RL) 

X2  ■  CMPLXiO  O. 2  ♦K*B> 

X3  ■  CMPlX'O  O  2  •K^A) 

X4  •  CMPlX<0  0.K*<A*D>> 

X9  •  RNL«<CexPiX2>^CEXP<X3)'CMPLX<2  0  )«CEXP<X4)) 

X6  ■  RL*«CEXP< -X2» -CEXP< -X3> •CMPLX<2  0  )*CEXP(-X4)> 

OCTF  •  Xl*tX9*X6) 

RETURN 

END 

C 

c*««««****««»**t*»«**«« 

c 

SuOROvTJnE  RADPP(i.Y  iz  ZR  VL  PR) 

C 

C  RADPR  Evaluates  tmc  raduted  acoustic  power  at  the  apertures 

c 

complex  VI 
real  IV  lZ.  PR 

PR  •  <CA05<vl>««2  )*ZR*LV*lZ 

RETURN 

END 
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c 

c 

c 


subroutine  QVTS<FflEO.  Ul.  TVSL.  ANCMAX,  FLAG-  FLAGS.  FLA03> 


C  OVTS  EVALUATES  DTME  FARFIELD  PRESSURE  DIRECTIVITY.  GIVEN  THE  AMPLITUDES 
C  OF  THE  APERTURE  PARTICLE  VELOCI’IES  AND  PMESES 

C  2»THE  transmitting  SENSITIVITIES  AT  ONE  METER  PER  UNIT  AMP  AND  VOLT 
C  TME  FaRFIElO  pressure  IS  CALCULATED  USING  C  H  SHERMANS  RESULTS 

COMMON  /ZIP/  A.  0.  L.  LY.LZ.  SRAD.  ROD-COD.  ROM.  COM,  PI.  BETA,  CP.  CD 

l.CON.  BM.FO.FF.FD 

DIMENSION  LO<  109).  LI  (  100.  100).  L2<  100.  100).  JH  100) 

1.  Y1  (  100).  WK(410)>  PS<  109) 

S.PSL<109).  ANC<109> 

3.PN<109) 

COMPLEX  Pi  1  <  109).  Cl  I.  C22.  OH.  C9.  SUMI.  SUMS.  SUKl  1.  SUMSS.  PI.  PS.  Al,  AS 
INTEGER  N.  NO.  1.  J 

REAL  LO  LI.LS.  MUl.MUS.KA.KR.  L.  LY.LZ.  Jl 

C  N  IS  TmE  number  of  TERMS  REQUIRED  FOR  CONVERGENCE  OF  THE  SERIES  EXPRESSION 
C  FOR  the  FARFIELO  PRESSURE  ND  IS  THE  NUMBER  OF  POINTS  CALULATED  ON  THE 
C  DIRECTIVITY  CIRCLE  R  IS  THE  RADIAL  DISTANCE  (METERS)  IN  THE  FARFIELO  AT  WHICH 
C  PRESSURE  calculation  IE  MADE 
C 

DATA  N.  NO.  R/99.  73.  I  / 

NPTS  -  NO 

RD  •  SQRT<LY#l2/PI) 
r  •  atan(Rd/srad> 

KR  •  2  •PI*FR1-Q/C0M*R 
KA  ■  2  •P I •FRtO/COM*SRAO 

Al  ■  CMPlx<0  0.  (ROM*COM«UI/R)*SORTISRAD»SRAD/(2  •PMKA))) 

A2  •  CMPlX'O  0.  <R0M«C0n«U2/K>»SQRT<SRAD«SRAD/<S  •Pl^KA))) 

C  initialize  the  LEGENDRE  POLYNOMIALS 
LOM  )  -00 
L0<2)  •  1  0 
L0<3)  •  COS<  T) 

C  CALCULATE  THE  POLYNOMIALS  OY  RECURSION 
00  1  1«4.N*2 
RI  •  FLOAT< I-3> 

LO<I)  -(S  ‘RI^l  )/<RI»l  )«L0(3)«L0<I-I)*RI/<RI*I  )*10(I-S> 

1  CONTINUE 

C  DA  IS  the  angle  S'^EP  size 
DA  -  2  #P:/riOAT<NO-l ) 

C  CALCULATE  The  a.yGLE  DEPENDENT  LECENDflE  POLYNOMIALS 
DO  2  J«1 

angle  •  PlOAT(U-1 )«0A 

c  angle  •  p:/2 

LMl.U)  •  ?  0 

L1<2.J)  •  '.OS(ANGLE) 

L2<1.  J)  •  ;  0 

l2<2. -j)  <  cas<Pi-ANGLe> 

DO  3  N 

RK  •  FL0A’i*\-2> 

Li<K.U)»<«;  •RK^l  )/«RK*l  )*LI(2.  J)*LI<H»I.  J>-RK/(RK»I  )«LMK-S.J) 
L2<K.U)-'2  )/<RK*I  )  •LSI  2.  J)  •L2<K-t .  J/-RK/ ( RK+ 1  )«L2(K-2,J) 

3  CONTINUE 

|F<FLAG3  ro  0  0)  GO  TO  20 
2  CONTINUE 

C  MM0SUR  AND  MMOSTN  CALCULATE  THE  OESSEL  FUNCTIONS  OF  TME  FIRST  AND  SECOND 
C  KINDS  RESPECTIVILY 
20  CALL  MMOSUR<KA.  9- N* I . Jl .  UK.  lER ) 

CALL  MMOSYN<KA.  9. N* I , Y 1 . lER ) 

Cll  ■  CEXP(vMPLr<0  0  -MUD* 

CSr  -  CEXP(CMPLX<0  0. -MUD) 

C  CALCULATE  THE  N  DEPENDENT  TERMS 
DO  4  I»3.N*2 

C3  •  <FL0AT<I)-3  >/KA*UHI-r)-JI(I-D 
C4  -  (FlOAT(I)-3  )/KA*YI<l-r)-YI(I-D 
DM  •  CMPLX<C3-  C4> 

C9  ■  C£XP(CMPLX<0  0.  -PI/2  •FLOATS  1-2) ) ) 

Pll<l-2>  ■  CMPLXt<t0(I-2>-L0<l)).0  0)*C9/DM 
4  continue 
C  SUMMATION  LOOPS 
DO  9  J>).ND 
SUMl 1  •  CMPLX(0  000) 

suMrr  »  CMPLX(0  o  o  o> 

DO  6  I>l  N 

SUMl  -  PI  K  D*CMPlx*lD  I.  J)- 0  O)  ♦  SVMIl 
SUMS  •  Pi  n  D»CMPLX<LrM  J>  0  O)  ♦  3UM22 
SUMll  ■  SUMl 
SUM2r  •  SUM’ 

6  continue 

C  calculate  The  FARFIElD  ANCLE 

ANGJJ)  •  <FlOAT(J>  -I  0>*DA«teo  /PI 
PI  •  Al -SUMl •Cll 
P2  •  A2«5UM2*C22 

;  Sum  THE  PRESSURE  FIELDS  OF  THE  SO»/RCeS  COHERENTLY  TO  FIND 
C  TmE  magnitude 

PS(J'  *  CABS<P1  ♦  P2> 
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PSL<J)  »  :o  •ALOClO<PS< J>/l  0£-06» 

C  WRITEI4  *»  PSL(J) 

!F«PLAG3  EQ  0  00  TO  30 

d  CONTINUE 

IF^FLAC  Eo  I  0>  CALL  01*PL TiNPTS.  ANC.  PSL.  1 0) 

C  IF^FLAC?  EO  I  0»  CALL  02PL T «NPTS> I . ANC. PS.  1 1 > 

C  LOCATE  OF  the  riAXN  OEAn  AXIS 
30  VMAX I  •  PS< I > 

IF<FlA03  eo  0  0>  CO  TO  10 

ANCMAX  *00 

DO  7  !■?  NPTS 

YMAX  «  AHAXl  <PS(  1  >.  YMAXn 

IF(YP.AX  CT  VriAXI)  ANCHAX  •  ANC<  I) 

YMAXl  «  YNAX 

7  continue 

C  CALCULATE  TmE  SENSITIVITIES 
10  TVS  ■  ymAX1/<CD*C0*LY«LY/<C0N*C0N) » 

TvSL  ■  20  *ALOC10<TVS) 

PETUPN 

END 

C 

c«  . . * . . . 


SUOROUTINE  P2PlT<NPTS.  XX  YV.FlACO* 

c 

C  PrPLT  IS  TmE  two  DIHENSIONAL  PLOTTING  ROUTINE  USING  ARLLlD 

c 

DIMENSION  XXUCOO.  YY(  |000> 

INTEGER  FLACa. NPTS 
xniNi  •  xx<  1 ) 

XMAX I  •  XX  <  I  ) 

VMINI  •  YY<  I  ) 

YMAXl  •  YV<1) 

00  100  I»2. NP1S 

xniN  «  AMINl < KX< I ). xniNi ) 

XMINl  •  XMIN 

XMAX  •  AMAXl <XX( I ). XMAXl » 

XMAXl  •  XMAX 

YMIN  •  AMlNl<Y><n.YMINl) 

YMINI  •  YMIN 

YMAX  •  AMAXt (YY< I ). YMAXl > 

YMAXl  •  YMAX 
too  CONTINUE 

IF<YMAX  EO  YMIN)  YMIN  •  0  0 
OX  ■  <XMAX  -  XK|N)/FLOAT(NPTS> 

XORG«I  0 
YORG»l  0 

XO€L-<XMAX-XM|N)/20  0 
VOEL-TYMAX-YMIN)/!©  0 
AYLSN-3  0 
AXLEN-5  0 
XTITL-AXLEN/J  0 
YTITL-AVLEN^O  S» 

0X>(XMAX-XM1N)/AXLEN 
0Y-(YMAX-YMIN)/AYLEN 
CALL  PlTLFNIL'PLOT-) 
call  PLTOlMdl  .6  3.  l.A  ) 

CALL  PLTORC<XORC.  YORC) 

GOTO  <l  2.3.4.5.6.7.0,9.10.11.12)  FLAGO 

1  CALL  PLTAXIStO  -0  .AYLEN.90  .  YMIN.  VMAX.  YDEL 
l,!3MMAG  REFL  C0eF.13.  2.  -  1.  -  I> 

CO  TO  200 

2  CALL  PLTAXISTO  .0  .AYLCN,90  . YMIN. YMAX .  YDEL 
I  -  rrMAPERATURE  PRESSUREIPA).  22.  2.  -  I.-  D 

CO  TO  200 

3  CALL  PLTAXISTO  .0  .  AVLEN.  90  .  YMlN.  YMAX .  YDEL 
I,25MAP£R  PRESTDO  RE  1XI0-6PA>  25. 2.  -  1.-  1) 

CO  TO  175 

4  call  PlTAXI3<0  0  AYLEN.90  YMIN.  YMAX.  YDEL 

1  26MAPER  PART  vEl < DO  RF  I  M/S»  26  2  -  1.-  I) 

CO  T->  )  75 

5  CALL  PcTAXJS<0  0  AYLCN  90  YMIN. YMAX  YDEL 

I  roNACOvs  pwr<DO  re  i  watd  23  2.  -  1.-  d 

CO  TO  175 

6  call  Pt»AXlS»0  0  AVltN  90  VMJN  YMAX  VDlL 

I  CJMPOWER  EFFICI€NCv»00  RE  I >  25  2  -  1  -  1 > 

CO  TO  l!0 

0  AVcEN  90  VMIN  YMAX  VOEL 

IPA/ <vOLT*vOlT)  )  20,2  -  1,-  |> 


call  Pk.TAXlS*0 
1  28M0vrs  (po  Re 
CO  TO  175 

call  PLTAX|b‘0 
1  24hOC^FR  «a«ps 
CC*  TO  150 
call  PlTAX|S«0 
1  leHAMC  MAX  RESP 
CO  TO  2C0 
CALL  PLTAXibtO  0 
1  14MPRCS  itvtLiDB 
GO  TO  250 
call  Pi  tax  IS(0  0 


0  AYlEN.  90 
' VOt  T •VOt  T  J  I 


0  AYLEN  90 
DeC>  18  2 


YMIN.  YM<'X  YDEl 


YMIN. YMAX  YDEl 
I  -  1  » 


AVtEN  90 

14  4  -  I 


VMIN  YMAX  Ypt  I 
~  I  ) 


AYLEN  »MIN  YMAX  YDEL 
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I.  19HPRES  AMPLITUDE  <PA),19,  A.-  1.-  1) 

CO  TO  250 

12  CALL  PLTAXISIO  .0  .  AYLEN.  90  ,  YMIN,  YMAX.  YDEL 
1.24HQCTFI  (AMPS/(VOLT*VOLT)  ).  24,  4.  -  1,-1) 

150  CALL  PLTAXIS(0  ,0  ,  AXLEN,  0  , XMIN, XMAX,  XDEL,  lOHINPUT  FREQ, 

1-10,  4,  -  1,  -  1  > 

GO  TO  300 

175  CALL  PLTAXISIO  ,0  ,AXLeN,  0  ,  XMIN,  XMAX,  XDEL,  14HL0G  INPUT  FREQ, 
1-14,  4,  -  1,  -  1  ) 

CO  TO  300 

200  CALL  PLTAXISIO  ,0  , AXLEN,  0  , XMIN,  XMAX.  XDEL,  17HAC0USTIC  FREQIHZ), 
1-17.  4,  -  1,  -  1  ) 

CO  TO  300 

250  CALL  PLTAXISCO  .0  . AXLEN.  0  , XMIN, XMAX,  XDEL 
l.llHANGEL  (DEG), -1  1, 4.  -  1,-  1) 

CO  TO  300 

300  CALL  PLTAXISCO  .  AYLEN,  AXLEN,  0  ,  XMIN,  XMAX,  XDEL,  LABX,  0,  0,  -  1,-1) 
CALL  PLTAXISC  AXLEN,  0  .  AYLEN,  90  ,  YMIN,  YMAX,  YDEL,  LADY,  0,  0,  1,  1) 

CALL  PLTDATAC  XX.  YY,  NPTS,  0.  0.  XMIN.  DX,  YMIN,  DY.  0  08) 

CALL  PLTLINE( XTITL. YTITL.  -  1) 

CALL  PLTEND( 11  0.  8  5) 

RETURN 

END 
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